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I WPROWADZENIE 


W literaturze światowej dotyczącej termodynamiki roztworów można 
znaleźć informacje o ponad tysiącu układów dwu- i trójskładnikowych wy- 
kazujących ograniczoną wzajemną rozpuszczalność. W większości są one 
jednak mało precyzyjne, uzyskane bez zwracania uwagi na dokładność po- 
miaru parametrów, czystość stosowanych substancji i.t.d. W rezultacie 
informacje o tym samym układzie mogą różnić się znacznie w zależności od 
źródła. Klasycznym przykładem jest układ (metanol(1) + cykloheksan(2)), 
dla którego podawana temperatura krytyczna mieszalności różni się o kil- 
kadziesiąt stopni. Danych podających np. precyzyjne temperatury krytycz- 
ne rozpuszczalności czy też składy krytyczne jest znacznie mniej. Ich 
liczba zawiera się w granicy kilkudziesięciu, natomiast danych dotycza- 
cych innych wielkości termodynamicznych, właściwości optycznych, elek- 
trycznych, akustycznych i innych danych fizykochemicznych ogranicza się 
do kilkunastu układów dwuskładnikowych, w których przeważająca liczba to 
układy wodne. 


W ostatnich latach obserwuje się wzrost zainteresowania w badaniach 
układów wykazujących niemieszalność. Wynika to z zapotrzebowania techno- 
logów (zwłaszcza dla układów wodnych ze względu na ochronę środowiska), 
rozwoju technologicznego pozwalającego mierzyć szybciej i z większą dok- 
ładnością, powstawania nowych teorii stanu ciekłego i roztworów, które 


wymagają weryfikacji doświadczalnej. 


Na tym tle badania równowag fazowych, a zwłaszcza równowag ciecz- 
ciecz-para w otoczeniu punktu krytycznego mieszalności, są rzadkością. 
Przyczyną tego jest najprawdopodobniej czasochłonność takich pomiarów. 
Układ, w którym współistnieją trzy fazy jednocześnie, ze swej natury 
długo dochodzi do stanu równowagi termodynamicznej. Stąd wynika koniecz- 
ność utrzymywania stabilnej i dokładnie zdefiniowanej temperatury przez 
dłuższe okresy czasu, co jeszcze do niedawna nie było łatwe. Rozwój 
elektroniki pozwolił na budowę kontrolerów temperatury, termometrów cyf- 


rowych i.t.d., co ułatwia prowadzenie długich doświadczeń. 


Punkt krytyczny, јако zjawisko przej$cia рагу w równowadze z ciecza 
w gaz, został zaobserwowany już w ubiegłym wieku. (Stanley 1971, Sengers 
i Levelt-Sengers, 1978) Istnienie punktu krytycznego mieszalności w dwu- 
składnikowych układach ciekłych zostało stwierdzone doświadczalnie przez 
Alexeva już w 1886 roku. Van der Waals w swojej pracy doktorskiej zało- 
żył ciągłość między stanem ciekłym i gazowym i przedstawił równanie sta- 
nu, którego forma matematyczna przewiduje występowanie punktu krytyczne- 
go czyli spełnienie warunku dP/dV=0. Odkrycie przejścia lambda w ciele 
stałym i podobnego mu przy przejściu helu w stan nadciekły z jednej 
strony, a ukazanie się prac teoretycznych Osangera nad modelem Isinga - 
z drugiej - wzmogły zainteresowanie punktami krytycznymi. Ponadto 
zwrócono uwagę na analogie pomiędzy przejściami fazowymi i ich punktami 
krytycznymi występującymi w ciałach stałych i w cieczach. 


Dotychczasowy kierunek badań równowag fazowych w aspekcie osobli- 
wości punktu krytycznego, zarówno eksperymentalny jak i teoretyczny, 
koncentrował się na układach jednoskładnikowych oraz dwuskładnikowych w 
pobliżu ich lini pseudokrytycznych. Lata ostatnie przyniosły kilka prac 
dotyczących zachowań się układów dwuskładnikowych w pobliżu krytycznej 
temperatury mieszalności, lecz w dalszym ciągu brak jest dokładnych da- 
nych o równowagach ciecz-ciecz-para z tego rejonu, które byłyby w stanie 
zweryfikować przewidywania teoretyczne. Co prawda, zasadność weryfikacji 
osiągnięć teoretycznych metodami doświadczalnymi jest poddawana w wat- 
pliwość (Levelt-Sengers, 1983). Niemniej, chęć empirycznego sprawdzenia 
przewidywań teoretycznych wymusza stosowanie coraz bardziej wyrafinowa- 


nych metod. 


II CEL PRACY 


Celem ргасу byta próba stworzenia konsystentnego opisu termodyna- 
micznego zjawiska rozpuszczalności ciekłych układów dwuskładnikowych w 
bezpośrednim otoczeniu górnego punktu krytycznego mieszalności. Wobec 
braku w literaturze dokładnych danych o równowadze ciecz-ciecz-para zbu- 
dowano stanowisko do pomiarów prężności par metodą statyczną nad roztwo- 
rem tworzącym dwie fazy ciekłe oraz zbadano wybrany układ. Na tej pods- 
tawie, oraz wykonując badania dodatkowe dokonano próby konsystentnego 


opisu układu. 


III PODSTAWY TERMODYNAMIKI ROZTWOROW О OGRANICZONEJ MIESZALNOSCI 


Szczegółowy opis termodynamiki roztworów można znaleźć w wielu mo- 
nografiach, podręcznikach akademickich i specjalistycznych (na przykład: 
McGlashan, 1979; Rowlinson i Swinton, 1982; Prausnitz i wsp. 1986; Mala- 
nowski i Anderko, 1991). Pojęcia stosowane w tej pracy, takie jak molowe 
wielkości cząstkowe, funkcje nadmiarowe i mieszania, współczynniki akty- 
wności i.t.d. nie będą dyskutowane. Ich definicje i zależności można 
znaleźć, w formie bardzo skondensowanej, w wielu pracach doktorskich do- 


tyczących termodynamiki roztworów (ostatnio Antosik, 1991). 


W rozdziale tym przedstawione zostały w sposób skrótowy termodyna- 
mikę (w ujęciu klasycznym) roztworów dwuskładnikowych posiadających na 
tyle duże odchylenia od doskonałości, że prowadzą one do rozpadu homoge- 


nicznego roztworu na dwie fazy. 


III.1. KLASYFIKACJA UKŁADÓW NIEMIESZAJĄCYCH SIĘ. 


Ze zjawiskiem niemieszania się dwóch cieczy spotykamy się bardzo 
często. W przybliżeniu zerowym (Eckert, 1976) zjawisko to tłumaczy się 
tym, że cząsteczki tego samego rodzaju oddziałują ze sobą mocniej niż z 
cząsteczkami drugiego rodzaju, co powoduje, że w pewnych zakresach рага- 


metrów P,T,x istnieją dwie fazy ciekłe. 


Na rysunku III-1 przedstawione są schematycznie podstawowe typy 
rówmowag fazowych ciecz-ciecz (tzw. luki mieszalności) w układach dwu- 
składnikowych znajdujące sie pod stałym ciśnieniem. Najczęściej spotyka- 
nym typem jest typ I. Krzywa współistnienia zamknięta jest, w miarę 
wzrostu temperatury, górną temperaturą krytyczną mieszalności (ang. 
UCST), natomiast obniżanie temperatury powoduje wypadanie fazy stałej w 
postaci roztworu stałego, eutektyku czy też składnika czystego. Mecha- 
nizm powstawania drugiej fazy ciekłej jest tłumaczony tym, że poniżej 


UCST oddziaływania pomiędzy cząsteczkami tego samego rodzaju przewa- 


zaja nad czynnikiem entropowym przeciwdziatajacemu tym oddziatywaniom. 


T I dwie П 
UCST fazy 
dwie 
fazy LCST 


x x 


Rys III-1. Podstawowe typy równowag ciecz-ciecz. (Schematycznie) 


Typ II, spotykany bardzo czesto w uktadach wodnych, posiada dolny 
punkt krytyczny (ang. LCST). Hirschfelder, Stevenson i Eyring (1937) 
pierwsi postawili hipoteze, Ze ten typ równowagi jest wynikiem istnienia 
silnych oddziaływań specyficznych np. typu wiązania wodorowego. W niz- 
szych temperaturach cząsteczki różnego rodzaju silnie oddziałują ze sobą 
właśnie poprzez te wiązania powodując wymieszanie cieczy. W miarę wzros- 
tu temperatury w pierwszej kolejności zrywane są te wiązania i znaczenia 
nabierają oddziaływania van der Waalsa molekuł podobnych. Dalszy wzrost 
temperatury powoduje nałożenie się równowagi ciecz-ciecz na linię kry- 
tyczną przejścia ciecz-gaz. Według Rowlinsona (1969) układy takie należą 
do III typu jego klasyfikacji wysokociśnieniowych równowag fazowych. 


Układ (nikotyna(1) + woda(2)) jest klasycznym przykładem równowagi 
typu III. Oba mechanizmy podane wyżej występują w tym przypadku. 
Najczęściej skład krytyczny UCST jest przesunięty względem LCST jak jest 
uwidocznione na rysunku. Wynika to najprawdopodobniej z różnicy między 
oboma mechanizmami rozwarstwienia. Z hipotezy Hirschfeldera wynika 
również wniosek, że pewnego rodzaju cząsteczki oddziałują ze sobą przy 
pomocy przynajmniej dwóch rodzajów sił, z których jedna jest silnie spe- 


cyficzna i skierowana w wybranym kierunku względem geometri cząsteczki. 


Typ IV spotykany jest pod wysokimi ciśnieniami. Przegląd tego typu 


rownowag przedstawiony jest w pracach Schneidera (1972, 1978). Zachowa- 
nie takie wykazuje np. uktad (butanol-2(1) + мода(2)) (Moriyoshi i wsp, 
1975) czy tez roztwory polimerów w niektórych niskoczasteczkowych 
weglowodorach ((polistyren(1) + cykloheksan(2)), Freeman i Rowlinson, 
1960) сту rozpuszczalnikach ((polistyren(1) + aceton(2)), Szydiowski i 
Van Hook, 1991). 


W przegladowej pracy Serensena i wsp. (1979) wystepuje, jako osobna 
grupa, tzw. luka wstegowa, dla której obszar niemieszalnošci rozciaga 
sie od równowagi ciecz-ciato state az до linii krytycznej przejšcia 
ciecz-gaz. W uktadach tych nie wystepuje punkt krytyczny mieszalnosci i 
dlatego ten typ luki nie jest dyskutowany w niniejszej pracy. 


III-2 TERMODYNAMIKA UKLADÓW Z LUKA MIESZALNOSCI 


Rozpatrując układ wykazujący ograniczoną mieszalność м warunkach 
izotermiczno-izobarycznych przebieg funkcji mieszania Gibbsa g w funk- 
cji składu ma analogiczny kształt jak funkcja Helmholtza dla układu jed- 
noskładnikowego w funkcji objętości (np. Rowlinson i Swinton, 1982), W 
obu przypadkach krzywa jest wklęsła. Należy podkreślić, że funkcja Gib- 
bsa (G) układu wieloskładnikowego w przestrzeni parametrów P i T jest 
funkcją wypukłą, analogicznie јак dla  utadu  jednoskładnikowego. 
Wklęsłość funkcji jest warunkiem stabilności materiałowej, która obok 
stabilności termicznej i mechanicznej warunkuje obecność fazy homege- 
nicznej. Na rysunku III-2 przedstawiony jest (schematycznie) przebieg 
funkcji G" dla krytycznej temperatury rozpuszczalno$ci To oraz dla tem- 
peratur powyżej i poniżej To Dla temperatury T«T. linia styczna wyzna- 
cza sktady równowagowe wspólistniejacych faz. Z rozwazah termodynamiki 
wiemy, że fazy znajdują się w równowadze gdy potencjały chemiczne 
składników w poszczególnych fazach są sobie równe. Jeśli dla roztworu 


dwuskładnikowego potencjały te przedstawimy w postaci: 


u 


Q7 8 7 х,(98/9х,), т III-1 


u 


Q7 8 + x, (8g/8x,) III-2 


2 РАТ 
gdzie g=G"/RT, to po odpowiednim przekształceniu 1 oznaczeniu 
poszczególnych faz ' i '' oraz odjęciu stronami powstałych równań 


uzysku jemy: 
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ЕЗ А - [ дв | III-3 


дх дх 
2 ^P,T 2 ”P,T 


Jest to termodynamiczne przedstawienie warunku równowagi fazowej. Można 
wykazać, że warunek ten jest równoważny równości aktywności składnika i 


w każdej fazie: 

x 7, = x, 7, III-4 
gdzie 7, jest współczynnikiem aktywności składnika i. Część krzywej za- 
warta między punktami a i b przedstawia zakres stężeń, dla których homo- 


geniczna faza jest bezwarunkowo niestabilna, natomiast zawarta między x' 


1 a oraz b i х"" zakres stężeń metastabilnych. 


Rys. 111-2. Schematyczny przebieg swobodnej entalpii mieszania dla układu 
dwuskładnikowego wykazującego niemieszalność. 


Aby spełnione były warunki stabilności materiałowej układu homoge- 
nicznego druga pochodna funkcji mieszania Gibbsa względem składu musi 
być dodatnia. 


ас ах > 0 III-5 


W miarę obniżania temperatury roztworu o składzie krytycznym, jeśli 
układ wykazuje tenedencję do rozdzielenia się na dwie fazy, wartość tej 
pochodnej będzie się zmniejszać, aż w punkcie krytycznym osiągnie zero. 
Dalsze obniżanie temperatury powoduje, że dla pewnego zakresu stężeń, 
krzywa с" staje sie wypukła i następuje rozpad na dwie fazy. Dla roztwo- 


rów rzéczywistych energia mieszania Gibbsa jest sumą dwóch udziałów: 
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energii mieszania roztworów idealnych i nadmiarowej energii Gibbsa (a), 
Układy z odpowiednio dużymi, dodatnimi odchyleniami od doskonałości będą 
wykazywały tendencję do rozpadu na dwie fazy. 


Tak więc warunkiem istnienia dwóch faz ciekłych jest aby: 
8%G'*/8x* + 8%G7/8x* = RT/(x,x,) + 8%G'/8x* < 0 III-6 
czyli 
BGO < - RI/(xix) < 0 III-7 


Jak powiedziano, w punkcie krytycznym 


аб! ах“ + 87G'/8x* = 0 111-8 
Najprostrza forma reprezentacji nadmiarowej energii Gibbsa jest wy- 
razenie uzyskane na podstwie teorii roztworów regularnych i które moze 


być łatwo wyprowadzone z modelu siatkowego cieczy (Guggenheim, 1952). 


G'- A RT x x, III-9 
Z rozniczkowania tego równania, па podstawie warunku stanu krytycznego 
III-8 wynika, Ze w punkcie krytycznym watrošé A wynosi 2. Dla wartošci 
nizszych mamy до czynienia 2 roztworem jednofazowym, a dla wartošci 
wyższych - dwufazowym. Jeśli stała А jest zależna od temperatury, i ma- 
leje ze jej wzrostem to mamy do czynienia z górnym punktem krytycznym 
mieszalno$ci, natomiast w przypadku odwrotnym - z dolnym. Mozna sobie 
również wyobrazić temperaturowa zależność stałej A w postaci пр. parabo- 
li. W tym przypadku może przyjąć dwukrotnie wartość krytyczną i mamy do 
czynienia z luką zamkniętą. Sczegółowo  tzkie zachowanie zostało 


rozpatrzone przez Buchowskiego (1982) w podręczniku "Chemia fizyczna". 


Z rozważań powyższych wynikają dwie ważne konsekwencje. W przypadku 
górnego punktu krytycznego układ będzie posiadał dodatnią entalpię mie- 
szania i nadmiarową entropię większą niż (-R/2), układy z dolnym punktem 
krytycznym - odwrotnie - ujemną entalpię mieszania i entropie nadmiarową 


mniejszą niż (-R/2). 


Ciepło właściwe pod stałym ciśnieniem w układzie dwuskładnikowym 
jest analogiem ciepła właściwego w stałej objętości układu jednosktad- 
nikowego Na jego wielkość składają sie dwa udziały: pierwszy, związany z 
poszczególnymi indywiduami tworzącymi roztwór i jego zachowanie jest 


ciągłe przy przejściu przez punkt krytyczny, oraz drugi, związany ze 


zmianami stezeñ w kazdej fazie i roznicy czastkowych ciepel molowych. 
Ten drugi sktadnik jest rózny w uktadzie jedno- i dwufazowym. Róznica 
ta, przy założeniu niezerowej wartości czwartej pochodnej са wzgledem 


składu, wynosi: 


sle == ге 111-10 
ste) 
ax4 © 
Ciepło właściwe oznaczone  subskryptem 2 oznacza układ  dwufazowy 
natomiast 1 - jednofazowy; ls oznacza wartość wpunkcie krytycznym. Jak 
widać różnica ta ma wartość skończoną. Podobnie wygląda różnica we 
współczynniku ściśliwości izotermiczne ј: 


2 
3 | E | 
ax2 © 
тер |- = —Á 1^ III-11 
| ele 
8x 
oraz dla współczynnika rozszerzalności izobarycznej: 
3 Е || 8s ) 
2'c 2'c 
Ge wid w ш л зок = 111-12 
2 (Tre 4 
8 g 
A mi 
8x 


Z powyższych wzorów wynika, że w przypadku ciepła właściwego i 
współczynnika ściśliwości różnica jest dodatnia, tzn. wielkość w 
układzie dwufazowym jest większa niż w jednofazowym. Dla współczynnika 


rozszerzalności zależy ona od znaków drugich pochodnych w liczniku. 


Jak wiadomo, na skutek wprowadzenia dodatkowego składnika, w 
przypadku roztworów punkt krytyczny będzie zależny od ciśnienia, pod 
którym znajduje się układ. Zależność tę można określić również na 
podstawie drugich pochodnych objętości i entropii w punkcie krytycznym. 


МЕЛЕ мана 


Constant Т 


xy 


Constant T 


100 
P/kPa 


Ху, У] 


Вуз. 111-3. Przykłady heteroazeotropii pierwszego (I) i drugiego (II) ro- 
dzaju (gwiazdka oznacza parametry azeotropu). 


III-3. RÓWNOWAGA CIECZ-CIECZ-PARA 


Układy posiadające lukę mieszalności wykazują zjawisko heteroazeo- 
tropii. Związane to jest z tym, że wrzący układ nie wykazuje zmiany tem- 
peratury wrzenia do czasu zaniku jednej z faz ciekłych. Zjawisko to jest 
zewnętrznie podobne do prawdziwej azeotropii, gdy brak zmiany temperatu- 
ry wrzenia spowodowany jest jednakowym składem fazy parowej i ciekłej. 
Van Ness i Abbott (1982) sklasyfikowali dwa rodzaje układów posiadają- 


cych tę własność ponieważ obok heteroazeotropii może występować azeotro- 
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pia homogeniczna dla stezenia, w którym istnieja tylko dwie fazy - 
ciekła i parowa. Badany w tej pracy układ wykazuje heteroazeotropie dru- 
giego odzaju. Przykładem dla pierwszego jest układ {тефапо1(1) + hek- 
san(2)). Ilustruje to rysunek 111-3. 


Pnieważ prężność pary roztworu w zakresie stężeń wewnątrz luki mie- 
szalności pozostaje stała (de facto, w miarę zmiany stężenia nominalnego 
rzeczywiste stężenia w obu fazach pozostają stałe, a zmieniają sie једу- 
nie względne ich ilości aż do zaniku jednej z faz) uzyskujemy cztery wy- 
rażenia wiążące ciśnienie nad roztworem z jego składem: 


o 


х zł” s 
P = Y i,j = 1,2 III-14 
i МАГ Y 


gdzie indeks i oznacza numer składnika, a j numer fazy. Z punktu widze- 
nia obliczeń równowagowych, wg. Van Nessa i Abbotta (1982), najlepsze 
jest obliczanie ciśnienia całkowitego z ciśnień cząstkowych dla każdego 


składnika w innej fazie. 


% У P° x y P 
P = 12 12 1 21 21 2 III-15 
Piz 9 Pa (x) 


Jest to szczególnie ważne, gdy mamy do czynienia z praktycznym niemie- 
szaniem się dwóch składników. Wtedy graniczny współczynnik aktywności 
pierwszego składnika rośnie do nieskończoności, a dla drugiego jest 


praktycznie równy jedności. 
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IV TERMODYNAMIKA ROZTWORÓW W OTOCZENIU PUNKTU KRYTYCZNEGO 


Występowanie osobliwości w otoczeniu punktu krytycznego przejścia 
ciecz-gaz i w otoczeniu punktu krytycznego mieszalności w układach dwu- 
składnikowych było znane już w końcu ubiegłego wieku (np. Alexev, 1886). 
Niestosowalność się danych doświadczalnych do równania van der Waalsa 
zostało eksperymentalnie potwierdzone przez Michaels'a dopiero w 1937 
roku. Ta rozbieżność pomiędzy teorią a eksparymentem przyjęto nazywać 
nieklasycznym zachowaniem się układów w otoczeniu punktów krytycznych. 
Od tamtego czasu zbadano i zmierzono wiele układów, ukazało się wiele 
prac prezentujących zjawiska krytyczne zarówno pod względem teoretycznym 
jak i doświadczalnym. Większość tych prac jednak dotyczy układów jedno- 
składnikowych (Levelt-Sengers, 1975; Sengers i Levelt-Sengers, 1978) lub 
zjawisk dotyczących innych układów niż ciekłe (Gopal, 1981). Ostatnio 
ukazały się także przeglądowe prace dotyczące punkyu krytycznego 
równowagi ciecz-ciecz układów dwu- i trójskładnikowych (Scott, 1978; Le- 
velt-Sengers, 1983; Levelt-Sengers i wsp., 1983a). 


Każda ze współistniejących faz scharakteryzowana jest przez wiel- 
kości intensywne i ekstensywne. Lecz istnieją takie wielkości intensyw- 
ne, jak na przykład cząstkowe wielkości molowe czy gęstośc, które dla 
każdej ilości fazy są jednakowe w całym jej obrębie, lecz są różne dla 
każdej ze współistniejących faz. Aby uniknąć związanych z tym niejednoz- 
naczności Griffiths i Wheeler (1970) nazwali je "gęstościami" natomiast 
rzeczywiste wielkości intensywne "polami". Nomenklatura ta jest szeroko 
stosowana w termodynamice punktów krytycznych. Punkt krytyczny w świetle 
powyższych definicji może być określony jako taki punkt w przestrzeni 
pole-gęstość, w którym zanikają różnice pomiędzy gęstościami poszcze- 


gólnych faz w równowadze. 


Innym pojęciem, wprowadzonym również przez Griffiths'a i Wheelera 
(1970), jest pojęcie słabej i silnej dywergenc ji" . О danej własności 


+ Pojęcie dywergencja jest anglikanizmem. Angielskie słowo diverge ozna- 


cza rozbiega się. Kojarzy sie to z brakiem zbieżności. Stosowane w га- 
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mówimy, że posiada silną dywergencje jeśli jej wartość zmienia się tak, 
jak odwrotność różnicy temperatury |Т - T |. 


W dalszych czešciach tego rozdziatu zostana porównane zachowania 
się układów jedno i dwusktadnikowych w otoczeniu ich punktów krytycz- 
nych, rozpatrywanych pod względem termodynamiki fenomenologicznej i nie- 


klasycznej, analogie i różnice pomiędzy nimi. 


IV.1. PUNKT KRYTYCZNY W UKŁADACH JEDNOSKŁADNIKIWYCH 


Rozważając układ jednoskładnikowy w stanie gazowym, podczas obni- 
żania temperatury do bliskiej krytycznej, zbliża sie on do stanu margi- 
nalnej stabilności. Zaburzenie parametrów stanu układu powoduje odpo- 
wiedź w postaci zaburzeń wielkości pochodnych jego parametrów. Zakła- 
dając, że funkcje stanu oraz ich pochodne są funkcjami analitycznymi 


można je rozwinąć w szeregi Taylora w otoczeniu punktu krytycznego. 


Powszechnie znane równanie stanu van der Waalsa jest klasycznym 
równaniem opartym na teorii pola uśrednionego, tzn. pod uwagę brane są 
tylko oddziaływania krótkiego zasięgu. Ponieważ to równanie jest anali- 
tyczne w punkcie krytycznym pod względem temperatury i objętości znako- 
micie nadaje się do demonstrowania dróg dojścia do stanu krytycznego i 
wpływu wspomnianych zaburzeń. Warunki stanu krytycznego, wyrażone jako 
pochodne energii swobodnej względem objętości wynoszą: 


2 3 2 
| 2А = 一 ӘР | = 0; | 一 = 一 [ 2r = IV-1 
av T OV ^T av T av T 


Równanie Gibbsa-Duhema dla jednej fazy po przeskalowaniu objętością 
i transformacj Legendra (Beegle i wsp. 1974) przybiera postaé: 


d(A/V) = -sdT + ydp IV-2 
Wprowadzenie zmiennych zredukowanych wzgledem parametrów krytycznych: 


A= AMVP Те T/T те ри P'-P/P p= up /P IV-3 
= c' à c' Р 3 P p я c H x: WP. © 


chunku wektorowym pojecie rozbiezno$ci nie wydaje sie adekwatne do zja- 


wiska. 
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pozwala rozwinąć w szereg Taylora zredukowaną energie Helmholtza, ciś- 
nienie i potencjał chemiczny. 


LN 1 * m * > 
X =) L -arar X (Ap )"(AT ) IV-4 


gdzie za x mogą być WE А, p lub Р" па X in га m- 
ta pochodna wzgledem p i n-ta pochodna wzgledem T . A oznacza T - 1 lub 
p 1. Ponieważ џи = (дА идр") i P°= up- A' to współczynniki 
rozwinięcia są między sobą współzależne: 


и =A m > 0 


mn ат 
= А + (m = ПА ; m > O IV-5 
,n mn 
"3 Pa T" Lan RE. ма 


Warunki stanu krytycznego (IV-1) narzucają współczynniki: 


A=A=0 м u =0; P =P =0 IV-6 
20 30 10 20 10 20 
natomiast warunki mechanicznej stabilności: 
Pu =А > Ü i Р 20 IV-7 
зо “зо 40 11 
a warunki stabilności termicznej: 
A «0 IV-8 


ог 
Z powyższych warunków wynikają zachowania się wielkości zredukowanych 
przy pewnych ustalonych parametrach: 


Izoterma krytyczna (AT °)=0 
Р-Р =и- р" = P. (Ар )°/6 + IV-9 
е К Be F, sa ФА 
Izochora krytyczna (Ap) = 0, vd = 0 
* * 
1/K_ =P АТ +... IV-10 
T 11 
* * 
СУТ = =Å. ... IV-11 
у 2 
Wzdłuż krzywej współistnienia 


Ap = +B јат |^? 


+ B IAT | Я IV-12 
(Ap + Ap )/2 = В ЈАТ | + IV-13 
AL АР, T "а z 


gdzie współczynniki В,, В... за kombinacjami współczynników rozwinięcia 
Taylora. Др oznacza róznice gesto$ci fazy parowej (indeks V) i ciektej 
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(indeks L). 
W punkcie krytycznym wyznacznik 


2 2 
au ш AE 
as V 8sav C, K. w 
= = D(S, V) IV-14 
2 2 Tap Cp 
^ [El B = 
950V av 
S 
2 T 
D(S, V) = = KC > 0 IV-15 
asav T V 


gdzie: T - temperatura, ар - współczynnik rozszerzalności termicznej pod 
stałym ciśnieniem, C ciepło właściwe w stałej objętości, Cp - ciepło 
właściwe pod stałym ciśnieniem, K. - współczynnik ściśliwości izoter- 
micznej, V - objętość. Dwie wielkości w powyższym wyznaczniku mają 
szczególne znaczenie w termodynamice punktu krytycznego i zostaną 


omówione w dalszych rozdziałach. 


IV.2. PUNKT KRYTYCZNY W UKŁADACH DWU- I WIELOSKŁADNIKOWYCH 


Układy dwu- i wieloskładnikowe wykazujące tendencję do rozwarstwie- 
nia mogą posiadać conajmniej dwa punkty krytyczne (dla mieszaniany o us- 
talonym składzie krytycznym), z których jeden jest punktem krytycznym 
przejścia gaz-ciecz, analogiczny do tego przejścia w układzie jedno- 
składnikowym. Pozostałe natomiast są punktami krytycznymi mieszalności 
dwóch faz ciekłych. Zachowanie się funkcji termodynamicznych w otoczeniu 
tego punktu jest podobne do wspomnianego przejścia ciecz-gaz i podlegają 
takim samym prawom. Warunkami stanu krytycznego w przypadku mieszanin 


są: 


2 3 2 
(=) = 一 JJ = 0; (=) - - [58] = о IV-16 
дх Т,Р Ox ^T,P дх T,,P 8x T,P 
Zastosowanie entalpii swobodnej jest bardziej odpowiednie gdyż pozosta- 


wia w dalszym ciągu zależność od jednej gęstości. 


Powyższe rozważania odnoszą sie do procesu izobaryczno - izoter- 
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micznego. М procesie ogólnym w przestrzeni P,T,x konsekwencją wprowadze- 
nia dodatkowego składnika jest przekształcenia się krzywej współistnie- 
nia w powierzchnię a punktu krytycznego w linię” . Ilustruje to rysunek 
IV-1. 


Rys. IV-1. Diagram fazowy układu dwuskładnikowego w przestrzeni P,T,x. 


W układach dwuskładnikowych w okolicy punktu krytycznego również musi 
być spełniony warunek zerowania się wyznacznika macierzy drugich pochod- 
nych entalpii swobodnej Gibbsa względem jej parametrów charakterystycz- 
nych T i -P, tzw. "macierzy podatności" (porównaj IV-14 i IV-15). 


8%G a<G C 
ue. = Va, 
ШЕ 298 ӘТӘР |, T 4 
, - = D(T,P) IV-17 
ас ас Vap VK, 
aTaP Әр? 
x T, x 
2 2 2 2 
D(T,P) = 25 一 » | ZĘ = VK C/T > 0 IV-18 
8T 8P 8T8P 
P,x Ге 


* Podobnie dzieje sie z punktem krytycznym przej$cia ciecz-gaz. Linia 
krytyczna łączy punkty krytyczne czystych składników. 
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Funkcje termodynamiczne w układzie k-składnikowym określone są, јак 
wiemy, w przestrzeni k+1 parametrów. Stąd, układ dwuskładnikowy posiada 
trzy parametry niezależne. Parametry te są bądź gestościami, bądź polami 
stosując nomenklaturę Griffiths'a i Wheelera (1970). W tablicy IV-1 
przedstawione jest za Scott'em (1978) zestawienie funkcji dla układów 
jedno i dwuskładnikowych. 


Tablica IV-1 


Podział funkcji termodynamicznych ze względu na stosowanie zmiennych 
niezależnych. 


pe Układ jednoskładnikowy Układ dwuskładnikowy 


A(T,V,x); H(S,P,x) 
P(T,V,x); T(S,P,x) 
S(T,V,x); V(S,P,x) 
C, G, V, x); K (S, P,x) 


A(T,V); H(S,P) С(Т,Р,х) 

P(T,V); T(S, P) A(T, P, x) 

S(T,V); V(S, P) S(T,P,x); V(T,P,x) 

c (T, V); K_(S,P) C, (T, P, x); K (T, P, x); op (T, P, x) 


G(T,P)=u(T,P) u. (T, P, A) 
S(T,P); V(T,P) SU,P,A); V(T,P,A); x(T, P, A) 
COPS R. CP) ќ СЪР) © (T; P.A); K,(T,P, A); а (T, P, A) 


Q(T, P, A) 


IV.3. NIEKLASYCZNE ZACHOWANIE W OTOCZENIU PUNKTU KRYTYCZNEGO 


Powszechnie znane jest tzw. nieklasyczne zachowanie się funkcji 
termodynamicznych w otoczeniu punktu krytycznego. Termin "nieklasyczne 
zachowanie" oznacza, że rzeczywisty przebieg pewnych funkcji termodyna- 
micznych nie jest zgodny z obliczonym czy przewidzianym przez klasyczne 
równanie stanu typu Van der Waalsa. Osobliwości tego punktu nie ograni- 
czają się do wielkości termodynamicznych w potocznym tego słowa znacze- 
niu. Tak więc obserwuje się również silny wzrost rozpraszania światła, 
wzrost lepkości, zwiększone pochłanianie dźwięku czy też zmiany 
własności elektrycznych. Ponadto, jak już wspomniano, nie dotyczy to 
tylko układów ciekłych, lecz podobne zachowanie obserwuje się w wielu 
innych układach, dalekich fizycznie od stanu ciekłego, jak np. przejście 


ferromagnetyk - paramagnetyk w temperaturze Curie. 
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T-T, 


š 


Rys.IV-2. Ilustracja pojęcia długości korelacji. 


Według współczesnej termodynamiki (Levelt-Sengers 1983) anomalie 
związane z punktem krytycznym wynikają ze wzrostu długości korelacji, co 
powoduje wzrost fluktuacji w układzie. Na rysunku IV-2 przedstawiony 
jest, na przykładzie rozkładu cząsteczek w prostej sieci, wzrost 
długości korelacji £ w miarę dochodzenia do temperatury krytycznej. W 
myśl teorii, funkcja charakteryzująca anomalie w punkcie krytycznym jest 
zależna od jednej zmiennej, która jest kombinacją temperatury i gęstości 
(składu). Takie zredukowanie ilości zmiennych nosi nazwę skalowania. Ро- 
nadto, ponieważ długości korelacji są dużo dłuższe niż zasięg sił 
międzycząsteczkowych, zachowanie się układu nie zależy od budowy 
cząsteczek wchodzących w jego skład. Zjawisko ma więc charakter uniwer- 
salny i zależy tylko od tego, do której klasy uniwersalności należy 
układ. Na przykład punkt krytyczny mieszalności i w przejściu fazowym 


ciecz-gaz zachowują się podobnie gdyż należą do tej samej klasy, nato- 
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Przejscie fazowe I rzedu 


a) b) 6) 


База 2 $ У 


| 
| 
| 
| 
! 
| 
! 


Przejscie fazowe П rzedu 


I 
I 
I 
I 
1 
I 
I 
' 
I 
1 
| 
| 
| 
| 
T, T T T T 


Rys.IV-3. Zachowanie się energii Gibbsa i jej pochodnych w temperaturze 
przejścia dla przejść fazowych pierwszego i drugiego rodzaju. 


miast przejścia fazowe w ciekłych kryształach mają inne charakterystyki. 


Dzisiaj obowiązuje teza, że układy należące do tej samej klasy uni- 
wersalności wyrażonej m.in. takim samym parametrem uporządkowania сћага- 


kteryzują się podobnym przebiegiem funkcji termodynamicznych czy też 
diagramów fazowych. 


Wszystkie wymienione anomalie wskazują na to, że punkt krytyczny 
mieszalności w układach ciekłych jest przejściem drugiego rodzaju. Ilu- 
struje to rys. IV-3. Rozpatrując pochodne entalpii swobodnej w funkcji 
temperatury, w przypadku przemiany pierwszego rodzaju nieciągłość obser- 
wowana jest już dla pierwszych pochodnych funkcji G (rys.IV-3 b i c). 
Dla przemiany drugiego rodzaju natomiast otrzymujemy zamiast skoku en- 
tropii - załamanie funkcji (rys. IV-3. b), a nieciągłość występuje w cie- 
ple właściwym (rys. IV-3 c), które jest drugą pochodną. Inną, ciekawą 
własnością punktu krytycznego jest to, że przejście zachodzi, w przeci- 


wieństwie do przejść pierwszego rodzaju, w wąskim zakresie temperatur 
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lecz za to w całej objętości układu. 


IV.3.1. Prawa wykładnicze 


Przedstawione w punkcie IV-2 równania od ТУ-9 do IV-13 definiują 
tzw. prawa wykładnicze (ang. power laws) dla klasyczneg punktu krytycz- 
nego. Wykładniki potęgowe zmiennych АТ" 1 др" oznaczane Sa literami а1- 
fabetu greckiego. Wymienione równania przyjma wiec postaé ogólna: 


Pe неи s P Соја IV-19 
е РОК = t'a 7 РЕН 
* * y 
1К = Р (AT? +... IV-20 
T 11 
* + * с 
СЛ” = -A (AT ) ... IV-21 
L 02 
Ap = +в, јат |? AE IV-22 
(др, + Ар") /2 = En E rcs IV-23 


Zgodnie z teoria uniwersalno$ci zjawisk krytycznych równania te mo- 
ва być przeniesione do opisu punktu krytycznego mieszalności. Wymaga to 
jednak zastąpienia pewnych wielkości charakteryzujących układ jedno- 
składnikowy przez ich odpowiedniki charakterystyczne dla układu wielo- 
składnikowego. Tak więc ciśnienie zastępowane jest przez różnicę poten- 
cjałów chemicznych składników, a gęstości faz, spełniające rolę рагате- 
tru uporządkowania, zastępowane są przez stężenia jednego ze składników 


w poszczególnych fazach. 


W przypadku układów dwuskładnikowych sytuacja jest bardziej skom- 
plikowana. Jak wynika z tablicy IV-1 poszczególne funkcje zależne są od 
trzech parametrów i ich określenie w funkcji jednego z nich wymaga u- 
trzymania stałości dwóch pozostałych. Konsekwencją tego jest różne га- 
chowanie się funkcji, a zatem różne wykładniki w prawach wykładniczych, 
w zależności od wzajemnej relacji między gęstościami a polami jako jej 


parametrami. 
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Tablica IV-2 
Porównanie niektórych wielkości termodynamicznych układu jedno i dwu- 


składnikowego 


jednosktadnikowy dwuskładnikowy ciekły 
P = const. 


Równanie fundamentalne |P(T, и) u, (T, A) 
Gęstość р = (8P/8u), х = "(ди Ида). p 
Krzywa współistnienia |р(Т) x(T) 
и(р, T) A(x, T) 

DRZE 2 _ аа 2 
K. = (д°Р/ди ds s (8 u, / aA РР 


Równanie stanu 
Wsp.ściśliwości 


Ciepło właściwe CY = Т(8 Р/ӘЇ ) | с = -T(8 и 781 j. 


РТ ,T 


Dla układu dwusktadnikowego, przykładowo dla niektórych własności, 
prawa wykładnicze przyjmują więc postać: 


k=l eh de y ow. T=T,P=P IV-24 
c с с 
= 
gdzie A = H м; Ах = |x = хо 
* -7 
К = С (AT)? +... PaP, x = x IV-25 
T с с 
* x 
С =A (АТ )©... рар. Хах IV-26 
P C c 
ЕЈ 2 “8 
x -x = B, [AT ЕА" Р = const. IV-27 
(х' СЕ B par" p + sa: P = const IV-28 


Współczynniki proporcjonalno$ci oznaczone literami noszą nazwę amplitud 
i są charakterystyczne dla układu. Wartości ich nie są jednak dowolne 
lecz są powiązane pomiędzy sobą. 


IV.3.2. Drogi dojścia do punktu krytycznego i wykładniki krytyczne 


Według Griffiths'a i Wheelera (1970), to czy dana własność będzie 
miała silną czy słabą dywergencję jest wynikiem kierunku zbliżania się 
do punktu krytycznego. Wielkości mierzone asymptotycznie równolegle do 
krzywej współistnienia będą wykazywały słabą dywergencję, natomiast pro- 
stopadle do niej - silną. 
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Jeśli punkt krytyczny obierzemy jako początek układu współrzędnych 
zdefiniowanych następująco: 
AT = T - Те IV-29 
AX = X - X. IV-30 
to badania zachowania sie układu w jego otoczeniu prowadzi sie wzdłuż 
wybranych kierunków przedstawionych na rysunku IV-4. Przedstawiony jest 
tutaj stan pożądany, który w praktyce jest niezmiernie trudno spełnić. 


Zagadnienie to jest omówione szerzej w następnym rozdziale. 


Rys. IV-4. Drogi dojscia do punktu krytycznego i związane z nimi wykładniki. 


Na rysunku VI-4 przedstawione są również wykładniki charakteryzujące dy- 
wergencje pewnych własności badanych wzdłuż wybranego kierunku. Jak wi- 
dać wykładniki « i у mogą być wyznaczone powyżej i poniżej temperatury 
krytycznej. Nie ma jednoznacznego, eksperymentalnego potwierdzenia rów- 
ności pomiędzy wykładnikami 2 i bez apostrofu. Np. «' może być wyznaczo- 
ny z pomiaru symetralnej krzywej współistnienia, natomiast с - na pod- 


stawie pomiarów ciepła właściwego. 


Wykładnik krytyczny powinien być obliczany w granicy gdy parametr 
zredukowany e zdąża do zera (Stanley, 1971; Scott, 1978). 
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айз. „4 = IV-31 
gdzie: Ф jest poszukiwanym wykładnikiem, y'= |? = Ж м zmiana wiel- 
kości sprzężonej, е - parametrem zredukowanym = = |E - E, |/Е Pod wiel- 
kość E może być podstawione stężenie, temperatura czy ciśnienie. Należy 
podkreślić, że jest to wartość graniczna i może silnie zależeć od odle- 
głości od punktu krytycznego mierzonej parametrem e. W przypadku tempe- 
ratury jako parametru e najmniejsza jego wartość osiągnięta eksperymen- 
talnie wynosi 10 со odpowiada róznicy pomiedzy temperatura pomiaru а 
Te rzedu kilku milikelwinów. 


М tablicy IV-3 przedstawione są wartości wykładników związanych z 
własnościami termodynamicznymi uzyskanymi eksperymentalnie, wynikającymi 
z modelu van der Waalsa oraz obliczonymi metodą renormalizacji grup z 


modelu Isinga. 


Tablica IV-3 
Wartości niektórych wykładników krytycznych. 


Wykł. | Własność Eksperym. Model Model 
у.а. М. Isinga 


Ciepło właściwe 


Krzywa współistnienia 
Współczynnik ściśliwości 
Potencjał chemiczny 


Z warunków stabilności oraz rozważań termodynamiki statystycznej 
wynika, że poszczególne wykładniki są ze sobą powiązane słabymi nierów- 
nościami. Zgodnie natomiast z teorią klasyczną odpowiednie wykładniki są 
sobie równe bez względu na kierunek pomiaru. Ponadto, tylko dwa z nich, 
a mianowicie B i 6 są od siebie niezależne. Pozostałe wynikają wprost z 


poniższych zależności. 


« + 28 ЖУ" x2 IV-32 
& +BS+ 1) £2 IV-33 
y Bio = 1) > O IV-34 
& ty ца + 1)Z(Š = 1) € 2 IV-35 


Eksperymentalnie nie potwierdzono jednoznacznie i nie obalono powyzszych 
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Rys.IV-5. Warto$ci zmiennej skalujacej x w otoczeniu punktu krytycznego. 


nierówno$ci. Dla klasycznego modelu jednak nierówno$ci te staja sie rów- 


nościami. 


IV.3.3. Skalowane równanie stanu 


Wspomniane trudności w zachowaniu kierunków dojścia do punktu kry- 
tycznego nakazywały opracować metodę pozwalającą stosować do weryfikacji 
teorii dane uzyskane w dowolnie wybranych kierunkach. Naturalnym narze- 
dziem jest równanie stanu. Griffiths (1964) zaproponował skalowane rów- 


nanie stanu stosując redukcję (skalowanie) zmiennych niezależnych. 


W termodynamice punktu krytycznego zakłada się, że na każdą 
wielkość termodynamiczną składają się dwa udziały: część regularna, 
dająca się opisać wyrażeniami analitycznymi, oraz anomalna - wynikająca 
z istoty samego punktu krytycznego. Np. zredukowana energia Helmholtza 


może być przedstawiona w postaci: 
А = А(Т') + pu (рт) + А’ (Ap , AT ) IV-36 
= o РИ Ре? sing Pp, 


gdzie A jest właśnie funkcją analityczną. Term anomalny można wyrazić 


24 


јако funkcje zmiennej skalujacej x. 


* 


А’ (Ар AT) 
sing Ps в, 


(1,х) IV-37 


|Ap” |? » A 


gdzie x = AT*/|Ap' | B. W monografii Sengers i Levelt-Sengers (1978) po- 
dane są wyrażenia jako funkcje x na takie wielkości jak ciśnienie. en- 
tropia, ciepło właściwe i.t.d. 


Praktyczną, szeroko stosowaną w NBS (Vicentini-Missoni i wsp., 
1969), forma równania stanu jest: 


др” = ap” Jap" | "E, [1 + z) Ë + Е, (1 + ы ЕЎ IV-38 


2-1)728 


E = 0/01 + Е) IV-39 


D jest amplituda opisujaca zachowanie sia potencjatu chemicznego wzdluz 


izotermy krytycznej, a E, jest dobieralnym parametrem. 


Wadą tego równania jest konieczność numerycznego całkowania danych 
w celu uzyskania potencjału chemicznego przed dopasowaniem do danych 
PTV. 


IV.3.4. Parametryczne równanie stanu 


Wobec wspomnianych niedogodności modelu NBS (Lentini i Vincentini- 
Missoni, 1973), Schofield zaproponował w 1969 roku parametryczne równa- 
nie stanu, które definiuje położenie punktu w przestrzeni Р,Т,х za рото- 
cą dwóch parametrów: r określającym odległość od punktu krytycznego i 6 
podającego kierunek do 'warstwy' o stałym r. Funkcje wiążące zredukowane 
zmienne doświadczalnei parametry równania , spełniające prawa skalowa- 


nia, są wybrane następująco: 


Ат = r т(е) T(6) = 1 - b202 IV-40 
др'= rM(6) ме) = ке IV-41 
Au = rËŠgç8o) H(6) = а e(1 - e?) IV-42 


Powyzej przedstawione funkcje T, M i H zdefiniowane sa dla modelu 
liniowego (Schofield, 1969). 
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Krzyw 
współistnienia 


Rys. IV-6. Współzależność parametrów parametrycznego równania stanu w 
otoczeniu punktu krytycznego. 


Jak widać na rysunku IV-6, parametr 6 zmienia swą wartość od 1 do 
-1 wzdłuż krzywej współistnienia. Na izotermie krytycznej 6 jest ujemną 
lub dodatnią odwrotnością parametru b, który jest obliczany na podstawie 


wykładników krytycznych 8 i à. 


W tablicy IV-4 przedstawione są wyrażenia uzyskane dla modelu 
liniowego, opisujące znane wielkości termodynamiczne za pomocą 
parametrów r i 6. Należy jednak pamiętać, że opisują one tylko udział 
anomalny. Dla uzyskania pełnego opisu konieczne jest dodanie 
składowych regularnych, które podane są w pełnej formie w pracy 
Levelt-Sengers (1983). 


Parametryczne równanie stanu jest obecnie szeroko stosowane nie 
tylko do opisu własności substancji czystych w otoczeniu punktu krytycz- 
nego, lecz również, przy zastosowaniu wykładnika korekcyjnego A, wprowa- 
dzonego przez Wegnera (1972), w szerokim zakresie temperatur, ciśnień i 
gęstości. Uzyskiwana jest doskonała zgodność obliczonych wielkości z da- 
nymi eksperymentalnymi np. dla wody (Kamgar-Parsi i wsp., 1983), 
ciężkiej wody (Levelt-Sengers i wsp., 1983), tlenu i gazów szlachetnych 
(Levelt-Sengers i wsp., 1976). Tam też można znaleźć kompletny zbiór wy- 


rażeń na części regularną i anomalną parametrycznego równania stanu. 
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Tablica IV-4 
Równania modelu liniowego dla układów dwuskładnikowych ciekłych 


Potencjał -u (T, P, A) 
"du, = SdT + xdA - VdP 


Ат = [T -T_(P)1/T_(P) г(1 - b?e?) 


x -x (P) кгбе 
"У A (T) агбде(1 = ө?) 
=u акг“ “p (0) 


2(anom) 
акг' 85 (0) 
akr2 "[6"(1 = 92) = p(0)] 
akr' “а(е) + akr” "[(1-b?e?)s(8) + 02(1-02)-р(0)] 


Е 2-а 1-а от В (аА 
(ak + а k,.)r2 "р(е) + akr se. ° kr ess). 


-  B8-38-b^uy _ 
2b* (2-0) (1– о 
sa 85—28-b^a(288-1) 
| 2b? (1-a)a 
р,= (288-3)/2« 
= (2-«)p, 


B5-36 
2b?a 


s = 
2 


IV.4. NIEKLASYCZNE MODELE SIATKOWE PUNKTU KRYTYCZNEGO 


Powszechnie wiadome jest, że modele siatkowe dają diagramy fazowe 
odbiegające od uzyskiwanych doświadczalnie. Przyczyny tego leżą w 
założeniach teoretycznych modeli. Poniżej podane są skrótowo własności 
dwóch modeli opracowanych początkowo dla układów wodnych, które opisane 


są szerzej (aczkolwiek również skrótowo) w załączniku A. 


27 


0 .2 4 х 5 5 10 


Rys.IV-7. Porównanie modelu sieci podwójnej i modelu Barkera-Focka 


IV.4.1. Model sieci podwójnej; "decorated lattice model" 


Idea "sieci w sieci" była znana od dawna (patrz artykuł przeglądowy 
Wheelera, 1977). Zasadniczą zaletą modelu jest możliwość jego rozwiaza- 
nia, gdyż łatwo jest go przetransformować do modelu gazu sieciowego (mo- 
del Isinga), który jest rozwiązywalny za pomocą rozwijania w szeregi czy 
przez renormalizacje grup. (Uproszczone rozwiązanie modelu  Isinga 


przedstawione jest w załączniku A.) 


Na rysunku IV-7. pokazana jest możliwość opisu za pomocą tego 
modelu realnego układu posiadającego lukę zamkniętą. Obliczenia 
przeprowadzono dla układu (gliceryna(1) + m-toluidyna(2)) stosując różne 
liczby koordynacyjne w. W przypadku zastosowania w = 6 (wartość 
najczęściej stosowana, również w rozwiązaniu modelu Barkera-Focka (1953) 
(patrz załącznik C) autorzy uzyskali krzywą (b). Powiększenie tej liczby 
do w = 5000 daje krzywą (с). opisującą dane eksperymentalne z 
zadowalającą dokładnością. Dla porównania krzywa (а) przedstawia 
rozwiązanie za pomocą modelu Barkera-Focka dla tego układu (czytaj, dla 
takiego samego stosunku górnej i dolnej temperatury krytycznej). Autorzy 
opisali również inne układy zawierające glicerynę. Z każdym razem stoso- 
wali tę samą liczbę koordynacyjną. Stwierdzili również, że obliczone 
wartości energii oddziaływania mają wielkości rzędu siły wiązania wodo- 


rowego, co według nich świadczy o fizycznym znaczeniu zarówno liczby ko- 
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ordynacyjnej jak i energii oddziaływań. 


Mozliwo$ci modelu Ъуїу badane przez wielu autorów w zagadnieniach 
związanych z przejściami fazowymi i zagadnieniami osobliwości punktu 
krytycznego. Tak więc przejście fazowe ciecz-para i punkt krytyczny 
ciecz-gaz badali Widom (1967) i Neece (1967). Rozszerzenia na układ dwu- 
składnikowy i azeotropię krytyczną dokonał Wheeler (1972). Układy dwu i 
trójskładnikowe były szeroko badane przez Widoma i Wheelera (Widom, 
1967; Wheeler i Widom, 1968, 1970) i wykazano analogię pomiędzy punktem 
krytycznym mieszalności i tzw. plait point'em. Clark (1968) badał му- 
kładniki krytyczne natomiast Zollweg i Mulholland (1974) wykazali możli- 
wość anomali w przebiegu symetralnej krzywej współistnienia. Bartis i 
Hall (1975) uzyskali diagram fazowy, w którym występuje równowaga "gaz- 


gaz' i związany z tym zjawiskiem podwó jny punkt krytyczny. 


T а) т A с) 
x x x 


Rys. IV-8. Krzywa wspótistnienia w okolicy temperatury krytycznej uzys- 
kanych za pomocą modelu sieci podwó jnej. 


Seria prac Wheelera i Andersena (Andersen i Wheeler, 1978a,b, 1979; 
Wheeler i Andersen, 1980) pokazała możliwość konstrukcji diagramów fazo- 
wych typu III i IV pokazanych w rozdz. III; zastosowali go do opisu 
krzywej współistnienia układów wodnych wykazujących krzywą zamkniętą. 
Ponadto (Wheeler i Anderson, 1980) model przewiduje kilka ciekawych 
przebiegów krzywej współistnienia w bezpośrednim otoczeniu temperatury 


krytycznej. 


Zachowanie się krzywej współistnienia tak, jak pokazano na rys. IV- 
8, możliwe jest w przypadku gdy nowy, skalowany wykładnik krytyczny в = 
3B przyjmuje wartosci: nieco mniejsza niz (а), równa (b) i nieco wieksza 
niż 1 (c). Według autorów zachowanie takie występuje nadzwyczaj blisko 
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temperatury krytycznej i dlatego nie może być stwierdzone eksperymental- 
nie. Jeśli wykładnik В=0.325, uzyskany teoretycznie z trójwymiarowego 
modelu Isinga jest prawdziwy, to wierzchołek krzywej współistnienia po- 
winien wyglądać jak przypadek (a) na rys.IV-8. Z drugiej strony, 
wykładnik 820.335 z tablicy IV-2, określony doświadczalnie, daje Вв > 1, 
a wtedy kształt krzywej współistnienia przewidywany przez model (c) jest 


zgodny ze stwierdzonym doświadczalnie. 


IV.4.2. Model Walker'a i Vause 


Prawie równocześnie z modelem sieci podwójnej Walker i Vause (1980, 
1983) zaproponowali nowe podejście do definicji energii oddziaływania 
międzycząsteczkowego poprzez hamiltonian. 


и - s, $, ) [K + (K - К IV-43 


NI 


gdzie K, i K, sa energiami oddzialywania czasteczek, s definiuje obec- 
ność w sąsiedniej komórce molekuły tego samego (i=j) lub różnego (i*j) 

rodzaju, natomiast 6 jest deltą Kroneckera, za pomocą której definiowana 
jest orientacja oddziałujących cząsteczek. Założenia i podstawowe 


własności modelu przedstawione są w załączniku A. 


Istotną cechą, różniącą go od modelu sieci podwójnej, jest włącze- 
nie modelu Potta dla stanów q-krotnie zdegenerowanych, czyli innymi sło- 
wy, wprowadzenie możliwości orientowania cząsteczek względem siebie w d- 
wymiarowej sieci. Model jest w stanie przewidzieć w zależności od doboru 
parametrów oddziaływania rozpad fazy ciekłej dwuskładnikowej i ponowną 
homogenizację w niższej temperaturze. Podstawowa różnica pomiędzy tym 
modelem, a omawianym poprzednio jest taka, że cząsteczki traktowane jed- 
nakowo bez względu na rozmieszczenie ich w sieci. W modelu sieci podwó j- 
nej oddziaływania międzycząsteczkowe są zależne od tego, do sieci które- 
go rzędu przynależy cząsteczka. Wartość parametru q, który można 
traktować jako liczbę koordynacyjną dla oddziaływań drugiego rzędu, 
wynosiła ok. 500, czyli jest o rząd mniejsza niż w modelu omawianym po- 
przednio. Autorzy interpretują tę wartość jako odwrotność kąta sferycz- 


nego w jakim oddziałuje wiązanie wodorowe. 
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Rys.IV-9. Zachowanie się ciepła właściwego obliczone modelem Walker'a i 
Vause 


Praca Goldsteina i Walkera (1983) prezentuje możliwości modelu w 
przewidywaniu zachowania sie pewnych wielkości termodynamicznych w funk- 
cji temperatury, takich jak entropia czy ciepło właściwe. Na rys. IV-9 
przedstawiona jest zależność temperaturowa ciepła właściwego dla luki 
wykazującej tzw. talię. Widać na nim "przygotowanie się" układu do dy- 
wergencji podobnej do widocznej po prawej stronie rysunku jeśli miałby 
pojawić się dolny punkt krytyczny. Niestety, brakuje danych mogących 
zweryfikować przewidywania modelu. 


Ponieważ opis zachowania się cząsteczek jest na poziomie mikroskopo- 
wym, możliwe jest, teoretycznie uzasadnione, przewidywanie efektu pods- 
tawienia izotopowego. Model ten ma również zastosowanie do predykcji 
równowag fazowych modyfikowanych obecnością elektrolitu. Autorzy 
porównali uzyskane rezultaty w postaci wielokrotnej luki zamkniętej ze 
stwierdzonym doświadczalnie takim zachowaniem się układu (ебапо1(1) + 
woda(2) + siarczan amonu(3)) pod wysokimi ciśnieniami (rys. IV-10). 


Ponadto przedstawiono zależność ciśnieniową modelu i wykazano dos- 
konałą zgodność uzyskanych równowag fazowych ze zmierzonymi 
doświadczalnie na przykładzie układu (2-butanol(1) + woda(2)) (Rys IV- 
11). 
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Rys.IV-10. Zanikanie luki mieszalno$ci pod wpływem ciśnienia dla układu 


(etanol(1) + woda(2) + 14% коСО„(3)} 


° 
LU 


° 
° 


° 


T/UCST (200 айт.) 


Rys.IV-11. Porównanie danych eksperymentalnych i obliczonych modelem 
Walker'a-Vause dla uktadu (2-butanol(1) + woda(2)) 


W pierwotnej wersji model był rozwiazywalny dla symetrycznej równo- 
wagi fazowej (x =1⁄2). Dalsze prace nad jego rozwojem pozwolity uzyskaé 
asymetrie wynikajaca z budowy molekularnej sktadników (Goldstein, 1983) 
oraz wprowadzenie (Walker i Goldstein, 1985) do hamiltonianu dodatkowego 
oddziaływania pomiędzy podobnymi cząsteczkami interpretowanego јако wia- 
zanie wodorowe pomiędzy cząsteczkami tego samego rodzaju. Rezultat 


przedstawiony jest na rysunku IV-12. 
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Rys.IV-12. Asymetria luki mieszalno$ci uzyskana przez wprowadzenie do- 
datkowej energii oddziaływania w hamiltonianie modelu W-V. 


Na rysunku tym, poza asymetryczna luka zamknieta przedstawione sa 
różnice entropii i fluktuacji gęstości pomiędzy współistniejącymi fazami 
oznaczonymi (-) i (+). Ponadto pokazany jest przebieg symetralnej luki, 
który wykazuje, podobnie jak w innych modelach, niezgodność z zasadą 
prostoliniowej symetralnej. 


Przedstawiony w tym punkcie model płynu wydaje się najpełniej pre- 
zentować zachowanie się cząsteczek w roztworach złożonych. W porównaniu 
z modelem Wheelera i Andersena jego możliwości predykcyjne są większe, a 
jednocześnie nie narzuca on specjalnych warunków na budowę i strukturę 
sieci. Tę "komplikację" autorzy zawarli w konstrukcji hamiltonianu. 
Wszystkie modele oparte na trójwymiarowym modelu Isinga wymagają skom- 
plikowanych rachunków i nie mogą być prosto przeniesione na inne 


własności termodynamiczne. 
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V DOŚWIADCZALNE BADANIA OTOCZENIA PUNKTU KRYTYCZNEGO MIESZALNOSCI 


V.1. SPECYFIKA BADAŃ W OTOCZENIU PUNKTU KRYTYCZNEGO 


Punkt krytyczny jako obiekt badań należy do trudniejszych. Wynika 
to z wrażliwości na wszelkiego rodzaju zaburzenia. Istnieje dosyć bogata 
literatura opisująca metody badań eksperymentalnych i związanych z tym 
trudności (Scott, 1978). Ponieważ w niniejszej pracy badanym układem 
jest ciekły układ dwuskładnikowy, rozważania na ten temat są ograniczone 
właśnie do tego typu badań. Szersze ujęcie tematu Czytelnik może znaleźć 
w "Experimental Thermodynamics, Vol.III" (Levelt Sengers, 1975). W przy- 
padku układów dwuskładnikowych np. wyznaczenie składu krytycznego przys- 
parza wielu kłopotów i niepewności związanych z czystością składników, 
braku jednoznacznych kryteriów oceny czy osiągnięty punkt na osi składów 
jest punktem krytycznym (kryterium jednakowej objętości dwóch współist- 
niejacych faz jest ciągle kwestionowane, szczególnie w bezpośrednim oto- 
czeniu temperatury krytycznej), trudnością oceny wizualnej istnienia 
drugiej fazy przy intensywnej opalescencji i.t.d. Innym czynnikiem za- 
kłócającym jest pole grawitacyjne utrudniające przechodzenie substancji 
przez granicę faz i przeciwstawiające się homogenizującemu działaniu en- 
tropii. Uzyskanie stabilnej i dobrze określonej temperatury było do nie- 
dawna barierą powodującą, że uzyskiwane dawniej wyniki są obarczone zna- 


cznym błędem. 


W rozdziale tym omówione są głównie doświadczalne aspekty pomiarów 
własności układów dwuskładnikowych w otoczeniu punktu krytcznego. W za- 
łączniku F przedstawione jest zestawienie układów dwuskładnikowych bada- 
nych w przeszłości pod wieloma względami. Nie jest to zestawienie pełne, 
ale daje obraz zakresu pomiarów nie tylko termodynamicznych sensu stric- 
to, lecz także innych własności fizykochemicznych, które ulegają wpły- 
wowi zjawisk związanych z przejsćiem ze stanu jedno- do dwufazowego. Po- 
nadto pokazuje, jak niewiele układów zostało zbadanych w tym aspekcie, w 
porównaniu do mnogości układów posiadających punkt krytyczny mieszalno$- 
ci (Francis, 1963). 
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V.2. KRZYWA WSPÓLISTNIENIA, NUKLEACJA 


Podstawową wielkością, w przypadku układu dwuskładnikowego wykazu- 
jącego tendencję do rozwarstwienia fazy ciekłej, jest krzywa współist- 
nienia. W literaturze opisano wiele metod pozwalających na jej wyznacze- 
nie, lecz za najdokładniejszą uznaje się wciąż metodę ampułkową (Alexev, 
1886). Polega ona na syntetycznym utworzeniu mieszaniny z naważek czys- 
tych składników, a następnie na oznaczeniu temperatury przejścia T 
czyli temperatury zaniku jednej z faz ciekłych. Przy napełnianiu ampułek 
należy zadbać о јак najmniejszą objętość fazy parowej. Jest to 
szczególnie ważne w przypadku, gdy oba składniki różnią się znacznie lo- 
tnościami. Doskonałym przykładem zastosowania tej metody jest oznaczenie 
krzywej współistnienia układu (cykloheksan(1) + bezwodnik octowy(2)). 
Nagarajan i wsp. (1980) określili w ten sposób T, dla ponad 60 próbek 
pokrywając zakres temperatur ok. 40 stopni. Dokładność tej metody rośnie 
w miarę dochodzenia do składu krytycznego. Jednocześnie wzrost opales- 
cencji utrudnia określenie T Osiagnieciu temperatury przej$cia towa- 
rzyszy zmetnienie próbki (ang. cloud point). Pomiar tego zjawiska, dla 
próbki o znanym składzie, jest również często stosowany w tym celu (Ochi 
i wsp., 1990). Jednak, gdy układ nie jest mieszany, jest to możliwe tyl- 
ko w bezpośrednim otoczeniu punktu krytycznego, bo tylko tam krzywa spi- 
nodalna pokrywa się z krzywą binodalną. 


Krzywą współistnienia uzyskuje się również poprzez określenie skła- 
du współistniejących faz. Tego typu pomiar powinien być preferowany ze 
względu na informację o równowadze, która może być wykorzystana w kore- 
lacjach i badaniach modeli termodynamicznych stanu ciekłego. Jednak o 
dokładności decyduje wybór metody określania składu. Ponieważ obie fazy 
różnią się pewnymi właściwościami, ich pomiar dla każdej z faz pozwala 
określić ich skład . Warunkiem jest, aby zależność mierzonej właściwości 
od składu była wystarczająco duża. Tak więc metoda ta wymaga skalowania. 
Takimi właściwościami mogą być np. gęstość (Greer, 1976), współczynnik 
załamania światła (Jacobs i wsp., 1977), natężenie pasm charakterystycz- 
nych grup funkcyjnych w widmie IR (Cannon i wsp. 1970), własności elek- 
tryczne roztworu (Śliwińska-Bartkowiak, 1977), a nawet lepkość (Dega- 
Dałkowska, 1980). Innym rodzajem oznaczania składu współistniejących faz 
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jest analiza chromatograficzna metoda GLC (Fabries i wsp., 1977). Wymie- 
nione sposoby są jednak zbyt mało dokładne do określenia krzywej współ- 


istnienia w bezpośrednim otoczeniu punktu krytycznego. 


Skonstruowanie dylatometru rozcieńczalnikowego, nie posiadającego 
kranów, przez Kumarana i McGlashana (1979) znacznie podniosło dokładność 
pomiarów vE, Dylatometr ten umożliwia zmierzenie dużej liczby punktów w 
funkcji stężenia dla stałej temperatury. Pomiar nadmiarowej objętości 
mieszania lub entalpii mieszania układu z luką mieszalności daje w ob- 
szarze dwufazowym liniową zależność od stężenia co umożliwia interpola- 
cję składów równowagowych (French i wsp., 1979) (metoda ta została wyko- 
rzystana w niniejszej pracy i jest opisana niżej), lecz w miarę docho- 
dzenia do temperatury krytycznej uzyskane wyniki obarczone są coraz 
większym błędem. Uwagi powyższe odnoszą się również do wykorzystania w 


tym celu pomiarów ciepeł mieszania. 


Dylatometr rozcieñczalnikowy może być również przyrządem do precy- 
zyjnego określenia samego wierzchołka krzywej współistnienia (Łuszczyk, 
1988). Sposób ten polega na określeniu x w pracujacym dylatometrze gdy 
pomiar prowadzony jest w izotermie, powiedzmy о dwa stopnie, poniżej 
temperatury krytycznej. Gdy w naczyniu następuje rozwarstwienie się roz- 
tworu, podnoszona jest temperatura aż do zaniku drugiej fazy, po czym 
następuje chłodzenie do ponownego pojawienia się drugiej fazy. Operacje 
te powtarza się do uzyskania dokładnej temperatury przejścia, po czym 
powraca się do temperatury pierwotnej, uzyskuje nowy skład i powtarza 
powyższą procedurę. Tak postępując przesuwamy się zwiększając stężenie 
jednego ze składników aż do obszaru homofazowego. W ten sposób, w jednym 
pomiarze uzyskujemy pomiar nadmiaru objętości i wierzchołka krzywej 
współistnienia wraz z dokładnym określeniem temperatury i składu kry- 
tycznego, gdyż skład roztworu w tym przyrządzie określane są bardzo 
dokładnie (Chyliński, 1987). Podstawowym warunkiem powodzenia tego spo- 
sobu jest pewność powrotu do temperatury bazowej, w której prowadzony 


jest cały pomiar. 


Wpływ zanieczyszczenia próbek na kształt krzywej współistnienia i 
temperaturę krytyczną był badany przez wielu autorów. Np. Tveekrem i Ja- 
cobs (1983) przeprowadzili takie badania dla układu (metanol(1) + cyklo- 
heksan(2)) dodając kontrolowane, niewielkie ilości wody. Stwierdzili 
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oni, ze zanieczyszczenie 0.1% obj. мода powoduje wzrost Те о 3.8 K. 
Zmienial sie również skład krytyczny, а wykładnik krytyczny В wzrastał. 
Innym, ciekawym podejściem do tego problemu jest wykorzystanie zjawiska 
wymiany izotopowej na przykładzie układu (kwas izomasłowy deuterowany (1) 
+ woda(2)) (Gulari i wsp., 1980; Knobler i Scott, 1982). Układ dwuskład- 
nikowy przechodzi w trój- a nawet czteroskładnikowy. Wpływ podstawienia 
izotopowego na temperaturę krytyczną jest znaczący. Pojawienie się deu- 
terowanej wody można traktować jako zanieczyszczenie układu. Knobler i 
Scott (1982) dają szczegółową analizę zachowania się krzywej współist- 
nienia w takim układzie. 


Diagramy fazowe układów binarnych pod wysokimi ciśnieniami stanowią 
częsty obiekt pomiarów. Wydaje się, że to, czy luka powiększa się, czy 
maleje ze wzrostem ciśnienia wynika z oddziaływań wiązaniem wodorowym 
pomiędzy cząsteczkami tego samego i różnego rodzaju. Na przykład układy 
(2-metoksyetanol(1) + alkany(2)) (Bijl i wsp., 1983) i {шефапо1(1) + cy- 
kloheksan(2)) (Singh i van Hook, 1987b) wykazują powiększanie się luki w 
miarę wzrostu ciśnienia natomiast (2-butanol(1) + woda(2)) (Moriyoshi i 
wsp., 1975) - odwrotnie. (Patrz przykład przy omawianiu modelu W-V, 
rozdz. IV-4.2). 


Nukleacja jest zjawiskiem polegającym na samorzutnym rozpadzie ho- 
mogenicznego roztworu na dwie fazy ciekłe. Pomiar temperatury tego roz- 
ради daje w rezultacie krzywą spinodalną. W połączeniu z krzywą współ- 
istnienia uzyskuje się obszar metastabilny. Istotnym problemem jest spo- 
sób wymuszenia "przegrzania" roztworu bez naruszenia mechanicznej sta- 
bilności roztworu. Stosowane są w tym celu: wyładowania elektryczne, 
wzrost ciśnienia, i.t.p. (Howland i wsp. 1980, Clerk i Sengers, 1983). 


Według teorii stosunek ST/AT powinien pozostawać stały wzdłuż krzy- 
wej współistnienia. ôT jest różnicą temperatur pomiędzy krzywą binodalną 
a krzywą spinodalną natomiast AT - różnicą pomiędzy temperaturą krytycz- 
ną i temperaturą przejścia na krzywej binodalnej. (Oczywiście dla danego 
składu. ) Badania realnych układów w bezpośrednim otoczeniu punktu kry- 
tycznego nie potwierdzają tej tezy. Howland i wsp., (1980) zrobili zes- 
tawienie dla kilku układów i stwierdzili, że w funkcji temperatury гге- 
dukowanej e zależność jest nieliniowa. Tłumaczone to jest utrudnioną dy- 
fuzją i powolnym powstawaniem kropel drugiej fazy w pobliżu punktu kry- 
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tycznego. 


V.2.1. Analiza i opis krzywej wspótistnienia 


Uzyskanie poprawnego opisu krzywej wspótistnienia ma podstawowe 
znaczenie w analizie gdyż w gruncie rzeczy celem jest wyznaczenie skła- 
dów faz w równowadze. Dopiero ta informacja jest przydatna w dalszych 
rozważaniach nad punktem krytycznym. W czasach, gdy komputer był narze- 
dziem mało dostępnym najczęściej składy równowagowe odczytywano z gra- 
ficznej prezentacji krzywej współistnienia. Na tak wyznaczonych składach 
Cox i Herington (1956), Stein i Allen (1973) prowadzili analizy obu 
gałęzi krzywej współistnienia osobno w otoczeniu punktu krytycznego. W 
międzyczasie pojawiły się jednak metody numeryczne pozwalające opisać 
krzywą binodalną uzyskaną metodą rozpuszczalnościową (Malesińska, 1960a; 
Malesińska i Malesiński, 1960b; Spinoli i Ferloni, 1974). Metody te 
opierały się na symetryzacji krzywej współistnienia i odłożenia jej w 
funkcji różnicy temperatury od temperatury krytycznej. Rogalski i Stry- 
jek (1980) zmodyfikowali równanie Ferloniego i Spinoliego tak, aby można 
było jednocześnie symetryzować oba ramiona krzywej współistnienia. 


Dla układów ciekłych prawa wykładnicze zostały zaadoptowane przez 
Ley-Koo i Green (1977), natomiast Wegner (1972) wprowadził dodatkowe 
człony pozwalające na rozszerzenie zakresu temperatur dla stosowania 
tych równań. Wykazał on, że wykładniki dalszych członów rozwinięcia są 
sumą wykładnika pierwotnego В i wielokrotności wykładnika korekcy jnego 
A, równego 1/2. 


A >» AS (B+iA) E 
|x | «pos V-1 


Prawa wykładnicze opisują również symetralna krzywej współistnienia z 
wykładnikiem krytycznym а tym samym, który jest odpowiedzialny za zacho- 
wanie się ciepła właściwego w otoczeniu punktu krytycznego mieszalności. 


(x' +ох', )/2 =a + ае + T E MUN V-2 


ї=2 


Powyzsze dwa rownania pozwalaja, poprzez proste zsumowanie, uzyskaé Баг- 


dzo dobre równanie do opisu krzywej wspótistnienia, pozwalajace korzys- 
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taé 2 danych rozpuszczalnošciowych, które sa najbardziej rozpowszechnio- 


ne w literaturze 


B <(B+iA) 


A е101-92+00-2001, Ета 


1-2 1-0 


х=А+Асє + 
0 1 


gdzie «,B - wykładniki krytyczne (mogą być parametrami dobieralnymi), a, 
- parametr będący jednocześnie estymatą składu krytycznego, f - mnożnik 
równy 1/2 dla x > Xo -1/2 dla x < xc. Ilość współczynników A i B zależy 
od zakresu temperatur w jakim znajdują się dane. Równanie powyższe sta- 
nowi dobre narzędzie interpolacyjne przy opracowywaniu danych rozpusz- 
czalnościowych jak również równowagowych. Porownano je z powszechnie 
stosowanymi równaniami Coxa (Cox i Herington, 1956) oraz Francosiniego 
(Spinoli i Ferloni, 1974) na około 50 układach różnego typu (Łuszczyk i 
Stryjek, 1984b) i uzyskano conajmniej dwukrotne zmniejszenie odchylenia 
standardowego oraz brak jakichkolwiek odchyleń systematycznych, które to 


są obserwowane przy stosowaniu wymienionych równań. 


Podobne równania przedstawili niezależnie Viner i Smith (1981) i 
Vnuk (1983a,b), z tym, Ze Vnuk podał również równanie do opisu luki zam- 
kniętej. Ich zdolność do opisu krzywej współistnienia jest podobna do 


zdolności omówionego równania. 


V.3. CIEPŁO WŁAŚCIWE 


Pomiary ciepła właściwego w otoczeniu punktu krytycznego są, rela- 
tywnie do innych wielkości termodynamicznych, dość często spotykane w 
literaturze. Wynika to z tego, że ciepło właściwe pod stałym ciśnieniem 
dla składu krytycznego 208. jest odpowiednikiem ciepła właściwego w 
stałej objętości (С) uładów jednoskładnikowych (Scott, 1978). Z teorii 
wynika, że logarytmiczna dywergenc ja Ср, х jest rządzona wykładnikiem «с, 
którego przewidywana teoretycznie wartość wynosi ok. 0.1. Wśród badanych 
układów są: {kwas izomasłowy(1) + woda(2)) (Pelger і wsp., 1981а,5; 
Klein i Woermann, 1975), (perfluorometylocykloheksan(1) + czterochlorek 
wegla(2)) (Pelger i wsp., 1981a; Schmidt i wsp., 1959), (trójetyloamina 
(1) + woda(2)) (Bloemen i wsp., 1980), 4dwusiarczek wegla(1) + acetoni- 
tryl(2)) (Govindarajan i wsp., 1973). Z powyzszego, niekompletnego zes- 
tawienia wynika, że badane pod tym względem układy reprezentują pełną 
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(1)Trojetylosmina - woda 

(2)гтајегуоаицаа - woda - ctanol(XE :0.07) 
(3) Lrójctyloamina - woda - ctanol(XE =0.10) 
(4)Trójetyloamina - 020 


= 
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Rys.V-1. Zależność ciepła właściwego od temperatury dla układu (trójety- 
loamina(1) + woda(2)) 


gamę oddziaływań międzycząsteczkowych. 


Na rysunku V-1 przedstawione jest typowe zachowanie się ciepła 
właściwego w funkcji temperatury w układach ciekłych na przykładzie po- 
miarów układu (trójetyloamina(1) + woda(2)). Ponadto rysunek ilustruje 
wpływ podstawienia izotopowego (0,0) na temperaturę krytyczną, a także 
destrukcyjny wpływ zanieczyszczeniem trzecim składnikiem (etanol) na za- 
leżność Ch, х od temperatury. 

W pracy Pelgera i wsp. (1981а) wykonano pomiary dla dwóch układów 
posiadających górny (c-C6F11CF3 + CC1,) i dolny punkt krytyczny (kw. 
i-masłowy + н,0) i stwierdzono zgodność, w granicy błędu eksperymental- 
nego, zachowania się układu z przewidywaniami teoretycznymi zarówno co 
do znaków amplitud (różnych dla górnego i dolnego punktu krytycznego) 
jak i wykładnika a. 


Bloemen i wsp. (1980) dla układu trójetyloaminy z ciężką wodą 
(LCST) stwierdzili zgodność z postulowaną stałością stosunku amplitud 
(Levelt-Sengers, 1983) obliczonych powyżej i poniżej To Pomiary dla te- 
go układu wykonano kalorymetrem skaningowym co dawało ok. 20000 punktów 
pomiarowych na serię. Osiągnięto również w ten sposób granicę możliwości 
podejścia do punktu krytycznego. (Wartość najbliższej temperatury pomia- 


ru wyrażonej w parametrze zredukowanym € wynosiła 7.10 *, со w skali 
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normalnej daje utamek milikelwina). 


V.4. ENTALPIA MIESZANIA 


Skonstruowanie skaningowego kalorymetru rozcieñczalnikowego, obok 
przepływowego, pozwoliło na wykonywanie dużej liczby dokładnych pomiarów 
izotermicznych. Ewing i wsp. (1985) wykonali, bodajże pierwsi, serię po- 
miarów entalpii mieszania (НЕ) w bezposrednim otoczniu punktu krytyczne- 
go układu (metanol(1) + cykloheksan(2)). Łącznie zmierzono 11 izoterm w 
zakresie 1.3 K, z których najbliższa temperaturze krytycznej była 6 mK 
powyżej Te: Istnieje pewne podobieństwo pomiędzy zachowaniem НЕ 1 nad- 
miarowej objeto$ci mieszania. Орје te funkcje sa pochodnymi energii swo- 
bodnej Gibbsa. Są uzyskiwane z pomiarów bezpośrednich i jeśli sa one wy- 
konane z dużą dokładnością, możliwe jest na ich podstawie wyliczenia 
drugiej pochodnej względem składu. Wielkość ta jest stosowana z kolei do 
weryfikacji zachowania się układu w otoczeniu punktu krytycznego. Rezul- 
taty swych pomiarów autorzy przeanalizowali z zastosowaniem parametrycz- 
nego równania stanu i wyznaczenia, drogą pośrednią, wykładników krytycz- 
nych В i у. (Ewing i wsp., 1986). 


Podobne pomiary, w tego samego typu kalorymetrze, przeprowadzili 
Marsh i Rogers (1989) dla układu ínitroetan(1) + cykloheksan(2)). Auto- 
rzy ci zmierzyli 6 izoterm w zakresie 2 K, z czego jedna była poniżej 
punktu krytycznego natomiast najbliższa T. - 12 mK powyzej. W pracy tej 
szczególną uwagę zwrócono na oczyszczenie składników. Np. zawartość моду 
w nitroetanie badano spektroskopią IR i oceniono jej stężenie na 107° 
ułamka molowego. Również tutaj przeprowadzono analizę wyników stosując 
parametryczne równanie stanu, lecz autorzy nie podali wyników numerycz- 


nych. 
V.5. NADMIAROWA OBJĘTOŚĆ MIESZANIA 

Nadmiarowa objętość mieszania (УЕ) jest, podobnie јак епфа1р1а 
mieszania,  wielko$cia  pochodna  entalpii swobodnej. Badania V w 
otoczeniu punktu krytycznego prowadzone były dla niewielu układów. Tak 


więc badane były: (metylocykloheksan(1) + perfluorometylocykloheksan(2)) 
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(Kumaran i McGlashan, 1980), (trójetyloamina(1) + мода(2)), (Letcher 1 
Spitieri, 1983), (heksan(1) + perfluoroheksan(2)) (Block i wsp., 1981). 
Wszystkie te uktady badane byty przy uzyciu dylatometrów rozcienczal- 
nikowych. Wcześniejsze badania tego ostaniego układu (Gaw i Scott, 1971) 
wykonane były metodą densymetryczną i wymagały powtórzenia, chociaż me- 
tody densymetryczne również są stosowane do badań w otoczeniu punktu 
krytycznego mieszalności (Widom i Khosla, 1980). 


Dla izotermy najbliższej temperatury krytycznej możliwa jest próba 
obliczenia wykładnika A. 


(УЕ) = Е +E (х-х ) + E |х-х ЈА V-4 
X. 0 1 c 2 c 


(ауё/ах) = E + Е (x-x )/Ix-x |·|х-х 10-1) ү-5 
х 3 4 C C C 


c 

Powyzsze równania zostaty wyprowadzone ze skalowanego równania sta- 
nu. Wartość wykładnika A=2.8 oraz jego zależność od temperatury i steze- 
nia jest szczegółowo dyskutowana przez Blocka i wsp. (1981). Dla układu 
badanego w tej pracy wartość tego wykładnika wynosi =2.84+0. 03 
(Łuszczyk, 1988), czyli jest zgodny, w zakresie błędu eksperymentalnego, 


z wartością teoretyczną. 


Niejako naturalnie związany z pomiarem VE jest wspótczynnik rozsze- 
rzalno$ci termicznej. Bada sie go za pomocą dylatometrów. W myśl teorii 
punktu krytycznego powinien on wykazywać dywergencję, której amplituda 
powiązana jest z linią krytyczną układu dwuskładnikowego. Według Kleina 
i Woermanna (1976) współczynnik ten przewiduje poprawnie zależność 
ciśnieniową temperatury krytycznej (dT_/dP) dla wielu układów. 


V.6. ZJAWISKA TRANSPORTOWE 


Do zjawisk transportowych powszechnie zalicza się lepkość, dyfuzję 
i przewodnictwo cieplne. Najczęściej badaną własnością spośród wymienio- 
nych jest lepkość, która jest mierzona wprost. Własność ta wykazuje bar- 
dzo duży wzrost w bezpośrednim otoczeniu punktu krytycznego. Najczęściej 
badana jest w funkcji temperatury dla składu bliskiego krytycznemu. Na 
rysunku V-2 przedstawiona jest przykładowa "lepkość  nadmiarowa" (јако 
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Rys.V-2. Ilustracja zaleznosci "nadmiarowej lepkošci" od temperatury i 
składu dla układu (3-MOPN(1) + woda(2)) 


różnica pomiędzy lepkością zmierzoną a udziałami wniesionymi przez czys- 
te składniki) układu {3-МОРМ(1) + woda(2)) (Łuszczyk 1 Stryjek, 1978). 
Widać tu wyraźne "przygotowywanie" się układu do przejścia przez punkt 
krytyczny dość daleko od temperatury krytycznej. Wzrost lepkości 
związany jest ze wzrostem długości korelacji i powstawaniem klusterów, 


dla których rośnie energia swobodna przepływu Є". 


М literaturze istnieje niewiele danych о takich pomiarach, а ist- 
niejace dotycza niewielkiej liczby uktadów. Мупіка to 2 trudnosci ekspe- 
rymentu. Dokładne pomiary można uzyskać właściwie м wiskozymetrach typu 
Oswalda, odpowiednio zmodyfikowanych, np. takich, które opisano w Book 
of ASTM Standards (1969). Ponieważ jest to układ otwarty, trzeba liczyć 
na utratę bardziej lotnego składnika. Stąd do pomiarów wybierane są 
układy posiadające temperaturę krytyczną bliską pokojowej, jak np. układ 
(anilina(1) + cykloheksan(2)) (D'Arrigo i wsp., 1977), którego tempera- 
tura krytyczna wynosi ok. 303 K. Aby uniknąć zmiany składu Sada i wsp. 
(1977) zbudowali wiskozymetr гемегѕујпу dla pomiarów układu {піїгођеп- 
zen(1) + pentan(2)) (Тт ok. 300 K; temperatura wrzenia репїапи ok. 310 
K). Po przygotowaniu sktadu krytycznego wiskozymetr byt zatapiany. Po- 
miar prowadzono obracając wiskozymetr i mierząc czas przepływu w obie 
strony kapilary. Innym sposobem na zapobieżenie utraty jednego ze 


składników jest stosowanie dużej ilościowo próbki połączonej z rezerwu- 
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агет wiskozymetru. 


Analiza danych lepko$ci kinematycznej pozwala wyznaczyé lepkošciowy 
wykładnik krytyczny v. Podstawową trudnością jest wydzielenie lepkości 
podstawowej, gdyż lepkość w otoczeniu punktu krytycznego jest traktowa- 
na, podobnie jak inne statyczne wielkości termodynamiczne, jako 
składająca się z dwóch udziałów: regularnego 1 anomalnego. Również z ро- 
miarów lepkości można uzyskać mieszany współczynnik dyfuzji. Najczęściej 
pomiary lepkości łączone są z pomiarami rozpraszania światła, gdyż na 
ich podstawie można określić długość korelacji i odpowiednie wykłdniki 
krytyczne y i у. Różnica między tymi wykładnikami daje z kolei wyktadnik 
rządzący dywergencją współczynnika dyfuzji. 


Innym, dość nietypowym zastosowaniem pomiarów lepkości w otoczeniu 
punktu krytycznego jest wykorzystanie ich do określenia krzywej 
współistnienia i obliczenia wykładnika krytycznego B (Dega-Dałkowska, 
1980). 


V.7. INNE WIELKOŚCI FIZYKOCHEMICZNE 


Poza wymienionymi wielkościami termodynamicznymi również inne wiel- 
kości fizykochemiczne wykazują osobliwości w otoczeniu punktu krytyczne- 
go. Są to własności pochodne, lub w jakikolwiek sposób zależne od 
ściśliwości lub ciepła właściwego. Ciekawe zależności obserwuje się dla 
tzw. nieliniowego efektu dielektrycznego (NDE). W otoczeniu punktu kry- 
tycznego następuje wzrost podatności dielektrycznej. Zjawisko to było 
badane szczegółowo przez Piekarę (1939) i badania tego efektu, co bardzo 
ciekawe, prowadzone są głównie w Polsce (Majgier-Baranowska, 1980; Maj- 
gier-Baranowska i wsp., 1981). Pnadto obserwuje się w otoczeniu punktu 
krytycznego wzrost stałej Kerra (Pyżuk i wsp., 1979, 1980). Badania 
własności elektrycznych pozwalają przewidzieć potencjalne wystąpienie 
punktu krytycznego, gdyż np. efekt NDE pojawia się dość daleko od tempe- 
ratury krytycznej” А 


* Np. na podstawie zależności МОЕ od składu w temperaturze pokojowej dla 


układu {nitrometan + acetonitryl) sam stwierdziłem rozpad na dwie fazy w 


temperaturze ok 275 K. 
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Zjawiskiem nierozerwalnie zwiazanym 2 punktem krytycznym jest roz- 
praszanie światła. Pomiar kątowego rozpraszania wymaga dość skomplikowa- 
nej 1 specjalizowanej aparatury, a jednak należy do częściej муКопума- 
nych w badaniach punktu krytycznego. Sam pomiar jest również skompliko- 
wany i należy uwzględnić wiele czynników zakłócających takich jak roz- 
praszanie wielokrotne, geometria próbki i.t.d. Przyczyną jest prawdopo- 
dobnie spora ilość informacji, jaka może być uzyskana z tych pomiarów 
(Chang i wsp. 1979). Do najważniejszych należy wzrost długości korelacji 
i wykładnik nią rządzący. Inną wielkością wykazującą nietypowe zachowa- 
nie jest pochłanianie ultradźwięków, który to efekt jest zależny od 
częstotliwości. Tłumaczone to jest anomalnym zachowaniem sie współczyn- 
nika ściśliwości i związanymi z tym fluktuacjami gęstości (Jaschul i 
wsp. 1984). 
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VI CZEŠC DOSWIADCZALNA 


VI.1. WYBÓR UKLADU. POMIARY WSTEPNE 


Јако przedmiot badań wybrano układ 3-metoksypropionitryl (3-MOPN) + 
мода. Temperatura krytyczna mieszania tego układu jest nieco wyższa od 
pokojowej; oba składniki mają zbliżone gęstości w otoczeniu punktu kry- 
tycznego, co redukuje znacznie efekt grawitacyjny; podobne współczynniki 
załamania światła powodują, że opalescencja krytyczna nie zaburza obser- 
wacji współistniejących faz, chociaż stwierdzenie występowania powierz- 


chni międzyfazowej jest utrudnione. 


Innym czynnikiem predestynującym układ do badań jest łatwa dostep- 
ność składników. Woda nie wymaga komentarza, natomiast 3-MOPN może być 
prosto zsyntetyzowany z łatwo dostępnych składników. Oczyszczanie jego 
jest również proste. Metoda syntezy i oczyszczania jest opisana w dal- 
szej części rozdziału. 3-MOPN jest związkiem mało znanym, tak więc jego 
budowa, własności i zachowanie się w innych układach są opisane we 


wstępie rozdziału VIII poświęconego dyskusji wyników pracy. 


VI.1.1. Krzywa współistnienia 


Krzywa rozpuszczalności tego układu była badana uprzednio (Stryjek 
i wsp., 1978b; Łuszczyk i Stolarska, 1984) metodą ampułkową. Wyniki 
przedstawiono na rys. VI-1. Ponadto, przy okazji pomiarów nadmiarowej 
objętości VE (Łuszczyk, 1988) wyznaczono składy równowagowe dla szeregu 
temperatur (na rys. VI-1. uzyskane składy równowagowe są połączone li- 
nią) oraz dokładnie określono temperaturę i skład krytyczny metodą opi- 
saną w rozdziale V-2. Rezultaty, dla porównania tych dwóch metod, podane 
są również na rys. VI-1. Dane Stryjka i wsp. (1978b) pokrywają większy 
zakres temperatur, lecz składy równowagowe uzyskane zostały przez inter- 
polację graficzną danych rozpuszczalnościowych. Co prawda istnieje pewna 


rozbieżność wyników w zakresie niższych temperatur, ale jest to zrozu- 
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miałe przy stosowaniu różnych metod. Istnieje jednak zupełnie dobra zgod- 


ność w zakresie temperatur istotnych z punktu widzenia tej pracy. 


***** |uszczyk,1988 

310.00 Luszczyk,1988 top 
+++++ Luszczyk,Stolarska, 1984 
+++++ Stryjek i wsp.,1978b 


300.00 


Temperatura (K) 


290.00 


280.00 
0.00 0.20 0.40 0.60 0.80 1.00 


x1 


Rys.VI-1. Krzywa współistnienia układu (3-MOPN(1) + woda(2)). 


Tablica VI-1. 


Składy równowagowe dwóch faz ciekłych układu (3-MOPN(1) + woda(2)) wyz- 
naczone na podstawie pomiarów nadmiarowej objętości mieszania. 


faza dolna faza górna 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


ооооооооо 
ооооооооо 
ооооооооо 
ооооооооо 


Krzywą współistnienia uzyskaną w trakcie pomiarów VE opisano 
równaniem V-3 omówionym w rozdz. V.2.1. W tablicy VI-1. zamieszczone sa 


numeryczne wartości tak uzyskanych składów równowagowych, jak i obliczo- 
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ne z wygładzenia danych. Te ostatnie wykorzystane będą w dalszych kore- 
lacjach. Parametr e jest temperaturą zredukowaną (e = |T = TL Zt któ- 
ry jest jednocześnie miarą "odległości" od punktu krytycznego. Przypom- 
nieć należy, że najlepsze pomiary metodą ampułkową pozwalaj otrzymać 
wartość tego parametru rzędu 107°. Ропїема2 krzywizna funkcji VE zdaza 
do 0, wyznaczone sktady w bezposrednim sasiedztwie temperatury krytycz- 
nej są obarczone błędem, dyskutowanym w pracy Łuszczyka (1988). Jednakże 
błąd ten w skali rysunku jest niewidoczny. 


0.00 


—0.20 


VE (cm3/mol) 


—0.30 


298.133 


311.146 
- 0.40 
0.00 0.20 0.40 0.60 0.80 1.00 


Rys.VI-2. Nadmiarowa objętość mieszania układu (3-MOPN(1) + woda(2)) w 
temperaturze 298.133 i 311.146 K. 


VI.1.2. Nadmiarowa objętość mieszania 


Nadmiarowa objętość mieszania (УЕ) układu była oznaczona uprzednio 
metodą dylatometryczną (Łuszczyk i Stryjek, 1978) dla temperatur powyżej 
temperatury krytycznej. Uzyskane wyniki są obarczone znacznym błędem 
spowodowanym konstrukcją zastosowanego dylatometru. Pomiary powtórzono w 
trakcie pracy stosując dylatometr rozcieńczalnikowy Kumarana i McGla- 
shana (1979) oraz metodykę pomiaru opracowaną przez Chylińskiego (1987). 


Wyniki numeryczne pomiarów zostały opublikowane (Łuszczyk 1988) i nie są 
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zamieszczone w tej pracy. Na rysunku VI-2. zaprezentowana jest zależność 
УЕ od składu w temperaturach 298.133 i 311.146 K, to jest znacznie 
poniżej i 3 mK powyżej temperatury krytycznej mieszalności. 


VI.2. PRZEGLĄD METOD POMIARÓW RÓWNOWAGI CIECZ-PARA 


W literaturze opisano szczegółowo wiele metod i aparatów stosowa- 
nych do pomiarów prężności pary (Malanowski, 1974; Aim i Hala, 1988). 
Generalnie podział metod pomiaru równowag ciecz-para wynika ze sposobu 


pracy aparatu i dzieli się na metody dynamiczne i statyczne. 


VI.2.1. Metody dynamiczne 


W metodzie dynamicznej następuje obieg fazy parowej i ciekłej w a- 
paracie na skutek wrzenia substancji. Ciśnienie układu jest utrzymywane 
na stałym poziomie i mierzona jest temperatura wrzenia. Metoda wymaga 
również oznaczenia przynajmniej składu fazy ciekłej pozostającej w гби- 
nowadze z parą, jeśli pomiar dotyczy mieszanin. Szczegółowy przegląd 
tych metod przedstawił Malanowski (1982). Metody tej grupy nie nadają 
się do pomiaru układów tworzących dwie fazy ciekłe. 


VI.2.2. Metody statyczne 


Zasada metody statycznej polega na pomiarze ciśnienia całkowitego 
nad roztworem w stałej temperaturze. Chociaż zasada jest bardzo prosta, 
to zasadnicze trudności wynikają z konieczności operowania substancjami 
całkowicie odgazowanymi. W literaturze opisano wiele sposobów odgazowa- 
nia, przechowywania i przenoszenia substancji do aparatury. Najbardziej 
typowym postępowaniem jest odgazowanie substancji poprzez wielokrotne 
zamrożenie, odpompowanie do próżni i rozmrożenie lub wręcz destylację 
substancji pod wysoką próżnią. Następnie próbki przemrażane są do 
ampułek, a przed pomiarem są przenoszone do aparatu również poprzez 
przemrożenie. Ten sposób postępowania jest względnie prosty w przypadku 
operowania substancjami łatwo wrzącymi. Trudności rosną w miarę wzrostu 


temperatury wrzenia i w celu ich ominięcia opracowano kilka metod prze- 
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noszenia Substancj do naczynia równowagowego np. za ротоса Systemu biu- 
ret (Murray i Martin, 1975) сту tez zastosowania specjalnych strzykawek 
dozujacych. Ta ostatnia metoda, opisana w wielu publikacjach przeglado- 
wych (Aim i Hala, 1988), wydaje sie szczególnie atrakcyjna ze wzgledu na 
możliwość uzyskania dużej liczby punktów eksperymentalnych w funkcji 
składu w jednym cyklu pomiarowym. Z perspektywy czasu wydaje się, że 
trudności techniczne tej metody są duże gdyż nie obserwuje się w litera- 
turze jej rozwoju. Poniekąd potwierdzają to obserwacje poczynione w 


trakcie tej pracy i opisane niżej. 


VI.3. WYBÓR I OPIS ZASTOSOWANEJ METODY 


Jak wiadomo, nie można zbudować uniwersalnej aparatury do pomiarów 
równowag ciecz-para. O wyborze metody decyduje cel, w jakim wykonuje się 
pomiary. W niniejszej pracy równowagę ciecz-para zmierzono metodą sta- 
tyczną. Pomiary tą metodą uznawane są za bardzo dokładne, pod warunkiem, 
że operuje się substancjami całkowicie odgazowanymi, a temperatura 
układu jest stabilna i dokładnie zdefiniowana. Dokładność wynika z ist- 
nienia prawdziwej równowagi w układzie pozbawionym przepływów mas i cie- 
pła. Ponadto metoda nie wymaga analizy współistniejących faz. Wielkości 
te są obliczane w trakcie opracowywania danych na podstawie informacji o 
masach wprowadzonych substancji. W przypadku badanego układu analiza 
składu którejkolwiek z faz byłaby bardzo trudna (podobne gęstości i 
współczynniki załamania światła, bardzo trudna analiza chromatograficz- 
na). Podstawową wadą metody statycznej jest jej czasochłonność, lecz 


znów, ze względu na cel pracy nie był to czynnik istotny. 


Celowe wydawało się zbudowanie urządzenia, z dozowaniem ciągłym 
składników czystych ze specjalnych dozowników, pracującego w warunkach 
izotermicznych. Niestety, próby wykazały, że nie ma gwarancji absolutnej 
szczelności aparatury. Wymagane przejścia szkło — metal z przetopami ko- 
varowymi oraz uszczelnienia elastomerami, konieczne ze względu na deli- 
katność aparatury, nie zdały egzaminu. Sprawdzono wiele rodzajów tworzyw 
"O"-ringów lecz za każdym razem długotrwałe pomiary wykazywały wnikanie 
par wody (która jest medium termostatującym) do naczynia równowagowego. 
Być może, inny system pomiaru ciśnienia pozwoliłby na wykonanie aparatu- 


ry z metalu, co umożliwiłoby zastosowanie dozowania ciągłego. 
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Wobec powyzszego zdecydowano sie zbudowaé aparat, wzorujac na opi- 
sanym w pracy Heady i Cahn'a (1973). Jest to wersja "klasyczna", opisana 
w literaturze przez wielu autorów. Wprowadzono једпак kilka istotnych 


ulepszeh, przedstawionych w dalszej cze$ci. 


VI. 4. OPIS STANOWISKA BADAWCZEGO 


Zestaw do pomiarów prężności par metodą statyczną zbudowany jest z 
trzech zasadniczych części: termostatowanego manometru wewnętrznego z 
dołączonym naczyniem równowagowym, manometru zewnętrznego oraz linii wy- 
sokiej próżni. Schemat układu przedstawiony jest na rysunku VI-3. Posz- 


czególne elementy układu omówione są w/g oznaczeń na rysunku. 


Tm - TERMOSTAT został szczegółowo opisany w pracy doktorskiej K. Chyliń- 
skiego (1987). Tutaj opisane są tylko modyfikacje i usprawnienia. Zbudo- 
wany jest on w formie prostopadłościanu ze szkła epoxydowego o wymiarach 
350x600x750 mm i pojemności ok. 150 dm. Grubość ścianek wynosi 22 mm. W 
ścianach przeciwległych o najdłuższym boku wstawione są okna z szyb 
płasko-równoległych, o grubości 6 mm, o wymiarach 150x400 mm. Ściany 
termostatu pokryto dodatkowo sklejką o grubości 20 mm w celu podwyzsze- 
nia stabilności temperatury. Termostat umieszczony jest na podnośniku 
umożliwiającym opuszczanie zbiornika przy nieruchomej pokrywie, do któ- 
rej umocowana jest również cała aparatura. Jako układ wymuszający obieg 
cieczy termostatującej zastosowano cztery turbiny, zamocowane do рокгу- 
wy, umieszczone w narożach. Prędkości obrotowe turbin są parami różne 
dla uzyskania przepływu wody termostatującej zarazem turbulentny i lami- 
narny co zapewnia maksymalną stabilność temperatury. Naroża dna 
zaokrąglono kształtkami z blachy aby zapewnić łatwiejszą cyrkulację wo- 
dy. Termostat wyposażony jest w grzałkę o mocy 1400 W, oraz w chłodnicę 
wodną. W przeciwległym końcu (od strony wolniejszych turbin) umieszczona 
jest czujka regulatora temperatury (UNIPAN typ 650), przy pomocy którego 


jest kontrolowana temperatura. 


Wnętrze termostatu oświetlone jest przez tylną szybę lampą jarze- 
niową o mocy 2x8 W z matowym kloszem. Położenie lampy jest regulowane w 


sposób płynny. 
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Rys.VI-3. Schemat aparatury. Objaśnienia oznaczeń w tekście. 


Mw - MANOMETR WEWNĘTRZNY zbudowany jest z kalibrowanych rur ze szkła py- 


reksowego o precyzyjnej średnicy wewnętrznej 22 mm. Ramię od strony na- 
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czynia pomiarowego Np jest krótsze w celu zmniejszenia objeto$ci fazy 
parowej. W miejscu zwężenia umieszczony jest znacznik wewnętrzny Mr wyz- 
naczający część kalibrowaną i niekalibrowaną. Objętość niekalibrowana 
jest wyznaczana każdorazowo po zmianie naczynia równowagowego, natomiast 
kalibrowana - z geometrii manometru. Dłuższe ramię manometru dołączone 
jest do kanału próżniowego PR1 poprzez kran К1. Dolna część manometru 
podłączona jest do rezerwuaru rtęci R o pojemności ok. 300 cm. Rezerwu- 
ar jest odcinany od manometru poprzez kran K2, natomiast odpowietrzanie 
i zapowietrzanie rezerwuaru zapewniają krany КЗ i K4 łączące go z 
kanatem prózniowym. Kran K2 zabezpiecza przed cofaniem sie rteci до re- 
zerwuaru pod wpływem ciśnienia panującego w układzie. Ma to szczególnie 
znaczenie przy pomiarach ciśnień zbliżonych do atmosferycznego. Oba ra- 
miona manometru wyposażone są w maski (nie przedstawione na rysunku) w 
postaci metalowych, cienkościennych i poczernionych walców miedzianych o 
długości ok. 100 mm ze szczelinami w części centralnej o szerokości 1 
mm, przesuwane wzdłuż osi manometru i sterowane z zewnątrz. Umożliwiają 
one uzyskanie ostrego obrazu granicy powierzchni rtęci i zapobiegają ге- 


fleksom świetlnym. 


PR1 - LINIA PRÓŻNIOWA części ruchomej aparatury wyposażona jest w wejś- 
cie do podłączenia pomiaru próżni Pc, złącza manometru zewnętrznego, za- 
silania gazowym azotem (K20) oraz urządzenia do tłuczenia ampułek (K13). 
Podłączenie do głównej linii próżniowej, odcinane kranem K5, poprzez 
złącze Z ze szlifami kulistymi pozwala na wyjęcie aparatu wraz z manome- 
trem (Mw) i rezerwuarem (R) z termostatu i zastąpienie go w razie po- 


trzeby innym urządzeniem pomiarowym. 


PR2 - GŁÓWNA LINIA PRÓŻNIOWA wykonana jest ze szklanej rury pyreksowej o 
średnicy ok. 15 mm, do której dołączone są, poprzez wymrażarkę W: dyfu- 
zyjna pompa próżniowa DPP oraz pompa olejowa Po; wejście do zasilania 
układu suchym azotem N2; manometr pomocniczy Mp o długości ok 600 mm 
(służący do kontroli napełniania zbiornika azotu ZN2); wyjście próżniowe 
jednego z ramion manometru zewnętrznego zamykana kranem K16; dwa wypro- 
wadzenia zamykane kranami K7 i K8, zakończone szlifami. Dołączone są one 
do części kanału ograniczonej z obu stron kranami K9 i K10; część ta 
jest wyposażona w wyjście do pomiaru próżni Pp i służy do odgazowania 


substancji czystych i przygotowania ampułek. 
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Rys.VI-4. Budowa naczynia pomiarowego (Np). 


Np - NACZYNIE POMIAROWE przedstawione jest na rysunku VI-4. Zbudowane 
jest w kształcie walca ze szkła pyreksowego o wymiarach orientacyjnych 
55 x 70 mm. Wewnątrz wbudowane jest mieszadło szklane MM z zatopionym 
rdzeniem magnetycznym w górnej jego części. Mieszanie fazy ciekłej i ра- 
rowej wymuszane jest zewnętrznym magnesem M poruszanym silnikiem elek- 
trycznym z regulowaną prędkością obrotową. Przy odpowiedniej intensyw- 
ności mieszania uzyskiwane jest bardzo szybkie (ok. 15 min.) dojście do 
stanu równowagi fazowej, nawet jeśli były obecne dwie fazy ciekłe. Połą- 
czenie naczynia z manometrem Mw o średnicy wewn. 10 mmm wyposażono w 
małą chłodniczkę wodną Ch, zasilaną z zewnątrz, do zabezpieczenia przed 
destylacją substancji z naczynia na powierzchnię rtęci. Jest to szcze- 
gólnie ważne w trakcie zmiany temperatury w termostacie, w czasie 
którego występują różnice temperatur. W czasie pomiaru chłodnica jest 


odłączona. 


Mz - MANOMETR ZEWNĘTRZNY zbudowany jest, podobnie jak manometr wew- 
nętrzny, z kalibrowanych rur ze szkła pyreksowego o średnicy wewnętrzne j 


22 mm. Oba ramiona o długości roboczej >500 mm zakończone są kranami 
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próżniowymi, z których К11 podłączony jest do głównego kanału próżnio- 
wego, a K12 do układu zasilania azotem. Cały manometr zamontowany jest w 
obudowie tekstolitowej z przezroczystymi ściankami. Obudowa umieszczona 
jest na trójkątnym stoliku z poziomnicą i regulowanymi nóżkami. Ramiona 
manometru wyposażone są również w opisane wcześniej maski ułatwiające 
pomiar poziomu rtęci. Wewnątrz obudowy znajduje się termometr rtęciowy T 
o dokładności 0.1 K pozwalający mierzyć aktualną temperaturę. Manometr 
oświetlony jest od tyłu za pomocą 8 W świetlówek z matowym kloszem. Nie 


powoduje to zbędnego ogrzewania manometru. 


Rys.VI-5. Przyrząd do odgazowywania składników mieszaniny. 


Wspomniany system zasilania azotem jest uruchamiany, jeśli w naczy- 
niu równowagowym panuje ciśnienie większe niż 50 kPa (350 mmHg). Dopływ 
azotu ze zbiornika ZN2 poprzez krany K19 i K6 (przy zamkniętym kranie K5 
i otwartym K1 i K12) powoduje rozkład ciśnienia na oba manometry. 
Zwiększa to zakres pomiarowy układu do ok. 100 kPa. 


Ponadto, w przygotowaniu próbek i w czasie pomiaru stosowane są: 
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- urządzenie do odgazowywania składające się z dwóch zbiorników o 
pojemności ok. 150 cn? każdy, połączone ze sobą rurką, do której 
dołączony jest kran zakończony szlifem (rys. VI-5); 

- urządzenie do tłuczenia ampułek (rys. VI-6.) ze szkła pyreksowego 
(dlugości ok. 400 mm) skonstruowane w ten sposób, aby uniemożliwić cofa- 
nie się ampułki w trakcie tłuczenia. Tłuczkiem jest rdzeń magnetyczny 
pokryty teflonem. Przed tłuczeniem koniec ampułki (jej wygląd również 
przedtawia rys. IV-6. ) jest nacinany aby ułatwić stłuczenie. Zrezygnowano 
z ampułek opisanych w literaturze (Malanowski 1974, Wieczorek 1976) ze 
względu na ich pracochłonność. System ten został sprawdzony wielokrotnie 
w niezależnych pomiarach. Podobny sposób wprowadzania substancji był 
stosowany przez Diaz Penę i McGlashana (1959) w pracy o pomiarach 
ściśliwości. 


Rys.VI-6. Przyrząd do tłuczenia ampułek oraz wygląd ampułki. A - ampułka 
przed wprowadzeniem próbki; B-B miejsce zatapiania próbki. 


Do pomiaru wymaganych wysokošci stosowany byt katetometr produkcji 
ZSRR typu KM-8 о długości roboczej 500 mm i rozdzielczości 0.01 mm z 
możliwością oszacowania odczytu do 0.005 mm. Katetometr ustawiony jest 
na niezależnym stoliku i umożliwiał odczyt dwóch manometrów przez obrót 


w osi pionowej bez konieczności ustawiania poziomu podstawy. 


Pomiar próżni określany był przy pomocy próżniomierza McLeoda, 
który umożliwiał oznaczenie próżni statycznej rzędu 1078 mmHg co odpo- 
wiada 1072 Ра. Konstrukcja systemu próżniowego pozwalała mierzyć próżnię 
w kilku punktach aparatury, w miejscu aktualnie wykorzystywanym (patrz 
rys. VI-3. ). 


VI.4.1. Metodyka pomiaru 


Pomiar prężności pary w prezentowanym aparacie składa się z kilku 


etapów i czynności przygotowawczych. 


1. Należy przygotować pewną ilość próbek czystych, odgazowanych składni- 
ków w dużych i małych ampułkach. Pierwsze stosowane są do określenia 
prężności pary czystych składników i jako ilość wyjściowa do dalszych 
pomiarów dla mieszanin. Małe natomiast - do tworzenia mieszanin w па- 


czyniu równowagowym. 


2. Naczynie równowagowe po skończonej serii pomiarów jest odcinane do 
mycia i dlatego przed rozpoczęciem następnych pomiarów należy окге5- 
lić ponownie jego objętość całkowitą tzn. obejmującą przestrzeń do 
znacznika kalibrowanej części manometru wewnętrznego (Mr). Objętość 
tę określano stosując równanie stanu gazu doskonałego. Do одротрома- 
nego naczynia równowagowego, przy opuszczonej rtęci w manometrze wew- 
nętrznym, wprowadzano suchy azot do ciśnienia około 30 kPa i zamykano 
rtęcią. Po ustaleniu temperatury mierzono ciśnienie panujące w ukła- 
dzie. Następnie podnoszono ciśnienie w układzie dołączając manometr 
zewnętrzny. Podniesienie ciśnienia powodowało zmniejszenie objętości 
zajmowanej przez azot. Ponownie układ termostatowano i określano pa- 
nujące ciśnienie. Z danych tych obliczano objętość całkowitą naczynia 


w/g wzoru: 
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V = (p.V;- р,У,)/(р„- р,) VI-1 


gdzie objetosci % i LA obliczone sa z geometrii manometru. Taki spo- 
sób pozwala określić objętość A z dokładnością 0.2 cm3. 


3. Wprowadzenie próbki wyjściowej do układu przebiega następująco: od- 
pompowujemy naczynie pomiarowe i rtęć w rezerwuarze przy otwartych 
kranach K1, K2, КЗ i K4. Po około 2-3 godzinach pompowania krany te 
zamyka się i  zapowietrza część PR1 kanału próżniowego (przy 
zamkniętym K5 i przez K6) i poprzez kran K4 i K2 podnosi się rtęć do 
poziomu takiego, aby uzyskać niewielki prześwit pomiędzy ramionami 
manometru. Następnie zamykamy kran K2, zrównuje się ciśnienie w re- 
zerwuarze z ciśnieniem zewnętrznym i zamyka się kran K4. Dołączamy do 
wejścia K13 przyrząd do tłuczenia ampułek z umieszczoną wewnątrz am- 
putka. Ponownie odciągamy próżnię poprzez K5 i pompujemy tak długo 
(ok. 2 godz.), aż uzyskamy stabilną próżnię rzędu T«10 = 7.10 Ра. 
Czas ten jest potrzebny do usunięcia powietrza z mieszadła teflonowe- 
go, które służy jako tłuczek. Pod naczynie pomiarowe podstawiamy па- 
czynie Dewara z ciekłym azotem. Zamykamy kran К5 odcinając układ od 
pompy i tłuczemy ampułkę. Otwieramy kran manometru K1 i rozpoczyna 
się proces destylowania substancji do naczynia. Proces ten można 
przyspieszyć ogrzewając zbiornik Z przyrządu do tłuczenia ampułek 
(rys. VI-6) np. ciepłą wodą. Po wdestylowaniu całej substancji zamyka- 
my kran K1 a otwieramy delikatnie kran K2 podnosząc rtęć w manometrze 
do odpowiedniego poziomu. W ten sposób substancja została zamknięta 
rtęcią w naczyniu równowagowym. Teraz odpompowujemy rezerwuar i po 
usunięciu ciekłego azotu pozwalamy ogrzać się naczyniu do temperatury 
pokojowej. Podnosimy zbiornik termostatu i stanowisko gotowe jest do 


pomiarów. 


W przypadku, gdy w naczyniu równowagowym znajduje się już substan- 
cja, postępowanie jest podobne, z tym, że po opuszczeniu temostatu naj- 
pierw wymraża się substancje w naczyniu pomiarowym dla zrównania się 
ciśnień w lewym i prawym ramieniu manometru. Następnie odpompowu jemy re- 
zerwuar poprzez krany K4 i K3, a następnie powoli, poprzez kran K2 o- 
puszczamy rtęć w manometrze do takiego poziomu, aby uzyskać spory prze- 
świt między obydwoma ramionami manometru. Kran K1 jest cały czas 


zamknięty. Dalsze postępowanie jest już analogiczne z opisanym dla subs- 
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tancji czystej. 


Pomiary są prowadzone dla każdego składu w funkcji temperatury. Po 
dojściu temperatury do żądanej, wyłącza się mieszanie w naczyniu ротја- 
rowym i po odczekaniu 30 - 45 min. na ustalenie się równowagi temodyna- 
micznej, odczytuje poziomy menisków rtęci w manometrze, poziom znacznika 
części kalibrowanej oraz poziom cieczy w naczyniu. Wykonuje się jeszcze 
dwa takie pomiary w odstępach półgodzinnych i zmienia temperaturę 
układu. Stosunkowo długi okres pomiędzy pomiarami daje pewność uzyskania 
stanu równowagi. Testy 12-to i 24-ro godzinne nie dawały zmian różnic 


poziomów menisków rtęci w manometrze. 


Na serię pomiarową składają się: 

- Wdestylowanie do naczynia pomiarowego znanej masy (ок. 16 ml obję- 
tości) czystego składnika i wykonanie pomiaru prężności pary w przewi- 
dywanym zakresie temperatur. 

- Kolejne dodestylowanie niewielkich, znanych ilości (rzędu 2-5 ml) 
czystego drugiego składnika i wykonanie pomiarów w funkcji temperatury 
dla każdego składu; dodestylowanie prowadzi sie do momentu przekrocze- 
nie objętości substancji dwukrotnie większej niż objętość uprzednio 
wprowadzonej substancji czystej. 

-~ Odcięcie naczynia równowagowego od aparatury na odcinku pomiędzy mano- 
metrem a chłodnicą. Naczynie było myte w kwasie azotowym, chromiance, 
a po wielokrotnym wypłukaniu wodą, płukane acetonem cz.d.a. i suszone 
strumieniem gazu. Ponieważ czyszczenie naczynia jest długotrwałe, po- 
sługiwano się dwoma bliźniaczymi egzemplarzami. 

- Opróżnienie aparatu z rtęci i wymiana jej na czystą; oczyszczanie rtę- 
ci podane jest niżej. 

Przyspawanie czystego naczynia i napełnienie rezerwuaru czystą rtęcią 


czyni aparat gotowy do następnej serii pomiarowej. 


VI.4.2. Algorytm obliczania prężności pary z danych pomiarowych 


Typowy pomiar prężności pary polega па zmierzeniu przy pomocy Ка- 
tetometru różnicy wysokości menisków rtęci w manometrze wewnętrznym (lub 
w obu manometrach - jeśli manometr zewnętrzny jest stosowany), poziom 


znacznika kalibrowanej części manometru wewnętrznego oraz poziom cieczy 
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w naczyniu równowagowym. И przypadku stosowania manometru zewnetrznego 
należy odczytać aktualną temperaturę w termostacie powietrznym. Podane 
wielkości plus temperatura w termostacie właściwym pozwalają na oblicze- 


nie równowagi ciecz-para panującej w układzie. 


Ponieważ znana jest całkowita objętość x naczynia równowagowego 
(patrz pkt VI-4.1.) objętość fazy parowej obliczamy przez odjęcie od 
niej objętości fazy ciekłej V， w temperaturze pomiaru obliczanej w/g 


WZOTU: 
у, = Y п,/р, (T) VI-2 


gdzie m, jest aktualną masa składnika i w naczyniu, a p (T) - jego ges- 


tošcia wyrazona w postaci równania p > а + b T. Wpiyw objetošci nadmia- 


гоме} jest pomijalny ze względu па to, Hi stanowi przewaznie mniej niz 
10 4 objetosci molowej mieszaniny, a ponad to pracujemy z utamkami moli. 
Objętość mieszania jest więc mniejsza niż błąd określenia objętości na- 
czynia. Do uzyskanej wartości należy dodać objętość ramienia manometru 


ponad rtęcią, którą oblicza się jako objętość walca. 


Ciśnienie panujące nad roztworem jest wynikiem sumowania sie posz- 


czególnych udziałów ciśnień przedstawionych na rys. VI-7. 


pi = p, * Р, + p; = Ds = Ds VI-3 
gdzie: 
h,p,g 
р = 2 VI-4 
Pna273 8, 


h, jest wysokością słupa medium wywierajacego ciśnienie; Pi jest јево 
gęstością w temperaturze T1 lub T2, ата gesto$cia rteci w 273 K, g 
jest lokalnym, a S srednim przyspieszeniem ziemskim rownym odpowiednio 
9.802 i 9.81 m sec 2. Gęstość rtęci w funkcji temperatury przedstawia 


WZÓT: 


= 1 VI-5 
Рн (1 + 0.0001818-t)0. 07355597 


który jest modyfikacja wzoru podanego przez Beattie i wsp. (1941) dla 
objętości rtęci w funkcji temperatury. 


Znaczenie poszczególnych udziałów jest następujące: 
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p. - ciénienie stupa rteci w manometrze zewnetrznym w temperaturze Т1 
jeśli stosowany jest manometr zewnętrzny; 

p. — ciśnienie słupa rtęci w manometrze wewnętrznym w temperaturze T2; 

p. - ciśnienie wywierane przez słup pary mieszaniny znajdującej się w 
naczyniu równowagowym; 

р - ciśnienie słupa azotu pod ciśnieniem P, i w temperaturze T1 jeśli 
stosowany jest manometr zewnętrzny; 

p. - ciśnienie słupa azotu pod ciśnieniem P, i w temperaturze T2 jeśli 


stosowany jest manometr zewnętrzny. 


do próżni 


< 


Rys.VI-7. Rozkład ciśnień w układzie pomiarowyn. 


Następnie, w procesie iteracyjnym oblicza się skład fazy parowej 
stosując wzór na ciśnienie całkowite nad roztworem i wzór na 
współczynnik lotności fazy parowej. Niedoskonałość fazy parowej jest га- 
wsze uwzględniana, a współczynniki wirialne estymowane są metodą Tsono- 
pulosa (1974) lub O'Connella i Prausnitza (1967) w zależności od badane- 


go układu. 


Osobnego omówienia wymaga sposób obliczania współczynników aktyw- 
ności składników w fazie ciekłej. Ponieważ wyznacza się je tylko z jed- 
nego punktu można zastosować jednoparametrowy model na nadmiarową епег- 


gię Gibbsa (GP). W algorytmie zastosowano dwa modele: roztworów regular- 
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nych (omówiony wczesniej) oraz procedure zaproponowana przez Bruina i 
Prausnitza (1971) zastosowaną do równania NRTL, która wymaga dobierania 
tylko jednego parametru energetycznego. Zgodnie z tą procedurą rezydual- 


na energia Gibbsa E jest proporcjonalna do entalpii parowania E, 
g, = -ВЕ, VI-6 
Entalpie parowania oblicza sie z równania Clausiusa-Clapeyrona: 


VI-7 


E = -z, RT [re 


| 
i 


dT 
Paremetr ściśliwości z, oblicza się z wirialnego równania stanu (Praus- 
nitz i wsp., 1967). 
s 
s 1 1 By” i 


^К W ko 


Pochodną logarytmu ciśnienia po temperaturze można obliczyć z równania 
wyrażającego jego zależność od temperatury, jak np. Antoine'a. Współ- 
czynnik proporcjonalności B, dobrany empirycznie metodą prób i błędów, 
posiada wartość ok. 0.2. 


Kolejnym etapem jest obliczenie poprawek związanych z wywieraniem 
ciśnienia przez słup fazy parowej z jednej strony oraz, ewentualnie w 
przypadku stosowania manometru zewnętrznego, słupa azotu z drugiej. Ро- 


prawki te są przedstawione schematycznie na rysunku VI-7. 


Ponieważ zarówno w przypadku pomiarów za pomocą jednego manometru, 
jak i dwóch - jedno z ramion dołączone jest do próżni nie gorszej niż 
0.7 Pa. Tak więc nie jest konieczne uwzględnianie poprawki na prężność 
pary rtęci, która i tak jest stosunkowo mała i dla temperatury 353 K му- 
nosi ok. 1 Pa (Brombacher i wsp., 1960). 


Ostatnim etapem każdej iteracji jest znormalizowanie obliczonego 
ciśnienia, wyrażonego w milimetrach słupa rtęci, do warunków normalnych 
tzn. 273 K poprzez pomnożenie przez stosunek gęstości rtęci w warunkach 
pomiaru do gęstości w temperaturze 273 K. Oczywiście, dla pomiarów dwoma 
manometrami, część słupa rtęci przypadającą na manometr zewnętrzny redu- 


kowano względem temperatury w termostacie powietrznym. 
Obliczony współczynnik aktywności składników jest modyfikowany w 
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każdej iteracji. Gdy liczba moli każdego ze składników w fazie ciekłej 


nie zmienia się więcej niż o 1076 uznaje się, że dystrybucja składników, 


a zatem i składy faz zostały ustalone. 


T 
WCZYTANIE DANYCH O CZYSTYCH 
SKLADNIKACH 


WCZYTANIE DANYCH POMIAROWYCH 
DLA POJEDYNCZEGO POMIARU 
Nie јслу ZEWNĘTRZNY MANOM? Iak 


Oblicz: ilosci то11, 
cisnienie P2: 


Oblicz.gest.N, 
W temp.T i Т, 


składy faz ciekłej i 
parowej dla 7 =1,ф,=1 


Oblicz ciśn.p,, 
poprawki РА Ре 


Oblicz 11056 moli w 
fazie рагомеј w obj.V 


Oblicz wsp.aktywn.w/g 


1 рорг. p 
3 
jednoparametr.modelu 
P - Pob1<10. Oblicz nowy sktad 
fazy ciekłej 
Oblicz skład fazy par. 


i wsp. lotności skł. 


Oblicz ilość moli w fazie 


parowej i odejmij od ilości 
moli w fazie ciekłej 


Zapamiętaj ilość moli Tak 
w fazie ciekłej 


CZY ILOŚĆ 


MOLI KTÓREGOŚ SKŁ. [Nie |CZY OSTATNI| Nie 
ZMIENIŁA SIĘ O PUNKT 


WIĘCEJ NIŻ 
10 5 Tak 


Rys. VI-8. Schemat obliczeń składu fazy ciekłej z danych doświadczalnych. 
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Wszystkie obliczenia pomiarów w tej pracy były wykonane w/g opisa- 


nego algorytmu opracowanym programem EXPVLE2. 


VI.5. POMIARY TESTOWE 


Nowozbudowany przyrząd wymaga wszechstronnego przetestowania. Nale- 
ży wybrać taki układ, aby wykazał on wszelkie słabe strony aparatu (w 
celu ich ewentualnego usunięcia), zweryfikował projektowany zakres po- 


miarowy oraz pozwolił ustalić metodykę pomiaru. 


VI.5.1. Wybór układu testowego 


Istnieje bogata literatura na temat pomiarów równowag ciecz-para 
omawiająca różnego typu przyrządy, metody pomiarów sposoby obliczeń i 
korelacji i.t.d. Brak natomiast informacji na temat zalecanych sposobów 
weryfikacji czy testowania budowanych urządzeń, badania ich dokładności 
i precyzji. O ile dla pomiarów np. objętości mieszania czy entalpii mie- 
szania są zalecane przez IUPAC układy do testowania, w pomiarach równo- 
wag ciecz-para istnieje zupełna dowolność w wyborze takiego układu. О- 
czywiście wybór układu testującego w głównej mierze zależy od przewidy- 
wanego zakresu ciśnień i temperatur, zastosowanych materiałów i.t.p. Nie 
mniej, istnieje konieczność wyboru takiego układu lub układów, aby można 
było porównać uzyskane wyniki z wynikami innych autorów. Taką pracą, 
która w pewnej mierze wypełnia tę lukę jest seria artykułów Smitha i 
współpracowników (Smith et al, 1982a,b) na temat układu (benzen(1) + cy- 
kloheksan(2)). Istnieje kilkadziesiąt prac podających rezultaty pomiarów 
równowag ciecz-para dla tego układu, wykonanych różnymi metodami, lecz 
autorzy wybrali tylko kilka jako naprawdę godnych uwagi. W tej pracy 
właśnie ten układ został użyty do testowania aparatu z następujących po- 
wodów: 

- łatwa dostępność składników, 

- względnie duża liczba opublikowanych danych innych autorów, 

- odpowiedni zakres prężości par, pokrywający przewidywany zakres pomia- 
rowy, 

- silnej wrażliwości na zawodnienie. 


Ostatni punkt jest szczególnie ważny ze względu na to, że jako medium 
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termostatujace używana jest woda. Układ jest w stanie wykazać ewentualne 
nieszczelności układu. Było to szczególnie widoczne w czasie prób skon- 


struowania systemu ciągłego dozowania składników. 


VI.5.2. Przygotowanie pomiarów testowych 


Oba składniki układu testowego zostały użyte bez specjalnego o- 
czyszczania. Zarówno benzen (POCh, cz.d.a.) jak i cykloheksan (Chemipan, 
standard chromatograficzny) zostały przedestylowane pod normalnym ciś- 
nieniem do jednego ze zbiorników naczynia do odgazowania (rys.VI-5. ), 
które podłączone zostało do kanału próżniowego. Procedura odgazowania 
polegała na wielokrotnym (10 - 15 razy) odciaganiu próżni znad ogrzewa- 
nej w jednym ramieniu zbiornika substancji; jednocześnie destylowała ona 
do drugiego ramienia. Już dziesięciokrotne powtórzenie tej operacji da- 
wało substancję wolną od rozpuszczonego powietrza, co zostało stwierdzo- 
ne pomiarami prężności pary w temperaturze 298 K. Dla pewności operację 


powtórzono jeszcze kilkakrotnie w obu przypadkach. 


Następnie do aparatury podłączano (do wejścia K7 lub K8), zważone z 
dokładnością 0.1 mg, ampułki (A, rys.VI-6.), które odpompowywano na 
gorąco przez co najmniej godzinę, po czym, przy zamkniętych kranach do 
próżni i innych części układu, wdestylowano substancję do ampułki ze 
zbiornika podłączonego do drugiego wejścia, wymrażając ją ciekłym azo- 
tem. Po odcięciu palnikiem, ogrzaniu do temperatury pokojowej i umyciu 
szlifu ampułki ponownie ważono określając masę substancji. Szacowany 
błąd tej metody przygotowywania próbek wynosi = 1 mg co stanowi poniżej 
0.1 % masy próbki. 


W ten sposób przygotowano, dla każdego ze składników, po 3 ampułki 
wyjściowe o pojemności ок. 16 ml, oraz 18 - 20 ampułek o pojemności ok. 
5 ml. Masy substancji zawierały się: dla dużych ampułek do 12 g, dla 
małych zaś nie mniej niż 0.5 g. 


VI.5.3. Wyniki testowania przyrządu 


Pomiary testowe przeprowadzono zgodnie z procedurą opisaną w roz- 


dziale IV-4.1. Część ampułek została zużyta na pomiary wstępne i ustala- 
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nie ргоседигу postepowania. W rezultacie przeprowadzono pomiary w catym, 
założonym (tzn. od 293 do 343 K) zakresie temperatur dla czystego benze- 
nu i cykloheksanu oraz 13 sktadów syntetycznych. Wykonano pomiary w 11 
temperaturach co 5 K po trzy pomiary w kazdej temperaturze. Obliczenia 
preznošci par substancji czystych i ich mieszanin z odczytów pomiarowych 
wykonano w/g algorytmu podanego w rozdziale VI-4.2. Jako model fazy cie- 
kłej zastosowano jednoparametrowy model roztworów regularnych. Założono 
niedoskonałość fazy parowej. Współczynniki wirialne obliczono metodą 
Tsonopulosa (1977).Wartości stałych i wielkości termodynamicznych ko- 
niecznych do tych obliczeń zestawiono w tabeli VI-2. 


Gęstość czystych składników obliczana była ze wzoru: 
p, = a, + b,*(T = 273.15) + c, (Т - 273.15)? VI-9 


Prężności pary czystych składników obliczne były równaniem Antoine'a w 
wersji SI: 
InP = А - B/(T-C) VI-10a 
i i i i 
lub oryginalnej: 
log Р, = А, 一 B /(C + t) VI-10b 


Tablica VI-2. 


Wartości stałych termodynamicznych i współczynniki równań opisujących 
wielkości termodynamiczne dla benzenu i cykloheksanu. 


Wielkość Cykloheksan 


Masa cząsteczkowa У 84. 156 ° 
Temperatura wrzenia 80.738 [ C] 
Temperatura krytyczna 554.1 [K] 
Ciénienie krytyczne š 40.6 [atm] 
Objętość krytyczna 309. 6 [ст] 
Wspólczynnik acentryczny 0.208 
Współcz. acentr. homom. 0. 208 
Moment dipolowy à 0.0 [D] 
Stałe równania na gęstość 
fazy ciekłej 
a 0. 90277 EN 0. 7974 > 
wzór VI-9 b -1.0011*10 ` -9.413•10 
C 0. 0. 
State równania Antoine'a 
6. 905660 6. 844980 
wzór VI-10b 1211. 033 1203. 526 
220. 790 222.863 
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Benzen (done wlosne) 
Boublik 1 wsp. 1973) 


** = - Benzen 
65.00 


Preznosc pory (kPo) 
№ ~ 

a л 

e e 

o o 


5.00 
290.00 300.00 310.00 320.00 330.00 340.00 350.00 


Temperatura (K) 


Rys.VI-9. Porównanie preznosci pary czystego benzenu z danymi literatu- 
rowymi. 


Cykloheksan (dane wlasne) 
s.s... Cykloheksan (Boublik i wsp. 1973) 


65.00 


45.00 


25.00 


Preznosc pary (kPa) 


5.00 
290.00 30000 310.00 320.00 330.00 340.00 350.00 


Temperatura (K) 


Rys.VI-10. Porównanie prężności pary czystego cykloheksanu z danymi 11- 
teraturowymi. 
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Ма rysunkach VI-9. i VI-10. porownano uzyskane preznosci рагу ben- 
zenu i cykloheksanu 2 danymi przedstawionymi w monografii Boublika i 
wsp. (1973) zawierającej rekomendowane prężności pary substancji czys- 


tych. 


Wyniki pomiarów prężności pary czystego benzenu i cykloheksanu jak 
i poszczególnych mieszanin przedstawione są w załączniku D. W nazwie 
próbki zawarty jest skład syntetyczny, który niekoniecznie jest równy 
rzeczywistemu w naczyniu równowagowym. W przypadku mieszanin, prężności 
pary mieszanin obliczane były z własnych stałych równania Antoine'a 


również zamieszczonych w tablicy D-1. 


Tablica VI-3. 


Temperatura: 313.136 K 


Współczynniki o 
C1= .442021 . 618828Е-03 
C2= -.242866E-01 .298491E-02 RMSD=0. 004 kPa 
C3= .391393E-02 . 320889E-02 
C4= -.455298E-03 . 806661E-02 


obl. Pexp AP 
[kPa] [kPa] 


кә кэ к» HB кэ кэ вэ ка ка ка Emp рэ на 
O > O O i> O O омон > 


0. 

0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
0. 0 
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Tablica VI-4. 


Temperatura: 333. 136 K 


Wspótczynniki 
С1= .386014 
C2= -.138488Е-01 
C3= -.535412E-03 
C4= -.282884Е-02 


ођ1. Рехр 
[кРа] 


ооооооооооооо 
ооооооооооооо 


Tablica VI-5. 


с 


.178654Е-03 
. 554046E-03 
. 905222E-03 
. 147168Е-02 


AP 
[kPa] 


Temperatura: 343.122 K 


Współczynniki 
C1= .367985 
C2= -.815266E-02 
C3= -. 334635E-02 
C4= -.309402E-02 


Pexp 
[kPa] 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


с 


. 195154E-03 
. 79407 2Е-03 
. 994417E-03 
. 210610E-02 


AP 
[kPa] 


КМ5р=0. 007 kPa 


GE 


7, 75 [J/mo1] 


ке кә кэ кэ кэ кэ кэ mB кэ кэ кэ кэ р> 
ч лоб Q ‹л&О Q O O (л(л оо 


RMSD=0. 010 kPa 


GE 


7, Ta [J/mol] 


к> кэ кэ кэ ка вэ вэ кэ кэ ка кэ вэ ә 
O O O. JE i> > = = > O. (О (л 
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2 26.00 
= 
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0 ***** Dane Scatcharda (Smith i wsp. 1982b) 
G 24,00 Dane wlasne (311.136 K) 

23.00 
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x(1) 


Rys.VI-11. Porównanie zmierzonej całkowitej prężności pary układu {Ъеп- 
zen(1) + cykloheksan(2)) z danymi literaturowymi. 


58.00 


57.00 . 


сл 
e 
o 
o 


55.00 + 


Cisnienie calkowite (kPa) 
r» 
© 


53.00 
52.00 
*+++* Done Boublika (Smith i wsp. 1982b) 
Dane wlasne (333.140 K) 
51.00 
50.00 
0.00 0.20 0.40 0.60 0.80 1.00 
x(1) 


Rys.VI-12. Porównanie zmierzonej całkowitej prężności pary układu {Ъеп- 
zen(1) + cykloheksan(2)) z danymi literaturowymi. 
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Rys.VI-13. Porównanie zmierzonej całkowitej prężności pary układu {Ъеп- 
zen(1) + cykloheksan(2)) z danymi literaturowymi. 


Na rysunkach od VI-11 do VI-13 porównane są izotermy interpolowane 
z danych zawartych w załączniku D z najlepszymi danymi literaturowymi 
rekomendowanymi przez Smith'a i wsp. (1982b). W tablicach od VI-3. do 
VI-5. podane są wartości numeryczne uzyskane z korelacji izoterm równa- 
niem Redlicha-Kistera, stałe tego równania i ich wariancję ст, oraz 
średnie odchylenie standardowe RMSD (wzór VII-11 w rozdz.VII.2). W ta- 
blicach przedstawiono również obliczone składy fazy parowej, współczyn- 
niki aktywności obu składników i nadmiar entalpii swobodnej Gibbsa. AP 
jest różnicą pomiędzy ciśnieniem eksperymentalnym a obliczonym. 


VI.5.4. Omówienie wyników testowych 


Uzyskane wyniki potwierdzają przydatność zbudowanej aparatury do 
pomiarów równowagi ciecz-para w przewidywanym zakresie temperatur tzn. 
od 293 do 343 K. Potwierdzony został założony zakres ciśnień, który poz- 
wala mierzyć aż do ciśnienia zbliżonego do atmosferycznego. Limit ten 
wynika z zastosowanych kranów próżniowych, które mogą poprawnie pracować 


pod minimalnym, ale koniecznym, podciśnieniem. Niewielkie rozbieżności z 
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danymi literaturowymi w przypadku mieszanin nie maja charakteru systema- 
tycznego i mogą wynikać z minimalnych zanieczyszczeń zawartych w benze- 
nie, który nie był badany na zawartość zanieczyszczeń. Destylacja prze- 
prowadzona została w celu usunięcia śladów wilgoci zawartych w benzenie. 
Cykloheksan jako źródło zanieczyszczeń należy raczej odrzucić ze względu 


na użycie standardu chromatograficznego. 


VI.6. POMIARY WŁAŚCIWE 
VI.6.1. Przygotowanie odczynników 


3-MOPN. Związek został uzyskany na drodze syntezy (Astle i Ethe- 
rington, 1952) z metanolu i akrylonitrylu w środowisku alkalicznym wg 


poniższej reakcji: 


CHAOH * CH,,-CH-CN e CH40-CR, -CH, -CN 


Wydajność syntezy wynosi ok 85 %. Produkt oczyszczony został droga des- 
tylacji frakcyjnej. Zbierano frakcję wrącą w temperaturze 435 + 0.1 K. 
Czystość sprawdzono metodą GLC. Szczegóły oczyszczania podane są w pracy 
doktorskiej M. Antosik (1991) oraz w pracy na temat pomiaru nadmiaru 


objętości (Łuszczyk, 1988). 


Woda. Wobec silnego zanieczyszczenia wody dostępnej w sieci miejs- 
kiej przetestowano skuteczność oczyszczania wody na trzech jej próbkach. 
Pobrano wodę z kranu, ze studni indywidualnej w Rembertowie, oraz z sie- 
ci wodociągowej na osiedlu w sąsiedztwie Akademii Stabu Generalnego w 
Rembertowie. Oczyszczanie przeprowadzono w identyczny sposób, а miano- 
wicie: 

- wstępna destylacja; 

- gotowanie przez 3 godz. z dodatkiem nadmanganiamu potasowego w 
środowisku zasadowym i następnie przedestylowanie; 

- destylacja ponowna po uprzednim zakwaszeniu kwasem siarkowym. 

Tak uzyskana woda została zweryfikowana na zawartość związków aromatycz- 

nych poprzez pomiar absorpcji w widmie promieniowania ultrafioletowego i 

widzialnego za pomocą spektrofotometru firmy Beckmann model III. Naj- 


mniejszą absorpcję wykazywała próbka pobrana na osiedlu w Rembertowie. 
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Та wiec мода zostata uzyta do dalszych pomiarów. Dalsza weryfikacja ро- 
legata na pomiarze gestošci w funkcji temperatury (Luszczyk, 1988) oraz 


pomiaru temperatury wrzenia metoda ebuliometryczna (Gregorowicz, 1986). 


Rteé. Do pomiarów stosowana była rtęć cz.d.a (P.O.Ch. - Gliwice). 
Po każdzej serii pomiarowej rtęć była oczyszczana w/g procedury podane j 
przez Luxa (1960) tzn. była przepuszczana przez 8 % kwas azotowy, płuka- 
na w wodzie destylowanej i suszona w acetonie. Dodatkowo, przed użyciem 


była destylowana pod zmniejszonym ciśnieniem. 


VI.6.2. Stałość i stabilność temperatury 


We wstępie do tej pracy powiedziano, że dokładność i stabilność 
temperatury ma kluczowe znaczenie w badaniu punktu krytycznego. Wszyst- 
kie pomiary wykonane w ramach tej pracy prowadzone były pod kontrolą re- 
gulatorów temperatury UNIPAN typ 650 z przystawką mocy typ 651. Zapew- 
niało to stabilność temperatury rzędu 1 mK w bezpośrednim otoczeniu na- 
czynia pomiarowego w czasie trwania całego pomiaru w danej temperaturze 
tzn. przez 2-3 godz. Sprawdzane to było wielokrotnie w następujący 
sposób: pomiar temperatury wykonywany był termometrem kwarcowym firmy 
Hewlett-Packard model 2801A. Rozdzielczość termometru wynosi 0.1 mK. Do 
termometru podłączona była drukarka ERD-100 firmy MERA-KFAP poprzez in- 
terfejs zbudowany w Zakładzie Elektroniki IChF PAN. Podana konfiguracja 
pozwalała zapisywać temepraturę w zadanym interwale czasu. Najczęściej 
temperaturę zapisywano со 30 sek. co dawało ok. 400 odczytów dla testo- 
wanej izotermy. Z tej liczby pomiarów wyciągano wartość średnią oraz od- 
powiednią wariancję. Dla temperatur nieco wyższych od pokojowej warian- 
cja wynosiła 0.7 mK natomiast dla temperatur rzędu 350 K - 1.6 mK. War- 
tość bezwzględna temperatury określana była za pomocą okresowego skalo- 
wania termometru względem termometru oporowego firmy Tinsley znajduja- 
cego się w Zakładzie III IChF w zakresie temperatur 232 - 353 K. Odczyt 
z termometru przeliczany był do wartości absolutnej w stopniach Kelvina 


za pomocą wzoru: 


T = 273.15 + (0.248 + 0.999 tp) VI-11 


gdzie t jest temperatura odczytana z termometru kwarcowego w stopniach 


Celsjusza. 
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VI.6.3. Stan równowagi termodynamiczne j 


Istotnym zagadnieniem w pomiarach równowag ciecz-para jest to, czy 
układ jest w stanie równowagi termodynamicznej. Wtedy, gdy nie prowadzi 
się analizy składów obu faz, jedynym kryterium osiągnięcia stanu 
równowagi jest niezmienność ciśnienia w funkcji czasu. Jak wiadomo, 
układ trójfazowy szczególnie trudno dochodzi do stanu równowagi. W celu 
określenia czasu dojścia wykonano kilka serii pomiarowych 12 i 24 
godzinnych, w których mierzono ciśnienie panujące w układzie co 30 min. 
Uzyskane wyniki wskazują, że w przypadku układu trójfazowego stałość 
ciśnienia osiągana była, przy stosowanym systemie mieszania faz ciekłych 
i fazy parowej, już po ok 2 godz. Wszystkie pomiary ciśnienia wykonywane 
więc były po upływie tego czasu od osiągnięcia przez układ żądanej tem- 


peratury. 


VI.6.4. Równowaga ciecz-para układu (3-MOPN(1) + woda(2)) 


Przygotowano zestaw ampułek z czystymi substancjami pozwalające u- 
zyskać 31 składów w całym zakresie stężeń (Tablica VI-7 i VI-8). Nies- 
tety, wobec dużej różnicy w masach cząsteczkowych trudno było uzyskać 
składy o niewielkiej zawartości wody bez utraty dokładności określenia 
masy wody w ampułce. Z drugiej strony, ten zakres stężeń nie był intere- 
sujący ze względu na jego homogeniczność w badanym zakresie temperatur i 


znaczne oddalenie od składu krytycznego. 


Konieczne do korelacji współczynniki wirialne obliczone zostały 
metodą O'Conell'a i Prausnitza (1967) jako bardziej odpowiednią dla 
substancji silnie polarnych. Wymagane do tego celu parametry krytyczne 
3-MOPN zostały obliczone programem CPEST opracowanym na podstawie pracy 
Somayajulu (1989). Jest to metoda udziałow grupowych, której średni błąd 
w określaniu parametrów krytycznych jest mniejszy niż 2 Ж. Stałe wymaga- 
ne w korelacjach dla obu składników zebrane są w tablicy VI-6. 
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VI.6.5. Rezultaty pomiarów 


Prężności рагу 3-MOPN były już uprzednio badane przez Rogalskiego 
metodą ebuliometryczną (Stryjek et al,1978b) w zakresie temperatur od 
323 K do temperatury wrzenia oraz przez Janaszewskiego dla temperatur 
293.15, 313.15 i 333.15 K (Stryjek i wsp. 1978c) metodą statyczną. МуКо- 
nano trzy niezależne pomiary, z których dla dwóch wartości liczbowe 
przedstawione są w załączniku E. Miało to na celu porównanie wzajemnej 
konsystencji wyników oraz porównanie ich z dostępnymi danymi. Uzyskane 
wyniki wskazują na bardzo dobrą zgodność  (rys.VI-14.) z wynikami 


pomiarów poprzednich. 


Tablica VI-6. 


Wartości stałych molekularnych i innych wielkości wymaganych w korela- 
cjach 


Masa cząsteczkowa š £ 18.015 ^ 
Temperatura wrzenia а 100.0 [ C] 
Temperatura krytyczna 3 647.4 [K] 
Ciénienie krytyczne š m 218.3 [atm] 
Objętość krytyczna š 55.2 [cm] 
Współczynnik acentryczny : 0.344 
Współcz. acentr. homom. Я 0.01 
Moment dipolowy 1.87 [D] 
Stałe równania na gęstość 
fazy ciekłej 
a 0. 95848 1. 000295 
wzór VI-9 b -9.543:10 * 3. 396-10 Š 
c 0. -5.018-10 Š 
State równania Antoine'a 
8.510396 8.029549 
wzór VI-10b 2425. 490 1708. 778 
270.274 231. 824 


Pomiary preznosci рагу wody wykonano w celu porownania z danymi 11- 
teraturowymi (sprawdzenie czysto$ci) i wykonano tylko jedna serie pomia- 
rowa. Rezultat prezentuje rysunek VI-15. Jako dane porównawcze przedsta- 
wiono dane ТВС (1973). 
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Rys.VI-14. Prężność pary czystego 3-MOPN w badanym zakresie temperatur. 
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Rys.VI-15. Porównanie prezno$ci рагу czystej моду z danymi literaturowy- 
mi w badanym zakresie temperatur 
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Rys.VI-16. Porównanie rezultatów własnych pomiaru УГЕ z danymi literatu- 
rowymi. 


Rysunek VI-16. przedstawia rezultat interpolacji dla temperatury 
323.270 K oraz dla porównania wyniki Rogalskiego (Stryjek i wsp. 1978b) 
dla tej samej temperatury. Widać bardzo dobrą zgodność uzyskanych 
wyników (dane numeryczne w tablicy VI-8). Należy podkreślić, że poprzed- 
nie wyniki uzyskane zostały metodą ebuliometryczną, a więc różną od sto- 
sowanej w tej pracy. Ebuliometr pracował w warunkach izotermicznym, tzn. 
ciśnienie wymuszające było tak dobierane, aby utrzymać stałość tempera- 
tury. 


W przypadku tego układu przeanalizowano wpływ użytego modelu do 
obliczania współczynników aktywności na skład faz ciekłej i parowej. Jak 
widać na rysunkach VI-17 i VI-18 wpływ doboru modelu jest praktycznie 
bez znaczenia dla składu fazy ciekłej, natomiast dla fazy parowej nie 
może być pominięty. Rysunek został zrobiony na podstawie obliczeń dla 
próbki nr 13, oznaczonej w tablicy VI-7 jako "mw2334", a więc dla obsza- 
ru najbardziej interesującego. Różnica w składach fazy parowej jest 
przyczyną nie zamieszczenia ich w tablicy VI-8 oraz w tablicy wyników w 
załączniku E. 


Prężności pary składników czystych jak i mieszanin (3-MOPN(1) + wo- 
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da(2)) nie wymagaty stosowania manometru zewnetrznego. Рефпе rezultaty 
pomiarów zawierajace temperatury i odpowiadajace im ci$nienia dla posz- 
czególnych próbek i składników czystych, wraz z odchyleniami od równania 
Antoine'a (wzór VI-10a) podane są w załączniku F. Otrzymane, w wyniku 
przetworzenia programem EXPVLE2, stałe równania Antoine'a przedstawione 


są w Tablicy VI-7. 
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Rys.VI-17. Porównanie obliczonych składów fazy ciekłej przy użyciu róż- 
nych modeli. 


Zarówno w tablicy VI-7 (kolumna 2) jak i w nagłówkach tablic za- 
łączników D i E, dla rozróżnienia próbek, stosowane są nazwy 
zawierające, obok pierwszych liter nazw tworzących je składników skład 
wyjściowy mieszanian wynikający z naważek czystych składników. Często 
różni się on od składu mieszaniny w naczyniu pomiarowym (szczególnie dla 


roztworów o dużej zawartości 3-MOPN). 
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Rys.VI-18. Porównanie obliczonych sktadów fazy parowej przy uzyciu róz- 
nych modeli. 


Dane uzyskane nie stanowią jeszcze danych opisujących równowagi 
ciecz-para. Wobec względnie dużej objętości fazy parowej znacząco zmie- 
nia się skład fazy ciekłej w stosunku do składu syntetycznego. Jest to 
wynikiem dużej różnicy w prężnościach czystych składników. Jak wspomnia- 
no, w procesie numerycznej obróbki danych, składy fazy ciekłej były 
obliczane dla każdej temperatury i przybliżane za pomocą wielomianu dru- 
giego stopnia w funkcji temperatury. Posiadając funkcyjną zależność 
prężności pary i składu fazy ciekłej można przeprowadzić proces interpo- 
lacji tychże dla poszczególnych temperatur. Wybrano 16 temperatur w oto- 
czeniu punktu krytycznego, takich samych jakie badane były w pracy na 
temat objętości nadmiarowej dla tego układu (Łuszczyk, 1988). Wyniki in- 
terpolacji stanowią bazę do prowadzenia dalszej obróbki numerycznej w 
celu uzykania danych o równowadze ciecz-ciecz-para. Rezultaty w postaci 
składu fazy ciekłej і odpowiadającemu mu ciśnieniu dla poszczególnych 
temperatur przedstawiono w Tablicy VI-8. Składy i ciśnienia wydrukowane 


kursywą mierzone były nad roztworem o dwóch fazach ciekłych 
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Tablica VI-7. 


współczynniki równania VI-10a opisujące prężności substancji czystych i 
poszczególnych próbek. W tablicy podane są również odchylenia standardo- 
we. NP - liczba punktów doświadczalnych; NU - liczba punktów użytych w 
korelacji. 


Współczynniki równania 
A 


woda . 473733 
mw0191 . 718708 
mw0465 . 597942 
mw0551 . 834051 
mwO0801 . 920353 
mw0897 . 073603 
mw1104 .689529 
mw1417 . 851747 
mw1772 .065217 
mw1933 . 927257 
mw2139 . 175023 
mw2163 . 012472 
mw2334 . 845025 
mw2374 . 037566 
mw2399 . 844201 
mw2556 . 949682 
mw2656 . 790391 
mw2823 . 948707 
mw2987 . 837743 
mw3146 . 712344 
mw3322 . 572853 
mw3533 .572909 
mw3672 . 726138 
mw3994 . 244526 
mw4445 . 185658 
mw5008 ‚ 130613 
mw5707 . 440675 
mw6601 . 971051 
mw6890 . 893330 
mw7614 . 107835 
mw8420 . 708396 
mw8882 . 438238 
3-торп . 580926 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
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Tablica VI-8. 


Preznosci pary [kPa] układu (3-MOPN(1) + woda(2)) dla temperatur w oto- 
czeniu temperatury krytycznej mieszalności. 
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Tablica VI-8. c.d. 
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Tablica VI-8. c.d. 
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Tablica VI-8. с.а. 


312.177 K 313.160 K 318.150 K 323.270 K 
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VII KORELACJE 


Termin "korelacja" oznacza znalezienie Zaleznosci pomiedzy pewnymi 
wielkościami (nie koniecznie bezpośrednich) przy zastosowaniu wyrażenia 
matematycznego zawierającego współczynniki dobieralne zwane parametra- 
ші," Bardzo czesto parametry maja mniejsze lub wieksze uzasadnienie 
teoretyczne. W przypadku równowag сіес2-рага i  ciecz-ciecz takim 
wyrażeniem pozwalającym np. obliczyć skład fazy parowej, która jest 
wielkością skorelowaną, może być np. równanie stanu lub wzór na nadmia- 
rową energię Gibbsa. Zastosowanie równań stanu nie jest omawiane w tej 
pracy. Szczegółowy przegląd tego zagadnienia przedstawiony został w рга- 
cach Anderki (1989, 1990). Tutaj zostaną omówione korelacje za pomocą me- 
tody у-ф (jako głównie stosowanej w tej pracy) czyli algorytm Barkera i 
równania korelacyjne oraz związane z tym zagadnienia doboru funkcji ce- 
lu, metody optymalizacji, liczby stałych i.t.d. Problemy te są bardziej 
szczegółowo opisane przez Hofmana (1982) w jego pracy doktorskiej. W 
dalszej części rozdziału omówione są trzy programy komputerowe stosowane 


do obliczeń w tej pracy oraz rezultaty tych obliczeń. 


VII-1 METODA BARKERA 


W 1953 roku Barker zaproponował algorytm do korelacji równowag 
ciecz-para przy wykorzystaniu niepełnych danych równowagowych, czyli na 
podstawie znajomości ciśnienia P nad roztworem o składzie x w temperatu- 
rze T. Oczywiście chodzi tu o zbiór tych informacji dla kilku składów. 
Zastosowana została metoda najmniejszych kwadratów, a funkcją minimali- 


zowaną jest suma kwadratów odchyleń ciśnienia obliczonego od zmierzone- 


*Bardzo często korelacja jest traktowana na równi z wygładzaniem da- 
nych, które polega właściwie na dobraniu współczynników wyrażenia mate- 
matycznego w celu opisania rezultatów pomiaru w zakresie zmierzonej 


zmiennej niezależnej. 
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go. W tej metodzie wymagana jest znajomość wyrażenia na nadmiarową 


energię Gibbsa. 


W literaturze opisano wiele równań korelacyjnych па B Ich 
szczegółowe przeglądy, włączając występujące modyfikacje, były wielo- 
krotnie publikowane (Hala, 1977, Serensen, 1979). Większość tych równań 
ma obecnie marginalne znaczenie (np. Margulesa czy Van Laara). Znaczenia 
nabrały te, których parametrom można przypisać pewną zależność tempera- 
turową i których stosowalnść może być rozciągnięta na układy więcej niż 


dwuskładnikowe. 


Poniżej przedstawione są te równania, które stosowane były w tej 


pracy. 


VII-1.1. Równanie Redlicha-Kistera 


Rownanie Redlicha-Kistera (1948) jest zależnością typu wielomiano- 
wego, liniową w funkcji jego współczynników. Jest bardzo często stosowa- 
ne do opisu wielkości nadmiarowych w funkcji stężenia układów dwusktad- 


nikowych w postaci: 


X (x) = x(1-x) ) а, (2x-1)' VII-1 
1=0 


gdzie: x jest stężeniem w postaci ułamka molowego, wagowego czy też 


objętościowego. 


Powyższa postać zapewnia zerowanie się funkcji XE w warunkach brze- 
gowych, wymagane dla funkcji nadmiarowych. Jest to równanie bardzo elas- 
tyczne, pozwalające opisać nawet bardzo skomplikowane przebiegi funkcji 
nadmiarowych. Jednocześnie, stosowane bezkrytycznie, może prowadzić do 
absurdalnych wyników. Sygnalizowane to jest w wielu pracach np. Gierycza 
(1982). Nadmienić należy, że opis bezpośrednio mierzalnych funkcji nad- 
miarowych takich jak entalpia mieszania czy nadmiarowa objętość miesza- 
nia jest typowym przykładem wygładzania danych, a nie korelacji. Istnie- 
je wiele odmian tego równania stosowanych do konkretnych układów i kon- 
kretnych wielkości. Modyfikacje polegają na dokonaniu pewnych transfor- 
macji zmiennej niezależnej (stężenia) [np. Neau (1972) zmodyfikowała 


zmienną stężeniową dodając parametr zależny od rodzaju korelowanego 
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układu], lub wykorzystaniu w korelacji wyrazów zawierających określone 
wykładniki wielomianu (np. tylko parzyste czy tylko ułamkowe i ich wie- 


lokrotności). 


Rozszerzenie tego równania па układy  wielosktadnikowe jest 
niezmiernie trudne ze względu na brak teoretycznych podstaw wyznaczenia 
współczynników równania charakterystycznych dla układu dwuskładnikowego 
i ich wzajemnej zależności w układzie wieloskładnikowym. Posiadając 
jednak wystarczającą ilość danych w funkcji dwóch parametrów, np. 
stężenia i temperatury, można pokusić się o znalezienie zależności 
funkcyjnej współczynników od, powiedzmy, temperatury. Pozwoliłoby to na 
możliwość interpolacji funkcji, lecz wyciąganie szerszych wniosków jest 


nie uzasadnione. 


Najprostsza forma równania R-K, tzn. z jednym współczynnikiem, jest 
równoważna wyrażeniu na p" wyprowadzonemu z teorii roztworów regularnych 
Guggenheima (1953). Wpływ dalszych członów równania VII-1 na kształt i 
symetrię GË byt ostatnio przedstawiony w podreczniku Prausnitza i wsp. 
(1986). 


W przypadku zastosowania równania R-K do korelacji równowag ciecz- 
para w bezpośredniej korelacji nie stosuje się formy VII-1 lecz jej po- 
chodne względem składu gdzie ХЕ = GE/RT. 


VII-1.2. Równanie NRTL i LEMF 


Równanie NRTL (Non-Random Two Liquid) opublikowane zostało przez 
Renona i Prausnitza (1968) jako pochodne i niejako uogólnienie wcze$- 
niejszych równań takich jak Wilsona (1964), Blacka (1959). Krótko potem 
Renon opublikował książkę (Renon i wsp.; 1971), w której opisano stoso- 
wanie tego równania do układów dwuskładnikowych z rozszerzeniem na inne 
wielkości termodynamiczne, takie jak entalpia mieszania, równowaga 
ciecz-ciecz i graniczny współczynnik aktywności. Wprowadzono jednocze$- 
nie liniową temperaturową zależność stałych. 
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GE/RT = у Ed VII-2 
1 


üdeler G w жыра iv N е“ By” E, ET Su. "al 


#,, jest energia oddzialywania czasteczek i i j, natomiast ET jest pa- 
rametrem wyrażającym, wg. autorów, odwrotność liczby koordynacyjnej. Pa- 


rametry т i с są dobierane podczas korelacji danych eksperymentalnych. 


Równanie NRTL należy do grupy równań opartych na tzw. koncepcji 
stężeń lokalnych. Koncepcja ta polega na tym, że przyjmuje się, Ze praw- 
dopodobieństwo znalezienia cząsteczki A w bezpośrednim sąsiedztwie cen- 
tralnej cząsteczki B jest proporcjonalne do różnicy energii oddziaływa- 
nia ВВ i AB. Czynnikiem ważącym to prawdopodobieństwo jest czynnik 
Boltzmana ехр(-от). Takie założenie było w przeszłości krytykowane 
(Flemr, 1975) jako nie konsystentne termodynamicznie, nie mniej jednak 
równanie to zdobyło sobie ogromną popularność. Wynika ona z dwóch pods- 
tawowych zalet: 

- jest wielokrotnie różniczkowalne względem składu (czynnik exp), co 
przy ograniczonej liczbie dobieralnych parametrów pozwala przewidzieć 
rozpad roztworu na dwie fazy ciekłe; 

- daje możliwość przewidywania własności roztworów wieloskładnikowych na 


podstawie parametrów uzyskanych z danych dla roztworów dwuskładnikowych. 


Istnieje bogata literatura opisująca zastosowania równania NRTL jak 
również jego modyfikacji. Oprócz, niejako klasycznego, zastosowania do 
korelacji równowag ciecz-para, równanie to jest stosowane do korelacji i 
przewidywania równowag ciecz-ciecz w układach dwu- i trójskładnikowych. 
Zdolności korelacyjne były badane przez wielu autorów. Ponieważ stwier- 
dzono, że równanie to może mieć wiele rozwiązań, Mattelin i Verhoeye 
(1975) zbadali w jakim zakresie zbiory parametrów uzyskane z korelacji 
danych ШЕ dają rezultaty konsystentne termodynamicznie. Stwierdzili 
oni, że dla układu symetrycznego najwyższą wartość jaką może przybrać 
parametr ск aby równanie pozwalało na rozwarstwienie się roztworu jest 
0.427. Dla wyższych wartości niemożliwe jest spełnienie warunków stanu 
krytycznego (druga i trzecia pochodna funkcji c względem składu równa 
zero). Zależność kształtu krzywej współistnienia od wielkości parametru 
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a badana była również ostanio przez Bittricha i wsp. (1989). Według ich 
badań, użycie parametru «œ bliskiego -1 (sprowadza się to do równania 
znanego jako LEMF (Marina i Tassios, 1973)) daje lepszy opis krzywej bi- 
nodalnej w bezpośredniej bliskości punktu krytycznego. Niestety, autorzy 
nic nie mówią na temat temperaturowej zależności stałych energetycznych. 
Problemy związane z odpowiednią korelacją równowag ciecz-ciecz wynikają 
z braku teoretycznie uzasadnionego wyrażenia wiążącego te dwie wielkoś— 
ci (tzn. z gj) = f(T)). Zawarta w czynniku Boltzmana zależność tem- 
peraturowa nie jest wystarczająca, gdy chcemy korelować w większym za- 
kresie temperatur i sami autorzy równania zalecają stosowanie liniowej 
zależności stałych od temperatury. W przypadku równowag ciecz-para jest 
to wystarczające, lecz dane LLE korelują się źle. W literaturze podane 
są przykłady zastosowań różnych wyrażeń, m.in. wielomianu kwadratowego, 
za pomocą którego Ochi i wsp. (1990) opisali nawet luki zamknięte. 


Problemy związane z jednoczesną korelacją równowag ciecz-para i 
ciecz-ciecz oraz wspólnej korelacji granicznego współczynnika aktywności 
były rozpatrywane w naszej pracowni. Autor niniejszej pracy zaproponował 
wspólną korelację równaniem NRTL jednego z granicznych współczynników 
aktywności i danych LLE (Stryjek i Łuszczyk, 1980) w celu uzyskania da- 
nych LLE w zakresie stężeń homogenicznego roztworu. Pozwala to obliczyć 
trzy parametry NRTL. W późniejszej pracy (Stryjek i wsp., 1981) do tego 
celu sprawdzono również inne równania (LEMF, UNIQUAC i LCG (Vera i wsp., 
1977)). Zastosowanie współczynnika aktywności w granicznym rozcieńczeniu 
jest szczególnie atrakcyjne, ponieważ można go stosunkowo łatwo zmierzyć 
za pomocą chromatografii gazowej (Conder i Young, 1979). Ponadto, 
najczęściej pomiar dotyczy układu, którego składniki różnią się bardzo 
prężnościami pary i pomiar równowagi ciecz-para jest bardzo trudny. Tego 
typu korelacje prowadzone są również w innych ośrodkach naukowych (patrz 
np. wcześniejsze i obecne prace Kikica i wsp., 1983). Ogólnie, problemy 
współnej korelacji ШЕ i МЕ przedyskutowane zostały w pracy Malanow- 
skiego i wsp. (1981). 


Ostatnio De Pablo i Prausnitz (1989), w celu poprawienia opisu 
krzywej współistnienia rówaniem NRTL, zastosowali procedurę tranformacji 
składów za pomocą funkcji skalującej, analogiczną do zastosowanej przez 
Foxa (1979) do opisu krzywej współistnienia ciecz-gaz w układach јейпо- 
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sktadnikowych. 


VII-1.3. Równanie UNIQUAC 


Równanie UNIQUAC jest kontynuacja idei stezen lokalnych, z tym, ze 
całkowity nadmiar energii Gibbsa podzielony jest na dwa udziały: jeden 
związany z dystrybucją cząsteczek tworzących roztwór spowodowanej różni- 
cami w kształcie i wielkości cząsteczek, oraz drugi - związany z energe- 
tycznym oddziaływaniem pomiędzy elementami cząsteczek. Wprowadzony przez 
Abramsa i Prausnitza (1978) zyskał sobie sporą popularność jako model 
choć częściowo  podbudowany teoria molekularną і w którym część 


parametrów może być określona na podstawie budowy molekularnej. 


Смеа ac VII-3 
comb res 
E $, = 1 
Gy RT = ) x ln = * 5 } q,x,1n E VII-4 
E = = - 
B „ИВТ = | an ) өт) VII-5 
x r x q 
Ф = Ж. e, = DE VII-6 
n 
) ыы } х 39, 
ў j 
* n expI[-(u,.- u,,)/RT] VII-7 


Parametry Ө (średni ułamek powierzchniowy) i Ф (średni ułamek обје о5- 
ciowy) obliczane są na podstawie molekularnych parametrów Р, i q, окге5- 
lonych ze struktury czasteczek i promieni van der Waalsa podanych przez 
Bondiego (1968) dla poszczególnych grup atomów. Parametr z jest liczba 
koordynacyjna sieci, na której rozmieszczone sa poszczególne zreby ato- 
mowe czasteczek, przyjeta przez autorów w oryginalnej pracy równa 10. 
Równanie UNIQUAC posiada tylko dwa dobieralne parametry energetyczne, w 
związku z czym jego zdolności korelacyjne sa mniejsze niż równania NRTL. 
W celu usunięcia tej niedogodności powstało wiele modyfikacji (Anderson 
i Prausnitz, 1975; prace Nagaty i wsp., 1981а, 1981b; oraz Nagaty, 1982, 
1985). Skjold-Jergensen i wsp. (1980, 1982) podali wyrażenie na tempera- 


turową zależność parametru z. Również parametry energetyczne podawane są 
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јако liniowe, kwadratowe czy hiperboliczne funkcje temperatury. Szcze- 
gółowy przegląd zastosowań i modyfikacji tego rówania podane zostały w 


pracach Serensena i wsp. (1979a,b,c) i Rasmussena (1983). 


VII-2. FUNKCJA CELU I ALGORYTMY OBLICZEŃ 


Korelacje wyników doświadczalnych wymagają nie tylko wyboru 
równania korelacyjnego, lecz także wyboru funkcji celu oraz algorytmu 
(sposobu) korelacji. W obliczeniach związanych z doborem parametrów 
równania korelacyjnego stosuje się metodę najmniejszych kwadratów opartą 
na warunku: 

У. ох. x (A)) = min VII-8 

1 
gdzie X oznacza dowolna zmienna zalezna, natomiast А wektor parametrów. 
Metoda ta jest opisana w wielu podręcznikach (Strzałkowski i Śliżyński, 
1978). Warunek VII-8 jest podany w najprostszej formie, lecz w pewnych 
przypadkach, gdy korelowane są równocześnie różne wielkości, wprowadzane 
są dodatkowe składniki. Wtedy, gdy udziały poszczególnych wielkości nie 
są równocenne, пр. z powodu różnej dokładności, czy też wielkośći abso- 
lutnej, konieczne jest wprowadzenie ważenia punktów, czyli nadawanie im 
różnego ciężaru statystycznego. 


Wartość obliczona Кеше jest nie tylko funkcja zmiennej (zmiennych) 


niezaleznej, lecz WYRAZY A. Najczęściej warunek VII-8 jest 
nieliniowy względem parametrów, co oznacza, że obliczenia muszą być рго- 
wadzone w sposób iteracyjny, a w kolejnych iteracjach uzyskujemy coraz 
lepsze przybliżenie parametrów równania. Wiąże się to z zagadnieniem wy- 
boru startowych parametrów równania, które pozwolą uzyskać zbieżny ciąg 
rozwiązań. Powszechnie stosowana jest metoda Newtona-Raphsona (Ralston, 
1975) jako najprostsza, lecz również metoda Marquardt'a (1963), dająca, 
przy umiarkowanej ilości parametrów dobieralnych, ciąg zawsze zbieżny. W 
tej pracy właśnie ta metoda była stosowana we wszystkich rozwiązywanych 


problemach. 


Powyższe uwagi dotyczą sytuacji, gdy liczba punktów eksperymental- 
nych znakomicie przewyższa liczbę dobieralnych parametrów. W przypadku, 
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рау liczby te sa równe, do znajdowania parametrów stosowana jest metoda 
Newtona polegająca na rozwiązaniu układu równań liniowych VII-9 względem 
parametrów. Równania te są rozwinięciem w szereg, obciętym po pierwszym 


wyrazie, wokół prawdziwych rozwiązań. 


ЗЕ Е, AF 

Pa WA UA. | ве ж эзер ДА Żaku 
1 2 n 

8F ЈА AF, 
от ад ВА dł: waw шый. wi dia VII-9 
. .1 .2 n . 

ӘР, BF ӘР, 
P s - ӘА 8А, + BA, 8А, + + 8A SA, = 0 


VII-3. KORELACJE RÓWNOWAGI CIECZ-PARA 


Dane zamieszczone w tablicy VI-8 byty korelowane metodami: Tao 
(1961), MGC (Mixon i wsp., 1965)i Barkera (1953) z wykorzystaniem równañ 
Redlicha-Kistera, МЕТ. a takze LEMF.Y Wszystkie te metody zostały za- 
implementowane w programie FITVLE, który zapewniał jednakowe procedury 
minimalizacji, zunifikowany sposób obliczania wielkości pomocniczych ta- 
kich jak np. prężności pary, współczynniki wirialne, jednakowy sposób 
wyprowadzania danych i.t.d. Program umożliwia również przetwarzanie re- 
zultatów obliczeń programem EXPVLE2 tzn. interpolacji całkowitej pręż- 
ności pary i składu fazy ciekłej dla zadanej temperatury. 


Próba zastosowania metody Тао lub МСС miała na celu uzyskanie obli- 
czenia równowagi termodynamicznej tzn. danych PTxy bez stosowania jakie- 
gokolwiek modelu termodynamicznego. Jak wiadomo, wybór modelu determinu- 
je wielkość i kształt nadmiarowej energii Gibbsa, wartości współczynni- 
ków aktywności i.t.d. Jak zostało powiedziane wcześniej, metoda Tao oraz 
МОС wymaga wygładzenia danych о prężności pary w funkcji stężenia. Zas- 


tosowano metodę funkcji sklejanych trzeciego stopnia (ang. cubic spli- 


Ы" ргасу tej nazwa ГЕМЕ jest ѕіоѕомапа м umownie м Когејасјасћ 2 ujem- 
nym dobieralnym parametrem а, w celu odróżnienia od korelacji równaniem 


NRTL z dodatnimi wartościami а. 
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nes) (Fortuna i wsp., 1982, Hofman, 1982). Јако drugi sposób wygładzania 
prężności pary nad roztworem zastosowano opis "nadmiarowego" ciśnienia w 


stosunku do ciśnienia wynikającego z prawa Raoulta: 


P + x, (1 = x Y y C (1 - ак)" VII-10 
1-0 
Sposób ten jest czesto stosowany przez grupe B.D.Smith'a (Khurma, 1983). 
Niestety, obie metody dawały rezultaty niezadowalające. W metodzie Tao 
obecność punktu azeotropowego powodowało trudność zbiegania sie funkcji 
minimalizowanej do zadanego kryterium. Następstwem tego było uzyskanie 
błędnej estymaty stosunku współczynników aktywności i błędnego policze- 
nia współczynnika w n*1 punkcie wektora stężeń. Dobitnie ilustruje to 


rysunek VII-1 na przykładzie obliczonego składu fazy parowej. 


1.00 


0.80 
~ 
Ф 
S 0.60 
ОГУ, —— Metoda Тао 
u 一 — Redlich-Kister ^ 
= NRTL (3 sale) „ 
a d 
©) 
0.40 
ZS 
= 
> 


0.20 
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Rys.VII-1. Porównanie obliczonego składu fazy parowej za pomocą metody 
Tao i równaniami Redlicha-Kistera i NRTL z trzema parametrami 
dobieralnymi dla izotermy krytycznej. 


Metoda MGC w gruncie rzeczy sprowadza sie do policzenia numerycznej 


pochodne j GE w punkcie n z wartości funkcji w punktach n-1 oraz n*1. 
Ponieważ układ jest bardzo niesymetryczny, trudno było założyć wyjściowe 
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przybliżenie nadmiaru energii Gibbsa. Układy o mniejszych odchyleniach 
od doskonałości mogły być korelowane przy startowej zerowej wartości сЕ 
W przypadku badanego uktadu, nawet zalozenie modelu roztworów regular- 
nych z odpowiednio dużą stałą ( >2 ) nie dawały zbieżnego szeregu roz- 
wiązań. Algorytmy powyższych metod przedstawione są w załączniku B. 


W związku z powyższym zdecydowano się na metodę Barkera i dwa 
równania korelacyjne: równanie Redlicha-Kistera i równanie NRTL. Model 


UNIQUAC jako zbyt mało elastyczny wyłączono z korelacji VLE. 


Jako procedury minimalizacyjne stosowane były: procedura ORT (Mala- 
nowski, 1974) rozwiązująca układ równań liniowych metodą ortogonalizacji 
wielomianów oraz procedura CURFIT opisana w książce Bevingtona (1968), w 
której zaimplementowany jest algorytm Marquardt'a (1963). Jako kryterium 
dopasowania równań korelacyjnych do danych doświadczalnych uwzględniane 


było odchylenie standardowe zdefiniowane: 


I [P = РА)? 
RMSD = "Tv Т VII-11 
gdzie n jest liczbą punktów doświadczalnych, k - liczbą dobieranych sta- 


łych. Drugim kryterium była wielkość x (chi-kwadrat ) 
ZE ров 2 
X = pci - : VII-12 
c 
1 1 


gdzie jest estymowanym błędem pomiaru i-tego punktu. W praktyce usta- 
la się średni błąd jakim obarczone są pomiary i ten stosuje się dla 
wszystkich punktów. W procedurze CURFIT, oprócz parametrów, obliczane są 
również ich wariancje c. Wartości ich jednak silnie zależą od uzyskanej 
wielkości x. Im wartość ta jest bliższa jedności, tym bardziej realne 
są wariancje współczynników. Potrojona wartość wariancji, gdy x 1, da- 
je пат 95% przedział ufności. W tablicach VII-1 i VII-2, obok wartości 
współczynników i ich wariancji, podane jest również RMSD i x. 


Równanie Redlicha-Kistera, z powodu swojej elastyczności wydawało 
się najodpowiedniejsze do korelacji tak trudnej zależności ciśnienia od 
składu. Wobec dużej liczby punktów eksperymentalnych istniała możliwość 
dobrania odpowiedniej liczby stałych równania. (Dozwolone było dobranie 


nawet 10 stałych, stosując kryterium co najmniej 3 punktów na jedną 
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stata.) W praktyce okazało sie, że nie było to konieczne. 


Na rysunku VII-2 przedstawione są rezultaty korelacji równaniem Re- 
dlicha-Kistera dla wybranych izoterm z tablicy VI-8; izotermy 298.133, 
303.153 i 309.15 K znajdują się w zakresie występowania dwóch faz 
ciekłych, izoterma krytyczna 311.146 К oraz izotermy 313.160, 318.150 1 
323.270 K, w których faza ciekła jest homogeniczna w całym zakresie ste- 
żeń. Dobór izoterm pozwalał pokryć zakres temperatur w miarę równomier- 


nie. 


W korelacjach zastosowano równanie NRTL w jego oryginalnej formie 
zaproponowanej przez Renona i Prausnitza (1968). Jak powiedziano, ist- 
nieje wiele modyfikacji tego rówania, szczególnie do polepszenia korela- 
cji równowag ciecz-ciecz, lecz w dalszym ciągu największą elastyczność 


posiada jego podstawowa forma. 


W tablicach VII-1 i VII-2 przedstawione są numeryczne rezultaty ko- 
relacji dla wspomnianych 7 izoterm za pomocą równań odpowiednio Redli- 
cha-Kistera i МЕТ. z dobieralnymi trzema parametrami. Korelacje równa- 


niem NRTL graficznie prezentuje rysunek VII-3. 


Z izoterm, w których występują dwie fazy ciekłe z korelacji 
usunięte zostały punkty, które tylko nominalnie reprezentują swoje 
składy. W rzeczywistości składy te są równe składom współistniejących 
faz. Punkty w tablicy VI-8 zaznaczone kursywą zostały właśnie usunięte z 


korelacji jako znajdujące się w obszarze dwufazowym. 


Dla trzech reprezentatywnych korelacji na rysunku VII-4 przedsta- 
wione są bezwzględne odchylenia obliczonych ciśnień w kPa. Linią przery- 
мапа zaznaczono estymowany z warunków przygotowania i prowadzenia pomia- 
rów błąd równy 20 Pa. Jak widać, rozrzuty przekraczające znacząco tę 
wielkość występują dla stężen większych od 0.5 ułamka molowego 3-MOPN. 
W zakresie tym składy roztworu tworzone są z małych naważek моду, nato- 
miast całkowite ciśnienie nad roztworem zdeterminowane jest właściwie 
prężnością cząstkową wody. Nawet niewielki błąd ważenia czy też przemra- 
2апіа w tym przypadku może spowodować przesunięcie w składzie rzeczywis- 


tym powodujący pozornie błąd w pomiarze ciśnienia. 
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Rys.VII-2. Wybrane izotermy równowagi ciecz-para. Korelacja równaniem 
Redlicha-Kistera. 
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Rys.VII-3. Wybrane izotermy równowagi ciecz-para. Korelacja równaniem 
NRTL. 
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Rys.VII-4. Bezwzgledne odchylenia cisnienia obliczonego od zmierzonego 
dla trzech izoterm. 
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Ponadto, dla wyższych temperatur występuje pewna systematyczność 
błędu i to bez względu na równanie korelacyjne, co pośrednio potwierdza 


powyższe wnioskowanie. 


Tablica VII-1. 


Wyniki korelacji równowagi ciecz-para układu (3-MOPN(1) + woda(2)) za 
pomocą równania Redlicha-Kistera dla wybranych izoterm. 


Tablica У11-1а. 


Temperatura 298.133 K 


R-K współcz. с 

A(1) 2.36161 0.00410 RMSD = 0.021 [kPa] 
А(2) 0.22853 0.00923 

A(3) 0.42858 0.01713 2 = 1.17 

А(4) 0.16984 0.02433 

A(5) 0.12251 0.02942 


Pexp AP G 
[kPa] [kPa] 7, [J/mol] 


3.167 20.000 27.415 
3.223 0.002 20.893 
3.240 -0.006 14.827 
3.243 -0.005 13.427 
3.248 0.001 10.367 
3.258 0.012 9.445 


luka mieszalności 


.237 -0.015 1.702 
.250 0.015 . 501 
.190 0.000 
.181 0.017 
.014 -0.040 
.814 0.037 
.445 -0.016 
.264 0.000 


"ooooooo 
mÍ"-"oOoooooo 
об злом IN I 
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Tablica VII-1b. 


Temperatura 303.153 K 


|Hooooooo 


R-K 

A(1) 
A(2) 
A(3) 
A(4) 
A(5) 
A(6) 


=ооооооо 


współcz. с 

2.29695 0.01062 
0.25623 0.02053 
0.06273 0. 07311 
-0.36378 . 0.10079 
0.37003 0.08377 
0.94093 0.16925 


obl. Pexp AP 
Y, [kPa] [kPa] 


0.000 
-0.014 
-0. 006 
-0.005 

0. 004 

0.016 
-0.000 


RMSD = 0.020 [kPa] 


х? = 0.989 


272 
504 


. 631 
.041 
. 090 
. 333 
6. 


181 


mieszalności 


0.007 
-0.015 
-0. 

0. 
-0. 

0. 

-0. 

0. 


| на на на MA а јез јез = 


.662 
. 450 
.273 
.229 
. 131 
. 047 
.016 
.000 


PZOPODE 


GE 


7, [J/mol] 


OFORBROODNO 
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Tablica VII-1c. 


Temperatura 309.150 K 


К-К współcz. 
A(1) 2.24638 
А(2) 0.31120 
A(3) 0.24255 
A(4) 0.23297 
A(5) 0.72062 
A(6) 0.24335 
A(7) -0.87400 


RMSD = 0.023 [kPa] 


x5 = 1.305 


Dual S um 


obl. Pexp AP d 


Yi [kPa] [kPa] Y2 [J/mol] 


22.716 

.017 19.275 
.003 14.655 
-0.001 13.377 
.003 10.327 
0.012 9.352 
-0.009 7:675 
-0.016 5.865 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


M кә кэ кэ кэ кэ нә ә 
N O O > омео 


mieszalności 


-0.014 .266 
. 007 .188 
. 009 . 022 
.013 .824 
.010 .619 
.016 .424 
.013 .252 
. 052 .209 
. 042 116 
. 021 . 034 
. 007 . 005 
. 000 . 000 


人 

олш m n с ооо 
M PR ра в» ра а = = NN DN 
АЕ НЕИМАР NWA 
O RANN O Q N н ооо 0 
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Tablica VII-1d. 


Temperatura 311.146 K 


RMSD = 0.021 [kPa] 


x5 = 1.138 


ob1. Pexp GE 


y, [kPa] [kPa] 7, [J/mol] 


Busonaaanananananasanananonanaana 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
da 


"OOOoooooooooooooooooooooooooooooo 


Ка Киа ра пе ке ода Ва Ко Ба баба Ко ВВ РТА WPI 
Оф Q Ñ) N N) = = къ на на на ка = = = вз вз = MR вз = PMP E = M = = M M = = 


| Очелеоолосоолооодмотом он =- мошо онолоо 
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Tablica VII-1e. 


Temperatura 313.160 K 


с 
.00714 RMSD = 0.024 [kPa] 
. 01074 А 
.01809 х = 1.414 
. 02565 
. 05559 
. 05315 
. 07890 


ођ1. Рехр AP 
Y, [kPa] [kPa] 


ngo dolci за аа ci cd d ДО A AWA AAA AN AA R A asa 
= = = 
оргререргреооооооовососоввеллвовмој) 


О\ ЛЬ Q NN NN) NN) кә к» к» къ S мз 挛 
O оосо S O OA K AUAN ос лю ш ON NO = Q N = @ G@ = оллоо 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


=осоосооооссоооооосоосооооооосооососооо 
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Tablica VII-1f. 


Temperatura 318.150 K 


R-K współcz. с 

А(1) ‚ 15782 .00713 RMSD = 0.026 [kPa] 
А(2) . 36916 . 00952 

А(3) 21760 01478 x^ = 1.705 

А(4) . 34191 02728 

А(5) . 42937 04296 

А(6) . 39935 . 07045 

А(7) 69756 . 06296 

А(8) . 96095 . 08322 


1 
ооооооом~ 
оооооооо 


obl. Pexp AP G° 


у, [kPa] [kPa] 7, [J/mol] 


о 10000 0 0 0 O O O O O O O O O O O 0 O O O O O O 0 O O ою 
= = = = = = = Ñ Ñ Ñ Ñ Ñ Ñ N о @ U W ооо ао 0 0! 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


н ооооооооооооооооооооооооооооорооо 
O O O O + == = Q Q N осо с со SŁ O N Q O Q O X N N O N O о lO QI Q i O 


NARUN мо МУ кә dO кэ кэ ка ка на на IH на рз ра рэ зэ EI кэ ва вэ кэ I| кэ I| | I| P M 


= 
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Tablica VII-1g. 


с 
.00925 RMSD = 0.034 [kPa] 
01084 А 
01776 x = 2.840 
03833 
.04564 
10087 


.07827 
.12700 


АР GE 


[kPa] 7, [J/mol] 


о @ OQ O O @ O O Ə O Ə O O O O O O O O 
=O юе он оомонњномоооо 
Фо N = Q 60 Q N O = ЈУ N í ртн м 


M 
тик ко Ри ва БӨГӘ КУ S ро ОУ 00 ораи Бал Вид q tO Ко L 


ооонооо 
O N сюе (Q b O 
о - с 00 ил 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


гоооооооооооооооооооооооооооооооо 
О -ч О о е O Ñ & Q O X Q O => Q O O N AEN ONA X O Q > = оњ N O 


= (Л > (Q) N) Ñ м = = къ къ e M M M M M M ыз M M а M e e = = == = = = 


— 
FEE". 
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Tablica VII-2. 


Wyniki korelacji rownowagi ciecz-para uktadu (3-MOPN(1) + woda(2)) za 
pomoca rownania NRTL z dobiralnymi trzema statymi dla wybranych izoterm. 


Tablica VII-2a. 


МЕТ. bezw.wsp. с 
C(21) 2.46912 0.02736 RMSD = 0.020 [kPa] 


C(12) 1.66494 0. 01093 > 


alfa 0. 4295 0. 0012 x = 1.035 


obl. Pexp AP G 


Y, [kPa] [kPa] 75 [J/mo1] 


3.167 0.000 26. 
3.223 0.005 20. 
3.240 -0.001 14. 
3.243 -0.000 13. 
3.248 0.006 10. 
3.258 0.017 9. 


luka mieszalnosci 


. 237 -0.026 
‚250 ‚007 
. 190 . 002 
. 181 ‚021 
. 014 -0.033 
.814 0.038 
. 445 . 022 
.264 0.000 


cooooooo 
"ooooooo 
M кэ кэ кэ а HEB 


O Q N Q O Q 0 0 
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Tablica VII-2b. 


Temperatura 303. 153 K 


NRTL bezw. wsp. с 

C(21) 2.45989 0.01988 RMSD = 0.025 [kPa] 
C(12) 1.62811 0.00819 е 

alfa 0. 4336 0. 0009 x = 1.500 


GE 
7, [J/mol] 


5 . 000 

. 322 .006 19. 
. 349 .001 14. 
.350 -0.000 12. 
.355 0.006 10. 
.365 0.018 9. 
.341 -0.001 7. 


luka mieszalnosci 


-0. 028 1. 
0.007 
-0. 000 
0.029 
-0.044 
0.049 
-0.027 
0.000 


-»=осоооооо 
"ooooooo 
O O O => N = BO 
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Tablica VII-2c. 


Temperatura 309.150 K 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


=оосоогооооооо 


МЕТ. bezw.wsp. 
C(21) 2.37800 
C(12) 1.53179 
alfa 0. 4305 


obl. Pexp 
Yi [kPa] 


. 948 
. 052 
. 091 
. 093 
. 098 
. 105 
. 083 
.084 


ови ао 


luka mieszalności 


. 084 
096 
081 
072 
049 
024 
878 
858 
497 
. 073 
. 370 


> m m m n 9 DD D D 


с 


0.00894 
0.00548 
0.0007 


AP 


[kPa ] 
.000 23. 


=0 


0. 
-0. 
-0. 
-0. 
-0. 
-0. 

0. 
-0. 

0. 
-0. 


0 


. 020 
.019 
.019 
. 027 
. 037 
. 023 
. 033 


. 006 
004 
016 
024 
032 
011 
031 
014 
091 
064 
021 
000 


: .520 ; 2 


RMSD = 0.035 [kPa] 


x 


18. 
13. 
12. 
9; 
8. 
7. 
6. 


= к» къ к» = кв» = = = М М М 


= 3.051 


ке на кэ кэ — вз 


GE 


7, [J/mol] 


олалом~омо 


O л 00 о О (л (л О > ол W 
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Tablica VII-2d. 


Temperature: 311.146 K 


NRTL bezw. мѕр. с 
C(21) 2.42638 0.00582 RMSD = 0.027 [kPa] 
C(12) 1.57257 0. 00523 " 
alfa 0. 4406 . 0006 X = 1.869 
obl. Pexp AP G 
Y, [kPa] [kPa] Ya [J/mol] 


оооооооооооооооо 


Lue СЕК 
ые и m o T > D а 


1 
o 


Bunnonananananananananananananaa 
° 


1 
> 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


ee 


пе посна На Бр бари бо mI бо ној BUS бо Ио 
O — Ñ O — O [Ñ O O. (n @ (л — O О O (n O юе = ою к= охо O (n i> юе O () x O 


m 
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Tablica VII-2e. 


Temperature: 313.160 K 


NRTL  Беги. wsp. с 

С(21) 2.42243 0.00524  RMSD = 0.030 [kPa] 
C(12) 1.55633 0.00469 , 

alfa 0.4422 0.0005 x? = 2.246 

obl.  Pexp АР g" 
у, [kPa] [kPa] 7, [J/mol] 


оооооооооооооооо 


IKAR Тик CIS К ШЫ | 
ооооооооооооо 


| 
o 


人 
e 


1 
x 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


гоооооооооооооооооооооооооооооооо 
=== 
= = = = = = = = = Ñ Ñ Ñ Ñ Ñ Ñ N Q Q @ W @ Q O > > Qm — @ o Ñ Q @ i 


O00 > о N AO O N O QA OA AO OQ O Q O N P.O = O G N N S P AO O G G O 


e 
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Tablica VII-2f. 


110 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


NRTL 
C(21) 
C(12) 
alfa 


| нон ааоововиооовоово оо 2 


bezw. wsp. 
2.. 41685 
1. 51684 
0. 4463 


с 
0.00408 
0.00361 
0.0004 


оооооооооооооооо 


кыз ОКА Де ЖА or 1. 
ооооооооообоо 


RMSD = 0.038 [КРа] 


x 


EP MP 
PRRRRNRRRRNNNNYNNNÓGCÓWOWORRUJWONÓG 


2 


= 3.722 
E 
7, [J/mol] 


| O N O + @ O. O JVOO + O i> QI O i> => O O + X @% ого = @ 0 00-50 


Tablica VII-2g. 


Temperature: 323.270 K 


МЕТ. bezw.wsp. с 

C(21) 2.49050 0.00334 RMSD = 0.057 [kPa] 
C(12) 1.49419 0.00282 А 

alfa 0. 4508 0. 0003 x = 8.263 

obl. Pexp AP ur 
Y, [kPa] [kPa] Ya [J/mol] 


оооооооооооооо 


"OOOoooooooooooooooooooooooooooooo 
ОЈ O л N о (Q @ G (n O N O @ G@ Q Ñ GQ GQ N G Z .O O. O O. QO O > X > O 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
1. 


Tablica VII-3. 


Wartości stałych równania NRTL na podstawie oddzielnych korelacji danych 
VLE poszczególnych izoterm. 


2. 63568 
2.46912 
2. 45989 
2.53021 
2.51673 
. 37800 
. 50448 
. 50147 
. 49971 
. 49884 
. 42638 
. 49817 
. 49580 
. 49329 
. 42243 
. 41685 
. 49050 


кә кэ кэ кэ HH рэ рэ рз рз рэ үз үз На | рэ ра 


2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 


Izoteme krytyczna zaznaczono gwiazdka. 


3.00 


ro 
сл 
° 


N 
© 
© 


(л 
© 


x C12 
+-+ +-+ + glfa 


° 
© 


О 
сл 
© 


Bezwymiarowe parametry МЕТ. 


.00 
290.00 300.00 310.00 320.00 330.00 
Temperatura (K) 


Rys.VII-5. State równania NRTL w funkcji temperatury na podstawie kore- 
lacji VLE. 
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Jzyskane z korelacji parametry dla wszystkich (16) izoterm układają sie 
liniowo w funkcji temperatury. Przedstawia to rysunek VII-5., a nume- 


ryczne warości (bez wariancji) zestawione są w tablicy VII-3. 


VII-4 KORELACJE RÓWNOWAG CIECZ-CIECZ 


Do korelacji równowag ciecz-ciecz użyto tylko własnych danych 
(Luszczyk, 1988), wygładzonych równaniem V-3. Korelacja polegała па 
znalezieniu stałych równania NRTL, ГЕМЕ i UNIQUAC ze składów równowago- 
wych dla poszczególnych (9) temperatur. Stałe równań korelacyjnych 2пај- 
cowano z rozwiązania układu równań VII-2, w których obie funkcje zdefi- 


riowane są następująco: 


8 7, 
F = In PT ln =F i=1,2 VII-13 
i x 7 
1 1 
gdzie ' i '' oznaczają wielkości w poszczególnych fazach. Taka forma za- 


pisu ma tę zaletę, Ze w procesie iteracyjnym obliczany jest tylko 
składnik zawierający współczynniki aktywności. W przypadku równania NRTL 
konieczne było założenie trzeciego parametru. Przyjęto « = 0.441, 
wartość zbliżoną do uzyskanej dla izotermy krytycznej w korelacjach VLE 
(zaznaczona gwiazdką w tablicy VII-3), a jednocześnie jeszcze dającą 
rozwiązanie warunków VII-14. Według pracy Mettelina i Verhoeye (1975) 
istnieje pewna graniczna wartość, powyżej której nie można uzyskać roz- 
wiązania warunku krytycznego. Wyniki korelacji przedstawia rysunek УТТ-6 
natomiast w tablicy VII-4 przedstawione są wyniki numeryczne w postaci 
bezwymiarowych stałych dla poszczególnych temperatur oraz sumy absolut- 
nych wartości funkcji VII-13. Podane są również składy współist- 
niejących faz. 


Tablica VII-4 zawiera również stałe (oznaczone gwiazdką) obliczone 
na podstawie składu mieszaniny w punkcie krytycznym oraz warunku stanu 


krytycznego w układzie ciekłym. 


=0, =0 VII-14 
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W załączniku С podane sa explicite wzory na powyższe pochodne wyprowa- 


dzone z obu modeli. 


Tablica VII-4. 


Wartości bezwymiarowych parametrów NRTL (« = 0.441) wyznaczonych z Коге- 
lacji równowag ciecz-ciecz. 


RS 
NNNNNNNNN 
popmpepuebp 


оооооооо 


2.50 


bon. 


юе Ci 


== Z 
alfa=0.441 » 
x 


Bezwymiarowe parametry NRTL 


2.00 
300.00 30200 304.00 306.00 308.00 310.00 312.00 
Temperatura (K) 


Rys.VII-6. Zależność temperaturowa stałych МЕТ. na podstawie izoterm 
równowagi ciecz-ciecz obliczonych dla stałej wartości рага- 
metru < = 0.441. 
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Podobne obliczenia przeprowadzono równiez dla orginalnego (« = -1) 
modelu LEMF i UNIQUAC, a rezultaty przedstawione sa, odpowiednio, w tab- 
licach VII-5 i VII-6 oraz na rysunku VII-7 i VII-8. Rowniez tutaj state 
wyznaczone z warunku VII-14 oznaczone są gwiazdką. 


Tablica VII-5. 
Parametry równania LEMF (« = -1) uzyskane z korelacji LLE. 


. 491942 * 0.952995 
. 399375 . 946886 
. 398737 . 949333 
. 403290 . 950866 
.415829 . 952176 
. 442763 . 953104 
. 466450 .953374 
. 507893 . 953865 
. 550764 . 955862 

. 964927 


ОООООООООО 


1.40 


Bezwymiarowe parametry ГЕМЕ 


0.20 
298.00 302.00 3506.00 310.00 


Temperatura (K) 


Rys.VII-7. Zależność temperaturowa stałych ГЕМЕ uzyskane z korelacji 
danych LLE. 
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Tablica VII-6. 
Parametry równania UNIQUAC uzyskane z korelacji danych LLE. 


. 270897E-01 
. 901275E-01 
. 673363E-01 
. 647758Е-01 


. 763430 
. 920260 
. 866509 
. 868672 
. 898278 
. 909081 
.954949 
.977011 
. 03033 

. 08829 


. 723780Е-01 
. 690049E-01 
. 800659E-01 
. 786705Е-01 
. 850637E-01 
. 892212E-01 


M-O—OOOOOOOE 
оооооооооо 


e 
сл 
о 


| 
© 
o 
о 


Bezwymiorowe parametry ЏМОЏАС"а 


—0.50 
295.00 300.00 305.00 310.00 315.00 


Тетрегаїига (К) 


Вуз. УТТ-8. Zależność temperaturowa parametrów równania UNIQUAC na pods- 
tawie korelacji równowagi ciecz-ciecz. 


Ponieważ dysponowano również  prezno$ciami рагу nad roztworami 
dwufazowymi możliwe było wyznaczenie na podstawie tych danych parametru 


« równania NRTL. Jako trzecią funkcję układu równań VII-9 zastosowano: 


Г =P - P(A ) VII-15 
3 eksp NRTL 


gdzie A jest zbiorem parametrów NRTL. 
NRTL 
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Tablica VII-7. 


Wartości bezwymiarowych paramerów МЕТ! wyznaczonych z korelacji równowag 
ciecz-ciecz i prężności pary nad roztworem dwufazowym (korelacja LLE+VP1). 


T/K  P/[kPa] х, 


оооооооооо 
юююююююююю 


овилаттто 
оооооооооо 
woRRRRORmwww 
оооооооооо 


4.00 


2.00 
*-+-+-+-+ alfa 
C21 


+ + — C12 


Bezwymiarowe parametry NRTL 


0.00 
298.00 302.00 306.00 310.00 
Temperatura (K) 


Rys. VII-9. Zależność temperaturowa stałych równania NTRL uzyskana ze 
wspólnej korelacji danych ЏЕ i prężności pary nad luką mie- 
szalności (VP1). 


Rysunek VII-9 i VII-10 przedstawia zależność temperaturowa uzyska- 
nych parametrów, natomiast tablica VII-7 i VII-8 ich wartości numeryczne 
oraz, podobnie jak w tablicy VII-4, sumę absolutnych wartości minimali- 
zowanych funkcji. Podane są również użyte do korelacji składy faz oraz 
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odpowiadajace im ci$nienia. Sa to wartosci srednie 2 pomiarow w obszarze 


dwufazowym. 


Tablica VII-8. 


Wartości bezwymiarowych parametrów NRTL uzyskanych z korelacji równowagi 
ciecz-ciecz i prężności pary nad roztworem dwufazowym (korelacja LLE+VP2). 


T/K  P[kPa] х, 


оооооооооо 
COPNNNNNND 


oooooooooo 
NNNNNNNNNN 
оооооооооо 


6. 
6. 
6. 
6. 
6. 
5. 
5. 
4. 
3. 


3.00 


*-* +-+ + alfa 


Bezwymiarowe parametry МКТ 


` 297.00 301.00 305.00 309.00 
Temperatura (K) 


Вуз. У11-10. Rezultaty korelacji LLE+VP2 dla drugiej gałęzi krzywej współ- 
istnienia za pomocą równania NRTL. 


Rezultaty z tablicy VII-7 uzyskane zostały z wykorzystaniem składów 
równowagowych lewej (VP1) gałęzi krzywej współistnienia (3 kolumna). Po- 
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dobne korelacje zostały przeprowadzone dla prawej (VP2) gałęzi krzywej 
używając tych samych ciśnień całkowitych. Rezultaty zawarte są w tablicy 
VII-8 i na rysunku VII-10. 


Tablica VII-9. 
Parametry równania ГЕМЕ uzyskane korelacji ШЕ i VP1. 


. 849105 
. 832230 
. 822390 
. 814991 


. 826579 
. 827947 
. 829727 
. 831855 
. 835494 
. 842079 
. 847049 
. 856249 
. 867950 
. 886420 


. 806399 
. 795124 
. 791503 
. 784580 
. 774362 
. 752558 


| ж» @ O: Q: = = Ж ооо 
оооооооооо 
оооооооооо 


= ГЕМЕ alfa 


+ += C21 


v. C12 


一 0.40 


—0.90 


Bezwymiarowe parametry LEMF 
© 
o 


- 1.40 
298.00 302.00 306.00 310.00 


Temperatura (K) 


Rys.VII-11. Zależność temperaturowa stałych LEMF uzyskana z korelacji 
LLE+VP1. 
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Podobnie jak dla modelu NRTL, wspólne korelacje LLE+VP1 i LLE+VP2 
przeprowadzono dla modelu ТЕМЕ. Należy przypomnieć o umownym znaczeniu 


nazwy ТЕМЕ, która oznacza dobieranie ujemnej wartości parametru «. 


Tablica VII-10. 
Parametry rówania LEMF uzyskane z korelacji LLE+VP2 


. 779073 
832230 
. 822390 
. 814991 


. 870382 
. 827947 
. 829727 
. 831855 
. 835494 
. 842079 
. 847049 
856249 
. 867950 
. 886420 


806399 
795124 
. 791503 
784580 
. 774362. 
. 752558 
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Rys.VII-12. Zależność temperaturowa stałych ГЕМЕ uzyskana z korelacji 
LLE+VP2. 


120 


VIII DYSKUSJA I WNIOSKI 


Przed omówieniem wyników korelacji, w celu uzyskania petniejszego 
obrazu badanego układu, omówione zostaną właściwości fizykochemiczne 
czystego  3-metoksypropionitrylu i jego zachowanie się w układach 
dwuskładnikowych. W dalszej części rozdziału omówione są, poza wnioskami 
wynikającymi z korelacji, inne aspekty wyników tej pracy. Wobec bogatego 
materiału doświadczalnego dostępnego w tej pracy możliwe jest przedysku- 
towanie zagadnienia zależności ciśnienia od składu w otoczeniu punktu 
krytycznego, czy też wpływu poprawki do nadmiarowej entalpii Gibbsa 


związanej z występowaniem objętości nadmiarowe ј. 


VII.1. 3-METOKSYPROPIONITRYL 


3-metoksypropionitryl (3-MOPN, C4HZNO. CAS index: [110-67-8]) jest 
rozpuszczalnikiem dwufunkcy пут, zawierającym grupę eterowa i nitrylowa, 
należący do tzw. rozpuszczalników aprotycznych takich jak np. N-metylo- 
pirolidon (NMP) czy dwumetyloformamid (DMF). Jego struktura chemiczna 


jest następująca: 


Lowe i Kewley (1976) na podstawie widm MMR stwierdzili, że 3-MOPN jako 
ciecz w przeważającej ilości występuje w postaci izomeru trans-. Posiada 
duży moment dipolowy (3.47 D) wynikający z obecności grupy nitrylowej. 
Temperatura wrzenia (ок. 438 K) wskazuje na znaczne oddziaływania mie- 
dzycząsteczkowe. Obecność grupy nitrylowej prawdopodobnie powoduje two- 
rzenie w fazie ciekłej dimerów powstałych przez antyrównoległe ułożenie 
się grup, analogicznie do stwierdzonych w innych nitrylach. Związek ten 
dobrze rozpuszcza się w alkoholach, acetonie, niskocząsteczkowych nitry- 


lach, jak również w węglowodorach aromatycznych (Antosik, 1991). W nas- 
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yconych weglowodorach tañcuchowych i cyklicznych (Antosik i wsp., 1990b) 
jak równiez w weglowodorach posiadajacych elektrony т wykazuje ograni- 
стопа rozpuszczalność z tym, że obecność wiązań podwójnych znacznie ob- 
niża temperaturę krytyczną mieszalności (Łuszczyk i Stryjek, 1984a). Z 
tego względu jest potencjalnym rozpuszczal-nikiem w destylacji ekstrak- 
cyjnej. (Davis i Makin, 1969; Łuszczyk i Stryjek, 1978; Stryjek i 
Łuszczyk, 1978a). 


VII-2 KORELACJE VLE I LLE 


Głównym celem pomiaru równowagi  ciecz-para, a następnie jej 
korelacji, szczególnie, gdy uzyskujemy w czasie pomiarów niepełne dane 
równowagowe PTx, jest: 

- obliczenie pełnej równowagi dla eksperymentalnego ciśnienia czy 
temperatury; 

- na podstawie uzyskanych wyników przewidzenie zachowania się 


układu w innych warunkach PTx. 


Pierwsza część zadania jest stosunkowo prosta; polega na doborze 
odpowiedniego równania korelacyjnego, takiego aby pozwoliło odtworzyć 
zmierzone wielkości z zadaną dokładnością (błędem). Wtedy, gdy zmierzone 
wielkości nie są odtwarzane należycie, konieczne jest poszukanie innego, 


bardziej adekwatnego równania. 


Rozwiązanie drugiego zagadnienia jest o wiele trudniejsze. Wymaga 
zastosowania równania opartego na pewnych, choćby półempirycznych, род5- 
tawach teorii stanu ciekłego. Jego parametry powinny posiadać określoną 
przynajmniej zależność od temperatury. W literaturze występuje wiele 
równań spełniających bardziej lub mniej taki warunek lecz chyba najpopu- 
larniejsze są, omawiane w rozdziale VII, rówania МЕТ. i UNIQUAC (Lafya- 
tis i wsp., 1989). 


Najczęściej problem korelacji ciecz-ciecz i ciecz-para rozpatrywane 
były rozdzielnie. Próby wspólnego ich traktowania pojawiły się wraz z 
równaniem NRTL, a praktycznym przykładem jest program NRTLB z książki 
Renona i wsp. (1971). Niestety, brak było danych LLE i VLE pokrywających 


ten sam zakres temperatur w całym zakresie stężeń. Z konieczności więc 
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można było badać własności równania NRTL i jego zdolność korelacyjna na 
ograniczonych zbiorach danych, takich jak częściowe dane VLE, graniczne 


współczynniki aktywności czy entalpie mieszania. 


W latach 80-tych Łuszczyk i Stryjek (1981) próbowali uzyskać kon- 
systentny opis równowagi VLE i LLE dla układów węglowodory (cz (1) + 3- 
MOPN(2)). Na podstawie badań literaturowych oraz własnych można powie- 
dzieć, że równania wyrażające nadmiarową energię Gibbsa oparte na kon- 
cepcji stężeń lokalnych mogą być z powodzeniem stosowane do jednoczesne- 
go opisu równowag ciecz-ciecz oraz ciecz-para tylko dla warunków даје- 
kich od stanu krytycznego. W pobliżu tego punktu korelacja prowadzona 
oddzielnie dla każdego typu równowagi, czy też jednocześnie, prowadzi do 
wyników niezadowalających. Obliczone prężności pary wykazują przegięcia 
i minima, natomiast obliczona temperatura krytyczna czy składy współist- 


niejących faz są niezgodne z obserwowanymi doświadczalnie. 


Zbadanie diagramów fazowych uzyskanych z równań korelacyjnych prze- 
prowadzone zostało z konieczności dla niewielkiej liczby układów dwu- 
składnikowych. Dane VLE i LLE uzyskane we wspólnym przedziale temperatur 
są niezwykle rzadkie, a takich, które uwzględniają jednocześnie punkt 
krytyczny mieszalności, praktycznie brak. Dysponowano danymi obejmujący- 
mi równowagi ciecz-ciecz (Łuszczyk i Stryjek, 1984a) i odpowiadające im 
ciecz-para (Stryjek i wsp., 1978c) dla układów dwuskładnikowych serii 


węglowodorów С, + 3-MOPN. Zbadane zostały luki mieszalności pentanu, 


izopentanu, EMO) cis- i trans- 2-pentenu, 2-metylo- i 3-metylo- 
butenu-1 z 3-MOPN. Dysponowano danymi o równowadze ciecz-para dla tych 
układów w trzech temperaturach: 293, 313 i 333 K. Mimo tego, że dane VLE 
nie są kompletne, tzn. nie pokrywają całego zakresu stężeń i nie są w 
bezpośrednim sąsiedztwie punktu krytycznego mieszalności, wykazały pewne 
prawidłowości w korelacjach przy użyciu wyrażeń na GE w ich oryginalnych 


formach. 


Skorelowano luki mieszalnošci równaniem NRTL z liniowa temperaturo- 
wą zależnością stałych. Z uzyskanych parametrów obliczono prężności par, 
które porównano 2 prężnościami doświadczalnymi. Okazało sie, że im izo- 
terma bliższa temperaturze krytycznej, tym gorsze odwzorowanie ciśnienia 


eksperymentalnego. Szczególnie jest to dobrze widoczne na przykładzie 
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izopentanu z 3-MOPN (Тт ok. 370 K). Izoterma 293 К wykazywala odchylenia 
rzędu 3 4, a dla 333 K już 15 X. Dla układów z olefinami, których tempe- 
ratury krytyczne leza w rejonie 293-318 K odchylenia ciśnienia 
całkowitego dochodzą do 60 % dla izoterm w 313 K. Jeśli odchylenia 
odłoży się w funkcji temperatury zredukowanej względem UCST badanych 
układów, to okazuje się, że leżą na wspólnej krzywej. Świadczy o jedna- 
kowym zachowaniu się obliczonych odchyleń, bez względu na układ. Dla 
wszystkich korelacji używano liniowej zależności parametrów od tempera- 
tury. Gdy z danych o rozpuszczalności wyinterpolowano numerycznie pro- 
gramem, opartym na równ. V-3, składy równowagowe, na podstawie których 
również można obliczyć parametry równania NRTL, to odłożone w funkcji 
temperatury dawały zależności raczej wykładnicze niż liniowe. W litera- 
turze spotyka się próby zastosowania innych zależności temperaturowych 

aby poprawić opis krzywej współistnienia, lecz wpływ takiej zmiany na 
obliczone prężności par lub jednoczesna korelacja LLE i VLE nie były ba- 


dane. 


8.00 


) 
° 
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Rys.VIII-1. Porównanie rezultatów korelacji LLE za pomocą różnych modeli 
z danymi eksperymentalnymi dla izotermy krytycznej. 


Celem niniejszej pracy było uzyskanie w miarę konsystentnego opisu 


układu wykazującego niemieszalność. Przeprowadzone obliczenia wykazały, 
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że zarówno równanie МЕТ. ze stałą о < 0.441, јак i UNIQUAC nie są w sta- 
nie odtworzyć ciśnienia nad roztworem dwufazowym i w pobliżu UCST z błę- 
dem eksperymentalnym. Ilustruje to rysunek VIII-1, na którym przedsta- 
wiono porównanie ciśnienia eksperymentalnego z obliczonym wymienionymi 
równaniami. Parametry energetyczne uzyskane zostały z warunku stanu kry- 
tycznego (wzory VII-14). O ile w przypadku równania NRTL istnieje 
możliwość doboru trzeciego parametru tak, aby obliczone ciśnienie było 
zgodne z doświadczalnym, to w przypadku równania UNIQUAC, jako modelu 
dwuparametrowego, takiej możliwości nie ma. Można próbować zmieniać 
liczbę koordynacyjną z, lecz jej wpływ na wielkość ciśnienia jest mini- 
malny. Z tego względu model ten nie był szerzej stosowany w tej pracy. 
Podobnie zachowuje się oryginalne równanie ГЕМЕ (а = -1), nie przedsta- 


wione na rysunku. 

W niniejszej pracy przeprowadzono cztery korelacje równaniem NRTL 
używając różne zestawy danych w tym samym zakresie temperatur. Uzyskane 
stałe tego równania zestawione sa na wspólnych wykresach dla poszczegól- 


nych parametrów na rys. VIII-2 do VIII-4. 
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Rys.VIII-2. Zestawienie uzyskanych stalych с równania NRTL 2 róznych 
wariantów korelacji. 
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Rys.VIII-3. Zestawienie uzyskanych statych C12 równania NRTL z róznych 
wariantów korelacji. 
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Rys.VIII-4. Zestawienie uzyskanych stałych C21 równania МЕТ. z różnych 
wariantów korelacji. 
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Na przedstawionych rysunkach widaé wzajemna konsystencje uzyskanych 
stałych. Przebieg parametru С uzyskanego 2 korelacji МЕ jest co praw- 
da rózny niz uzyskanego z korelacji LLE, ale wynika to z zakrzywienia i 
asymetrii krzywej binodalnej w poblizu punktu krytycznego. Nalezy pod- 
kreślić, że podczas korelacji równowag ciecz-para do równania korelacy j- 
nego nie jest "przenoszona" informacja o tym czy roztwór jest jedno- czy 
dwu fazowy. Jak powiedziano, ciśnienie ma charakter pola w sensie termi- 
nologii Griffiths'a i nie wykazuje żadnych anomalii w otoczeniu punktu 
krytycznego mieszalności. Przeciwnie jest w przypadku korelacji LLE, 
gdzie ta informacja przekazywana jest w składach wspólistniejących faz. 


Z punktu widzenia korelacji równowag ciecz-ciecz entalpia swobodna 
mieszania g" jest funkcja dostarczajaca najwiecej informacji. Obraz jej 
przebiegu w funkcji stężenia pozwala ocenić zastosowany do korelacji mo- 
del. Chociaż przy określaniu parametrów modelu na podstawie danych PTx 
ciśnienie całkowite stosujemy właściwie jako kryterium ich doboru (po- 
przez odpowiednie wyrażenia termodynamiczne), to o przydatności konkret- 
nego modelu można sądzić dopiero po dogłębnej analizie innych wielkości 
termodynamicznych, które mogą być na jego podstawie obliczone (np. e. 
Termodynamika precyzyjnie okre$la jakie sytuacje sa dozwolone, a јакје 
nie i jakiekolwiek odstępstwa od tych prawideł mogą zdyskwalifikować mo- 


del, lub co najmniej ograniczyć jego stosowalność. 


Na poniższych rysunkach przedstawione są obliczone G* na podstawie 
stalych równania Redlicha-Kistera (R-K) i NRTL podanych w tablicach VII- 
1 do VII-6. Rysunki VIII-5 i VIII-6 pokazuja przebieg с“ na podstawie 
korelacji МЕ. Zależność temperaturowa i stezeniowa ma charakter raczej 
monotoniczny dla obu równań, poza izotermą 303.153 K dla równania R-K. 
Związane jest to prawdopodobnie z użycie nieodpowiedniej liczby stałych. 
Potwierdza to zależność obliczonego ciśnienia (rys.VII-2) dla tej izo- 


termy. 


Na wszystkich rysunkach obrazujących zależności prezentowanych 
wielkości w funkcji składu pionową, przerywaną kreską zaznaczony jest 
skład krytyczny „= 0.2392 (podobnie jak w rozdz. VII). Ułatwia to 


umiejscowienie tego punktu na analizowanych krzywych. 
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Rys.VIII-5. Obliczona funkcja GM na podstawie parametrow rownania 
Redlicha-Kistera. 
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Rys.VIII-6. Obliczona funkcja GM na podstawie parametrów NRTL uzyskanych 
z korelacji danych VLE. 
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Rysunek VIII-7 przedstawia te sama funkcje obliczona ze statych 
МЕТІ. z korelacji LLE dla ustalonej wartości parametru «œ = 0.441. Widać 
tutaj jakościową zmianę stężeniowej zależności. Brak informacji o 
własnościach homogenicznych części roztworu powoduje, że przebieg 
funkcji wynika tylko z formy algebraicznej równania korelacyjnego i nie 
ma nic wspólnego z jej rzeczywistym przebiegiem. Np. na podstawie tego 
rysunku można wysnuć wniosek, że obok górnego punktu krytycznego 
istnieje jeszcze dolny punkt krytyczny mieszalności w zakresie wyższych 
stężeń 3-MOPN. Oczywiście to zjawisko nie jest obserwowane. 


I 
| 
| 一 — izoterma krytyczna 


50.00 


0.00 


—50.00 


—100.00 


—150.00 


GM (J/mol) 


311.143. 


—200.00 


—250.00 


NRTL, ЦЕ, alfa=0.441 
.00 0.20 0.40 0.60 0.80 1.00 


x(1) 


— 500.00 
0 


Rys.VIII-7. Obliczona funkcja GM na podstawie parametrów NRTL uzyskanych 
z korelacji danych LLE. 


Szczególnie interesujace jest porównanie obliczonej funkcji G" dla 
korelacji LLE+VP1 i LLE+VP2 (rozdzial VII-3). Z punktu widzenia termody- 
namiki nie powinno mieé znaczenia, czy cisnienia nad roztworem dwufazo- 
wym brane sa wraz z lewa czy prawa galezia krzywej wspótistnienia (Van 
Ness i Abbot, 1982). Na podstawie rysunków VIII-8 i VIII-9 widać, że tak 
nie jest, i to, stezenia której gatezi zastosujemy, w przypadku równania 


NRTL, ma znaczenie. 
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Korelacja (LLE+VP1) dla stezeñ wzietych z lewej cze$ci Кггумеј 
współistnienia (niskie stężenia 3-MOPN) powoduje powstanie niespodziewa- 
nych ekstremów, zupełnie nieuzasadnionych. Chociaź spełnione są kryteria 
zerowania się funkcji, co można ocenić na rysunku VIII-8 i widać w ta- 
blicy VII-7 (ostatnia kolumna) rezultaty nie są zadowalające. Można dla 
każdej izotermy przeprowadzić kilka linii stycznych wyznaczających 
składy równowagowe niezgodne z obserwowanymi doświadczalnie. Zupełnie 
inny charakter posiadają przebiegi gh jeśli do korelacji weźmiemy 
stężenia dla prawej części krzywej współistnienia (korelacja LLE+VP2). 
Zależność ta ma charakter monotonicznej zmiany wraz z temperaturą, dla 
niższych temperatur podobną do uzyskanej z pomiarów VLE. Około 3 stopni 
poniżej izotermy krytycznej i w tym przypadku zaczynają pojawiać się 
przegięcia sugerujące powstawanie dolnego punktu krytycznego mieszal- 


ności w zakresie wyższych stężeń 3-MOPN. 


Inną więlkością pozwalającą ocenić rezultaty jest oczywiście 
ciśnienie całkowite nad roztworem. Na rysunkach VIII-10 i VIII-11 przed- 
stawiono ciśnienia obliczone ze stałych uzyskanych z powyższych Когеја- 
cji” . Widoczna jest znaczna rozbieżność dla ciśnień obliczonych i eks- 
perymentalnych we wszystkich temperaturach w przypadku korelacji VLE+ 
УР1, czego należało oczekiwać na podstawie przebiegów б". Inaczej sytua- 
cja wyglada w przypadku korelacji LLE+VP2. Tutaj ciénienia zmieniaja sie 
monotonicznie ze wzrostem temperatury, chociaz w miare dochodzenia do 
temperatury krytycznej ich przebieg staje coraz bardziej niezgodny z ob- 
serwowanym. Również dla tego zakresu temperatur przebieg parametrów a i 
тс równania korelacyjnego jest silnie nieliniowy, mocno odbiegajacy od 


uzyskanego dla parametrów uzyskanych z korelacji VLE. 


* Na rysunkach przedstawiono równieź ciśnienia eksperymentalne dla 


wyższych izoterm w celu zachowania skali z następnymi rysunkami. 
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Rys VIII-8. Temperaturowa zależność GM па podstawie stałych równania 
NRTL uzyskanych z korelacji LLE+VP1. 
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Rys.VIII-9. Temperaturowa zależność GM na podstawie stałych równania 
NRTL uzyskanych z korelacji LLE+VP2. 
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Rys.VIII-10. Zestawienie ciśnień eksperymentalnych z obliczonymi na pod- 
stawie korelacji LLE+VP1 równaniem NRTL. 
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Rys.VIII-11. Zestawienie ciśnień eksperymentalnych z obliczonymi na pod- 
stawie korelacji LLE+VP2 równaniem NRIL. 
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Јак powiedziano w rozdziale poprzednim, do korelacji zastosowano 
równiez uogólnione równanie LEMF 2 trzema dobieralnymi parametrami. Ko- 


relacje tym równaniem wymagają osobnego omówienia. 


Przydatność stosowania ujemnego parametru « do korelacji równowag 
ciecz-para i ciecz-ciecz była wielokrotnie analizowana w przeszłości 
(Marina i Tassios, 1973; Stryjek i wsp., 1981) najczęściej w kontekcie 
korelacji ze granicznym współczynnikiem aktywności. Ostatnio Bittrich i 
wsp. (1989) badali możliwości korelacyjne równań NRTL i ТЕМЕ w bezpo- 
średnim sąsiedztwie punktu krytycznego mieszalności. Stwierdzili oni, że 
lepsze jest równanie LEMF gdyż daje krzywą współistnienia opisywaną 
wykładnikiem krytycznym В o wartości nieklasycznej, w przeciwieństwie do 
równania NRTL, którje daje klasyczną wartość В=0.5. Badania swoje prowa- 
dzili dla zakresu temperatur zredukowanych (wzór IV-31) od 0 do e=10 *. 


Powyższe badania potwierdza zachowanie się parametrów w funkcji 
temperatury w korelacjach przeprowadzonych w tej pracy. O ile w przypad- 
ku nawet najlepszej korelacji LLE+VP2 równaniem NRTL obserwujemy silne 
zakrzywienie niektórych parametrów w bezpośrednim sąsiedztwie T (rysun- 
ki VIII-2 i VIII-3 ), to w tym przypadku parametry dają wykres prawie 
linii prostej ze słabym zakrzywieniem, niewidocznym w skali rysunku. Co 
więcej, korelacje LLE+VP1 i LLE+VP2 dają dokładnie takie same wartości 
parametrów. Rozbieżności występują dla parametrów wyliczonych z warunków 
stanu krytycznego. Powodem jest brak zbieżności dla tego punktu chociaż 
wartości funkcji celu (tablice VII-9 i VII-10) prawie spełniają 
kryterium wyjścia z iteracji. Te drobne różnice nie są istotne z punktu 
widzenia korelacji. Nierozróżnialność parametrów powoduje, że te korela- 


cje noszą dalej wspólną nazwę LLE+VP. 


Aby uzyskać pełniejszy obraz zachowania się parametrów uogólnionego 
równania LEMF, przeprowadzono dodatkowo korelacje przedstawionych w roz- 
dz. VII siedmiu izoterm równowag ciecz-para. Wyliczone parametry podane 
są w tablicy VIII-1. Mają one charakter liniowy, podobny do uzyskanego w 
korelacjach równaniem NRTL. Jakość dopasowania równania do danych VLE 
jest podobna (ostatnia kolumna tablicy). Parametry te w funkcji tempera- 


tury zaprezentowane są graficznie na rysunku VIII-12 (linie przerywane) 
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wraz z uzyskanymi z korelacji LLE+VP (linie ciągłe). Jak widać, ich za- 
leżność temperaturowa jest bardzo podobna. W przypadku korelacji LLE+VP 
współczynniki temperaturowe poszczególnych parametrów różnią się, lecz 
są tego samego rzędu. Tablica VIII-2 przedstawia równania przy pomocy 


których aproksymowano poszczególne parametry. 


Tablica VIII-1. 
Parametry równania LEMF uzyskane z korelacji VLE. 


. 735551 . 928410 
. 721544 . 921269 
. 709497 . 913489 


. 702354 . 907026 
. 696775 . 904949 
. 682602 . 898990 
. 664708 . 893731 


N @ со 000-0 


Obliczone preznošci рагу dla temperatur ponizej temperatury kry- 
tycznej i ich zgodność z uzyskanymi doświadczalnie nie różnią sie zasad- 
niczo od uzyskanych z korelacji LLE+VP2 równaniem NRTL. Widoczne to jest 
z porównania rysunków VIII-11 i VIII-13. W przypadku równania NRTL zgod- 


ność z doświadczeniem wydaje się nawet lepsza. 


Wobec prawie liniowej zależności stałych od temperatury spróbowano 
porównać zmierzonene prężności pary trzech izoterm powyżej UCST z obli- 
czonymi z parametrów uzyskanych na podstawie prostej ekstrapolacji do 
tych temperatur. Wynik przedstawiono na rysunku VIII-13 (linie przerywa- 


ne). 


Uzyskane rezultaty wykazują powiększanie się różnic pomiędzy 
ciśnieniem obliczonym i doświadczalnym, szczególnie dla większych stężeń 
3-MOPN. Jest to zrozumiałe, gdyż powyżej temperatury krytycznej 
następuje zamiana wartości parametrów energetycznych i trend wzrostu 
ciśnienia obserwowany w temperaturach poniżej Те dla tych stezeñ jest 
zachowany. W $wietle tego nie można spodziewać się uzyskania poprawne j 
równowagi ciecz-para w temperaturach wyższych od krytycznej, obliczonej 
na podstawie parametrów uzyskanych z równowagi ciecz-ciecz, nawet jeśli 


w tej korelacji użyje się dodatkowej danej w postaci prężności pary nad 
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roztworem dwufazowym. Podobnych rezultatów należy spodziewać się м przy- 


padku zastosowania granicznego współczynnika aktywności. 


1.50 


c12 


2 
Cn 
о 


= LLE+VP 
*— — VLE 
+++-++ ekstrapolowane 


Bezwymiarowe parametry LEMF 
° 
сл 
[e] 


-1.50 
295.00 300.00 305.00 310.00 315.00 320.00 325.00 
Temperatura (K) 


Rys.VIII-12. Temperaturowa zależność stałych równania ГЕМЕ. 


Tablica VIII-2. 


Parametry równania ГЕМЕ zastosowane obliczenia ciśnień powyżej tempera- 


tury krytycznej. 


Zależność temperaturowa parametrów LEMF 


1.37931 - 0.0085959*T/K 
-0.84550 + 0.0053415*T/K 
2.22276 - 0.0044731*T/K 


T/K 311.143 - T T 


313.160 -2.017 309.126| -1.27791 0. 805695 
318.150 -7.007 304.136| -1.23501 0. 779004 
323.270 -12.127 299.016| -1. 19100 0.751692 
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Rys.VIII-13. Porównanie eksparymentalnych prężności рагу z obliczonymi 
na podstawie ektrapolacji parametrów do wyzszych tempera- 
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Rys.VIII-14. Zależność ciśnienia od temperatury i składu obliczonego na 
podstawie parametrów uzyskanych z "lustrzanego odbicia" 
względem temperatury krytycznej parametrów z korelacji 
LLE+VP. 
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W zwiazku 2 powyzszym zalozono, Ze w temperaturze krytycznej wszys- 
tkie parametry przyjmują swoje ekstremalne wartości, а dalsze zwięk- 
szanie temperatury wymaga zmniejszania ich wartości. Przeprowadzono nas- 
tepujaca estymację parametrów dla trzech temperatur wyższych od krytycz- 
nej. Obliczono różnice pomiędzy Те (311. 143 K) а tymi temperaturami (ko- 
lumna 2 tablicy VIII-2), а nastepnie dodano је do Te uzyskujac nowe tem- 
peratury T' (kolumna 3 tabl.VIII-2),odpowiadajace temperaturom ponizej 
krytycznej. Dla nich obliczono parametry (oznaczone ' w tablicy VIII-2) 
stosując wyrażenia podane również w tablicy VIII-2 (na rysunku VIII-12 
przebieg nowych parametrów zaznaczony jest gęstą linią przerywaną). Z 
koleji, na podstawie tych parametrów obliczono ciśnienia dla postulowa- 
nych temperatur powyżej Te: Rezultaty obliczeń przedstawione sa na гу- 
sunku VIII-14. 


Widoczne jest znaczne zmniejszenie różnic pomiędzy ciśnieniami 
doświadczalnymi i obliczonymi w porównaniu z rysunkiem VIII-13. Ponadto, 
w miarę wzrostu temperatury, również widoczne jest zmniejszanie się tych 
różnic, co sugeruje zanikanie wpływu otoczenia punktu krytycznego na ра- 
rametry równania korelacyjnego, a co za tym idzie, na obliczone 
prężności pary. Na podstawie uzyskanych rezultatów obszar ten należy 
oceniać na około 20 stopni powyżej i poniżej temperatury krytycznej roz- 


puszczalności. 


Powrót do klasycznego zachowania się układu, w miarę wzrostu tem- 
peratury powyżej То przy zastosowaniu opisanej "ekstrapolacji" ilustru- 
je również temperaturowo-stezeniowa zależność swobodnej entalpii miesza- 
nia przedstawiona na rysunku VIII-15 jako linie przerywane. Widoczna tam 
jest tendencja do powrotu kształtu oraz wartości tej funkcji uzyskiwane j 


na podstawie korelacji VLE. 
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Rys.VIII-15. Zależność temperaturowa GM powyżej i poniżej temperatury 
krytycznej uzyskana z parametrów korelacji LLE+VP. 


Takie modelowanie równowag i diagramów fazowych wynika z zastosowa- 
nych przy konstruowaniu równań założeń teorii pola uśrednionego Van der 
Waalsa. Na podstawie teorii krzywa współistnienia dwóch faz ciekłych 
jest stopnia drugiego, a zależność ciśnienia od składu (gęstości) - sto- 
pnia trzeciego. Kilka lat temu ukazała się praca Singha i Van Hooka 
(1987a) na temat zakresu temperatur, w którym obserwowane są anomalie 
punktu krytycznego mieszalności. Stwierdzili oni, że dla temperatury 
zredukowanej e > 10 3 można z powodzeniem stosować wyrażenia klasyczne 
na G. (np. równanie R-K) do opisu krzywej wspótistnienia. Zakres ten od- 
powiada ok. 0.3 K w skali bezwzglednej. W $wietle wyników tej pracy wy- 
даје sie, Ze zakres ten jest duzo wiekszy i chyba nalezy wyrozróznié dwa 
zakresy temperatur: wiekszy, rzedu 20 K, nazwany przedkrytycznym i dru- 
gi, właśnie rzędu 0.3 К manifestujący sie anomaliami punktu krytycznego 
w sensie zjawisk opisanych w rozdz. IV i V. Na potwierdzenie powyższej 
hipotezy można przytoczyć np. zachowanie się stałej Kerra dla układu 
(alkohol benzylowy(1) + metylocykloheksan(2)) (Pyżuk i wsp. 1980), która 
odłożona w funkcji temperatury zaczyna odchylać się od normalnego zacho- 


wania już ok. 20 K przed temperaturą krytyczną. Podobnie zachowującą się 
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własnością jest lepkość (patrz rys. V-2). 


VIII-3 KSZTAŁT KRZYWEJ PRĘŻNOŚCI PARY DLA IZOTERMY KRYTYCZNEJ 


Z rozważań termodynamicznych wynika, że odpowiednikiem ciśnienia w 
układzie jednoskładnikowym w przypadku dwuskładnikowego układu ciekłego 
jest różnica pochodnych potencjału chemicznego po składzie i tejże 
pochodnej w punkcie krytycznym. Jej zachowanie było badane przez Simona 
i wsp. (1972) dla układu (CF, (1) + CH, (2) }. Potwierdzili oni, ze 
krzywizna jej przebiegu jest zgodna ze znalezioną dla ciśnienia 
w układach jednoskładnikowych. Wyznaczony przez nich wykładnik krytyczny 
G wynosi 4.6. Poza tym przypadkiem nie znaleziono podobnej analizy dla 
innego układu dwuskładnikowego. 


Z drugiej strony, kilkadziesiąt lat temu Rice (1950) postawił hipo- 
tezę, że przebieg ciśnienia cząstkowego jednego ze składników w układzie 
dwuskładnikowym w temperaturze krytycznej mieszalności powinien mieć 
analogiczny prze 
bieg jak izoterma krytyczna w układzie jednoskładnikowym. W tamtej 
pracy nacisk położony był na zagadnienie, czy krzywa współistnienia ma 
płaski odcinek, taki jaki występuje w układach jednoskładnikowych. Po- 
nieważ nie jest powiedziane, którego ze składników dotyczy ta hippoteza, 
należy domniemać, że takie samo zachowanie będzie wykazywało ciśnienia 
całkowite nad roztworem. Ta praca dostarczyła odpowiedniego materiału 
doświadczalnego pozwalającego potwierdzić hipotezę Rice'a. Na rysunku 
VIII-16 przedstawiono obliczone ciśnienie gdy zastosowano następującą 
zaleźność: 


|P- P | = К. |х = x |? VIII-1 


gdzie K jest stałą proporcjonalności. Zależność ta jest słuszna dla 
układów jednoskładnikowych, co jest przedstawione w rozdziale IV. Jak to 
widać na rysunku odtworzenie ciśnienia całkowitego jest całkiem zadowa- 
lające dla wykładnika krytycznego 6-5, a więc zgodnego z tzw. niekla- 
sycznym, natomiast użycie wykładnika klasycznego 8=3 nie spełnia 
zależności VIII-1. 
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Rys.VIII-16. Porównanie eksperymentalnego ci$nienia zredukowanego do 
obliczonego z wykładnikiem 8=5 i 6-3. 


200.00 
Uklad: СЕ4(1) + CH4(2); Т=98.013 K 


***** eksperyment 


©150.00 — — delta=3 
工 delta=5 
Е delta=4.6 
ABS(P — P.) = K * ABS(x — x) 

O 
к 100.00 
B. 
Л 
m 
< 

50.00 

0.00 Е 

0.00 0.20 0.40 0.60 0.80 1.00 


х(1) 


Вуз. У111-17. Porównanie ciśnienia zredukowanego z obliczonym z wykład- 
nikiem 8=5, 4.6 i 3 dla układu (CF,(1) + CH, (2)). 
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Aby potwierdzić zależność stężeniową ciśnienia całkowitego w tem- 
peraturze krytycznej mieszalności podobne obliczenia przeprowadzono dla 
układu (CF4(1) + CH, (2)) i uzyskano potwierdzenie dla układu o zupełnie 
odmiennym charakterze. Co więcej, użycie wykładnika 6=4.6 wyznaczonego 
przez autorów dało jeszcze lepszy opis ciśnienia. Oczywiście wartości 
stałej proporcjonalności są dla każdego układu inne. Niestety brak одро- 
wiednich danych dla innych układów nie pozwala w pełni potwierdzić uni- 
wersalnego charakteru takiego zachowania, ale w świetle tego, co wiadome 
jest na temat punktu krytycznego należy przypuszczać, że dla innych 
układów dwuskładnikowych ciśnienie całkowite będzie miało ten sam prze- 
bieg. 


VII-3. WPŁYW POPRAWKI OBJĘTOŚCIOWEJ NA KORELACJE IZOTERMICZNYCH RÓWNO- 
WAG CIECZ-PARA 


W wyprowadzeniach wyrażeń na nadmiarową energię Gibbsa zakłada się, 
że układ nie wykazuje nadmiarowej objętości mieszania (VF=0). Pozwala to 
na przyrównanie GÉ do nadmiarowej funkcji Helmholtza. Pełne równanie 


Gibbsa-Duhema redukuje sie dla pomiarów izotermicznych do postaci: 


dP = ) x, dln 7, VIII-2 


Wyniki niniejszej pracy pozwalaja obliczyé роргамке do pochodne ј GE 
ро składzie jaką wnosi występowanie УЕ w roztworze rzeczywistym. 
Obliczenia te są zestawione w tablicy VIII-3. 


Jak powiedziano wcześniej, stosując metodę Tao konieczne jest poli- 
czenia pochodnej ciśnienia całkowitego względem składu. Wobec dużej 
liczby punktów doświadczalnych w izotermie i stosowaniu wygładzania да- 
nych ciśnienia za pomocą funkcji opisującej "nadmiarowe" ciśnienie, po- 
danej w rozdziale VII (wzór VII-10) należy przypuszczać, że policzone na 
podstawie tej funkcji pochodne dla temperatury krytycznej mieszalności 
(kolumna 3 tablicy VIII-3) mają wartości zbliżone do rzeczywistych. War- 
tości VE dla tej temperatury zostaty obliczone przez interpolacje za ро- 
moca równania (3) w pracy na temat VE (Luszczyk， 1988). Wspótczynniki 
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aktywnosci obliczone zostaty na podstawie stalych równania МЕТ. dla tej 
izotermy (311.146 K) zawartych w tablicy VII-2d. Wartość stałej gazowej 
dla stosowanych jednostek wynosi 8.3141-10 [cm -kPa-K '-mol 1]. 


Tablica VIII-3. 


Poprawka do pochodne j G" zwiazana z wystepowaniem VE obliczona dla izo- 
termy 311.146 K uktadu (3-MOPN(1) + woda(2)). 


юж лл юго 


. 0561332 
. 2666156 
. 7011618 
. 2545791 
· 1337569 
‚5036980 
. 8856307 
. 3085034 
. 8066959 


= m кә кэ => NN b 00 


OPO N [МЮ = PR M M 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 


Jak widaé poprawka nadmiaru objetosci (ostatnia kolumna tablicy 
VIII-3) dla rozpatrywanego układu, jak powszechnie jest uznawane, nie ma 
praktycznego znaczenia na wartość pochodne j e, W układzie rzeczywistym 
wartość tej pochodnej jest blisko liniową funkcją stężenia (testy kon- 
systencji) o względnie dużych wartościach w stosunku do obliczonej po- 
prawki. Jej wpływ może być znaczący dla układów o dużej różnicy w 
prężnościach par czystych składników i jednocześnie posiadających duży 
bezwzględny УЕ. Ponadto, w realnych układach wartości pochodnych 
ciśnienia są największe dla stężeń skrajnych, w których z kolei wartości 
vt są najmniejsze. W rezultacie następuje znoszenie się tych dwóch 


wpływów, co jest widoczne w tablicy VIII-3. 
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X ZESTAWIENIE WYNIKÓW 


Zbudowano i sprawdzono stanowisko do pomiarów równowagi ciecz-para 


metodą statyczną w obszarze ciśnień od 0.1 do 100 kPa. 


Zaprojektowano i zbudowano orginalne naczynie równowagowe z jednocze- 
snym mieszaniem fazy ciekłej i parowej co znacznie przyspiesza uzys- 
kanie stanu równowagi termodynamicznej. Naczynie wyposażone jest w 
chłodnicę zabezpieczającą przed destylacją składników na powierzchnię 


rtęci w manometrze podczas zmian temperatury. 


Zaprojektowano i zbudowano prosty i sprawny system przygotowania pró- 


bek o znanej masie i ich wprowadzania do aparatury. 


Zmierzono prężności pary dla substancji czystych oraz dla serii 31 
próbek o różnych składach mieszanin 3-MOPN z wodą. Badany układ сћа- 
rakteryzuje się dużą różnicą prężności par czystych składników i 


tworzy dwie fazy ciekłe w temperaturze poniżej 311.143 K. 


Zaproponowano równanie do opisu krzywej współistnienia, pozwalające 
wygładzać dane zarówno równowagowe jak i rozpuszczalności. Sprawdzono 


przydatność tego równania. 


Zmierzono 16 izoterm nadmiarowej objętości mieszania badanego układu 


w zakresie temperatur od 293 do 313 К. 


Przeprowadzono dyskusję zakresu temperatur, w którym zaczynają mani- 
festować się tzw. "zjawiska przedkrytyczne" i anomalie stanu krytycz- 
nego. Stwierdzono, że zjawiska te występują ok. 20 K powyżej i poni- 


żej krytycznej temperatury mieszalności. 


. Wykazano, że zależność ciśnienia całkowitego od stężenia jest stopnia 


piątego (tzw. wykładnik krytyczny 8) również dla układów dwuskład- 


nikowych w temperaturze krytycznej mieszalności. 
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9. Przeprowadzono dyskusję wpływu poprawki objętościowej na korelacje 


izotermicznych pomiarów równowagi ciecz para. 
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ZALACZNIK A 


Model Isinga 


W 1925 roku Ising zaproponowat model do rozwiazania problemu ferro- 
magnetyzmu. Zatozenia i sposób rozwiazania jest dla nas o tyle interesu- 
jacy, Ze moga byé przeniesione na inne zagadnienia zwiazane z wystepo- 
waniem stanów alternatywnych (Smirnowa, 1980). W zagadnieniach tych spo- 
sób rozumowania, poszukiwane wielkošci, terminologia sa podobne i w 


zwiazku 2 tym przedstawiono go w tym zalaczniku w zarysie. 


Energia elektronowa zwiazana 2 orientacja spinów w weztach sieci 

może być zapisana w postaci: 
Ее I же A-1 
как 
1 <k 

gdzie zmienna s, moze przyjmowaé warosci +1, w zalezno$ci od orientacji 
spinu м i-tym węźle, powiedzmy, "w górę" lub "w dół". J jest tzw. całką 
wymiany i jeśli przyjmiemy, że jest jednakowa dla wszystkicj węzłów, a 


jej oddziaływanie maleje szybko wraz z odległością to energia wynosi: 


E = -Ј 55 А-2. 


Stosując przybliżenie Bragga-Williamsa (1936) można wykazać, że wystepo- 


wanie domen magnetycznych jest przemianą drugiego rodza ји. 


Jeśli przez NĄ oznaczymy liczbę spinów, dla których s = +1, a przez 


NB - liczbę Spinow о wartości B -1, to całkowita liczba węzłów w 


układzie wynosi: 


N=N+N A-3 
A В 
Termodynamiczne prawdopodobieństwo pozwala nam obliczyć entropię układu 
spinów: 
N! 
јан к ME z 
A B 
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N N 
S-kinW--k [М —А+м £ A-5 
A N B N 


Parametr uporzadkowania jest zdefiniowany nastepujaco: 
т = M/M = (N - МУ A-6 


Przy catkowitym uporzadkowaniu, czyli przy ustawieniu wszystkich spinów 
w jedną stronę m = 1 natomiast dla całkowitego nieporzadku czyli dla N= 
N, parametr ten wynosi zero. Stad N /N = 1/2(1 + m) i N/N = 1/2(1 - m). 


Entropia wyrażona poprzez parametr uporządkowania wynosi więc: 


S(n) = -kN Ë (1+n) 1n [(1+n)/2] + Я (1-7) 1n La-921] A-7 


Widać stad, że wraz z zanikiem uporządkowania (m > 0) entropia układu 


rośnie. 


Energia układu spinów wynosi: 


Е=-Ј (СМ +N_-N ) = - JNz/2 A-8 
AA BB АВ 


gdzie N oznacza liczbę najbliższych sąsiadów о wartości spinu I lub J, 
a Z jest liczba koordynacyjna sieci. Liczby zgodnie i przeciwnie skiero- 
wanych spinów wynosza odpowiednio: 


н, = ZN ( 1+1)2/8 А-9 
一 = 2 СӘ 

М в = ZN (1-0) /8 A-10 
— per 2 p 

Np = ZN (1-n)7/4 A-11 


Jeśli w rozpatrywanym węźle znajduje się spin o pewnym kierunku, to w 
przypadku ferromagnetyku, istnieje duże prawdopodobieństwo, że sąsiedni 
z nim spin będzie skierowany w tym samym kierunku. Jest to tzw. uporząd- 
kowanie bliskiego zasięgu. Metoda Bragga-Williamsa uwzględnia tylko da- 
lekie uporządkowanie i stąd parametr m nazywa się parametrem dalekiego 


uporządkowania. 


Energia swobodna związana z daną konfiguracją spinów, wyrażona jako 
funkcja parametru uporządkowania wynosi: 
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f(L) 


f(L)=tgh(7EL/KT) 


Rys.A-1. Rozwiązanie graficzne modelu Isinga. 


F(n) = Eln) - TS(n) = 
-ŻNZJN + NkT Ë (1*3) ln [(1+n)/2] + А (1-n) ln сати) A-12 


W stanie równowagi statystycznej (дЕ/дп) = 0. Poprzez zrózniczkowanie 


otrzymujemy: 
zJ 1 1+m - 
ЕТ = 5 ln ‚=ч A-13 
а zJ Е 
"= teh [7] A-14 
zJ 
Podstawiajac | КТ )" = x otrzymujemy: 
КТ Е 
т * tgh(x) A-15 


Rozwiazanie przedstawia rysunek A-1. Jak powiedziano, dla temperatury 
krytycznej wartość parametru uporządkowania wynosi 0, a więc wartość 


tangensa hiperbolicznego wynosi 1. Stąd: 


T= zJ/k A-16 


Jak widać, temperatura krytyczna jest wyrażona za pomocą całki wymiany. 
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(b) Magnetyk 


Rys.A-2. Krzywa wspótistnienia dla modelu siatkowego (a) i modelu Isinga 
(b). 


Krzywa namagnesowania wyrażona w postaci parametru uporządkowania w fun- 


kcji temperatury zredukowanej ma uniwersalny charakter. 


Analogia pomiędzy przedstawionym wyżej rozwiązaniem modelu Isinga a 
modelem siatkowym cieczy Guggenheima (1952) jest oczywista. Rozwiązanie, 
określające warunki powstawania krzywej współistnienia 1 punktu 
krytycznego jest takie same w obu modelach (wzór 66, Buchowski, 1961) po 


odpowiednim podstawieniu zmiennej x. Podobieństwo to ilustruje rys.B-2. 


Model Barkera-Focka 


Historycznie, pierwszym, który podał model termodynamiczny przewi- 
dujący istnienie luki zamkniętej, a więc i dolnego punktu krytycznego 
byli Barker i Fock (1953). Założyli oni istnienie dwóch rodzjów oddzia- 
ływań międzycząsteczkowych w układzie dwuskładnikowym. Cząsteczki roz- 
tworu umieścili w węzłach sieci o liczbie koordynacyjnej z=6. Każda 
cząsteczka posiada z kontaktów z czego jeden oddziałuje z sąsiednią mo- 
lekułą w inny sposób niż pozostałe (z-1) kontaktów. Jeśli cząsteczki 
różnego typu oznaczymy odpowiednio A i B, a kontakty odpowiednio Ap i Aq 
oraz Bp i Ва to zaproponowane przez autora oddziaływania i odpowiada jące 


im energie przedstawione są w tablicy A-1. 
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Tablica А-1 


Rozkład i rodzaj energii w modelu Barkera. 


Oddziaływanie Energia 


Ap-Ap, Ap-Aq, Aq-Aq 


Bp-Bp, Bp-Bq, Bq-Bq 0 
Ap-Bp U >0 (odpychanie) 
Ap-Bq, Aq-Bp, Aq-Bq U <0 (przyciaganie) 


Potencjaty chemiczne sktadników wyrazone sa w postaci: 
ш RT [(z-1)1n(x, /х,х, |) + ln(x Xx, а) + т xl A-17 


gdzie i - A,B. Parametry Xip i Xia znajdowane są w wyniku rozwiązania 
układów równań uzyskanych przez Barkera (1952) metodą przybliżenia 
quasi- chemicznego. Ich forma matemetyczna jest dość skomplikowana i nie 


jest omawiana. 


Z warunków równowagi termodynamicznej dla dwóch współistniejących 
faz: 


и = и, oraz u, = u , A-18 


gdzie subskrypty oznaczaja poszczególne fazy, oraz warunku stanu kry- 


tycznego: 


d^F/dx? = dsFvdx3 = 0 A-19 


gdzie F jest swobodna energia mieszania. 


Autorzy, dla warunków symetrycznych tzn. x = 1⁄2, uzyskali linie 
krytyczna јако funkcje zastosowanych energii 9, 1 U,. Jeśli temperature 
wyrazimy w postaci U,/KT, natomiast stosunek energii przyciągającej i 
odpychającej w postaci parametru m = |], to ich wzajemna zależność 
daje linię krytyczną przedstawioną na rysunku A-3(c). Rozmiar krzywej 
jak i jej kształt jest również zależny od liczby koordynacyjnej sieci co 
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Rys.A-3. Zamkniete luki mieszalnosci obliczone modelem Barkera-Focka dla 


2 = 6 (a) 1 2=4 (b), oraz linia krytyczna (c) dla 2=6. 


również przedstawione jest na rysunku. Autorzy również podali kształt i 
rozmiar luki mieszalności w przestrzeni T-x w zależności od wybranego 
parametru m. Cięcia lini krytycznej i odpowiadające im luki pokazuje 


rys. A-3. 
Model Wheelera i Andersena 

Każda cząsteczka jest tak pomyślana, że posiada jeden wyróżniony 
punkt kontaktowy A, i B i (w - 1) zwykłych punktów styku A, i B。 bez 


wzgledu сту znajduje sie w komórce sieci podstawowej czy wtórnej. Od- 


158 


działywania występują tylko pomiędzy cząsteczkami pierwszego i drugiego 
rzędu. Przyciąganie, czyli ujemna energia oddziaływania występuje tylko 
pomiędzy różnymi cząsteczkami w obu rodzajach komórek, lecz tylko gdy 
możliwe jest oddziaływanie wyróżnione dla cząsteczki drugiego rzędu. 
Przedstawia to tablica D-1. 


Tablica A-2 


Definicje aktywności energetycznych molekuł w modelu Wheelera i Ander- 
sena 


odpychanie 


przyciąganie 


Jeżeli w komórce drugiego rzędu znajduje się cząsteczka A lub B i 
posiada najbliższych sąsiadów w dwóch komórkach pierwszego rzędu zarówno 


A jak i B, wtedy oddzialuje przez specyficzny kontakt energią U nato- 


2* 


miast pozostale kontakty energia U, czyli odpychajaco. Gdy najblizszymi 


sąsiadami sa cząsteczki innego жЕ rozwazana molekuta, to wtedy po- 
siada ona dwa kontakty o energii U,- U,» oraz (w - 2) kontaktów o ener- 
gii Е Należy zwrócić uwagę, ze energia oddziaływania pomiędzy 
cząsteczkami w komórkach pierwszorzędowych jest niezależna od orientacji 
cząsteczek. Istotne jest tylko to, czy jest to stan A czy B. Takie upro- 
szczenie pozwala na odniesienie bezpośrednie do modelu Isinga, a przez 
to do rozwiązania zagadnienia jak dla trójwymiarowego modelu Isinga dla 
ferromagnetyków. Funkcja podziału jest odpowiednio przekształcona i ak- 
tywności składników i temperatura T spełniają rolę pola magnetycznego i 
temperatury Ti odpowiednio. Równowazne wielkošci rozpatrywanego modelu i 


modelu odniesienia wyrażają sie następująco: 
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Rys.A-4. Podstawowa sieé modelu Wheelera i Andersena (a) oraz obraz 


sieci I i II rzedu (b); obok wyglad komorki II rzedu (c). 


Co + y qe -2H E 
А [m+ (о - 1)т, 1 (1 + č) а exo[ -2J | "T 
[(o + СУ) (бо + y)11⁄2 x 


gdzie H jest polem magnetycznym, J - całką nakładania modelu Isinga, С 
jest zdefiniowane jako stosunek aktywno$ci sladnika A i B: 


С = ZA p A-22 
natomiast 
一 p 2 — 
у = 211+ (о 2)n, A-23 
przy czym 
n7 exp(-U, /kT) i n= exp(-U /kT). A-24 


Przej$cie fazowe w modelu isinga ma miejsce przy zerowym polu mag- 
netycznym H = 0. Wielkość A odgrywa role pola magnetycznego i musi po- 
siadać wartość 1. Dla tego warunku wyrazeni па £ przyjmuje postać: 


160 


2.0 


о 0.2 0.4 0.6 0.8 
(kT/U,) 


Rys.A-5. Zależność parametru £ od temperatury. 


1 + (е - 1) 1,7, 
Ё = 21 — А-25 
а gł 2nm+ (о - 2)т2 
12 1 
Dla prostej, trójwymiarowej sieci stwierdzono, Ze istnieje pewna kry- 
tyczna wartość parametru S = 0.64183, ponizej której nastepuje rozpad na 
dwie fazy. Ta niższa wartość występuje w pewnym zakresie temperatur dla 
odpowiednio dobranych wartości energii U, 1 U,. Przez analogie do magne- 
tyzacji jako parametru uporządkowania w modelu Isinga, składy równowago- 


we oblicza się wykorzystując jej nieciągłość. 


m (2) 
1 o 1 о > 
<= ж; t у а [1 * 24577) А-26 
gdzie 
x= NĄ/ (NA+ Ng): А-27 


Model Ма1кега i Vause 


W tym modelu cząsteczki umieszczone sa w d wymiarowej sieci. Analo- 


gicznie do modelu Isinga zmienna Si przyjmuje wartości +1, w zależności 
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od tego jakiego rodzaju molekuła znajduje sie м i-tym węźle. Ponadto 
każdą cząsteczkę charakteryzuje druga zmienna o, która może przyjmować 
wartości 1,2,3...q. Zmienna ta informuje o orientacji cząsteczki. Wynika 
stąd, że cząsteczka może mieć q wyróżnionych kierunków oddziaływania. 
Cząsteczki oddziałują ze sobą poprzez następujący hamiltonian: 


Я = (1755s) IE + (K- K)8 ] A-28 


> 
стр 2 сїт] 


Кин 


gdzie 8 jest delta Kroneckera przyjmującą wartość 1 dla jednakowych war- 
tości vi i O dla pozostałych. Taka konstrukcja hamiltonianu pozwala 
określić w prosty sposób warynki oddziaływania najbliższych sąsiadów w 
sieci. Jeśli sąsiednie cząsteczki są tego samego typu, wartość czynnika 
(1 - 8,3, ) wynosi O i cząsteczki nie oddziałują. W przypadku różnych mo- 
lekuł, elcz gdy ich względne orientacje są różne, co wyraża się różnymi 
warościami zmiennej c, cząsteczki oddziałują ze sobą siłą odpycha jącą 
К. Natomiast w przypadku różnych warości zmiennej spinowej i takiego 


2 
samego stanu wewnętrznego molekuł oddziaływanie jest przyciąga jące K,. 


Model ten jest rozwiązany pod względem mechaniki statystycznej 
przybliżeniem Migdal-Kadanoff'a z zastosowaniem metody renormalizacji 
grup, w której współczynnik przeskalowania sieci wynosi 3. Uzyskany dia- 


gram fazowy w przestrzeni parametrów pokazany jest na rysunku А-6. 


Obszar, którego centalnym punktem jest HT obejmuje układ jednofazo- 
wy czyli stan całkowitego wymieszania. Na lewo od krzywej AB znajduje 
się obszar ciekły dwufazowy, najbardziej dla nas interesujący. Przedsta- 
wiony diagram fazowy pozwala określić występowanie przejścia fazowego 
dla pewnych stałych watrości parametrów K, i K>, które sa scharakteryzo- 


1 
wane przy pomocy parametru R = -K,/K,. 
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Rys.A-6. Diagram fazowy w przestrzenie parametrów modelu Walkera-Vause. 


Punkty charakterystyczne oznacza ја: 

HT - punkt wysokotemperaturowy; początek układu odpowiadający temepratu- 
rze nieskończone j, 

- punkt krytyczny przejścia ferromagnetycznego 

- jak wyżej dla antyferromagnetyków 

2q-stanowy model Pott'a 

- punkt multikrytyczny 

- górny punkt krytyczny mieszalności 


КЕ jj ЗА 
1 


- dolny punkt krytyczny mieszalnosci 


Przedstawiony hamiltonian (А-28) ma pewne ogóle własności: 


1. К = Ко“ K 


W tym przypadku redukuje sie on do: 
Я = 1⁄2 K >' (1 - 5,5, ) А-29 


czyli sprowadza sie do modelu Isinga, omówionego wcześniej. 
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K = 1⁄2 K, >'(1 - ss )ó A-30 
1 1j isj 


W punkcie P nastepuje krystalizacja mieszaniny czasteczek A i B 


3. Куж LAE lecz К< 0 
Ten przypadek wykorzystuje pełną postać hamiltonianu. Określa on miejsce 
na diagramie fazowym, w którym występuje i zanika przejście fazowe 


ciecz-ciecz. 


Przedstawiona na rys.A-6 linia kropkowana (oznaczająca wspomnianą, 
stałą wartość parametru R) wyznacza na lini krytycznej  A-F-B 
występowanie dolnego i górnego punktu krytycznego mieszalności. Jej 
przecięcie się z linią D-P-Q oznacza wejście w obszar występowania mie- 


szaniny cząsteczek A i B w fazie stałej. 


Rozwiązanie tego, czyli znalezienie współistniejących składów 
sprowadza sią, podobnie jak w modelu Wheelera i Andersena do przypisania 
wielkości obliczonych za pomocą hamiltonianu wielkościom 


charakterystycznym (magnetyzacji) modelu Isinga. 


+ 
ш iR A-31 
2 
(8q? + 20) 
Am(K , K,) Am (K^ Ке) A-32 
p dd 1*^2 


Czynnik 2, oraz pochodne wielkości energetyczne oblicza się przy pomocy 
wzorów rekurencyjnych podanych przez autorów. q pełni rolę liczby 


koorrdynacyjnej, jak to powiedziano w rozdz. IV-4.2. 
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ZALACZNIK B 
Metoda Tao 
Metoda  obliczania  równowagi  ciecz-para  stosujaca  numeryczne 
calkowanie równania Gibbsa-Duhema dla procesu izobaryczno-izotermicznego 
nosi nazwę metody Tao (1961). Próba zastosowania jej w tej pracy 
wynikała z tego, że nie wymaga ona ztosowania konkretnego wyrażenia na 


u^ a przez to, konkretnego modelu. 


Dla układu dwuskładnikowego równanie 111-14 może być zapisane nas- 


o o 
Supp Xp. 1251; a B-1 
AE 72 ФР 7; 


Równanie to zawiera tylko dwie niewiadome 7, i T, Z drugiej strony wia- 


tępująco: 


domo, że: 
LGE/RT =x lny +x lny B-2 
1 1 2 2 
n+1 
x 
E E E f 7, 
G /RT| n-G /RT| (n-1) = AG /ВТ = | [ ln — +r Jax, B-3 
*, "n n 72 
х 
1 
V*apP 
gdzie Г = dla warunków izotermicznych; subskrypt n oznacza, że war- 


RTdx 
tość obliczana jest w punkcie n o składzie x. Stosując całkowanie nume- 


ryczne metodą trapezów uzysku jemy 


7, 7, 
e(n-1⁄2)1n|—| = (177) | + €(n-1/2) = = 
2 2 
ne (п-1)є 


та | (п-1)е* (T| )(e/2) B-4 


пе“ T (а-у 


Odpowiednia transformacja i sumowanie prowadzi do wzorów na składniki 
równania  b-4, które rozwiązuje się ze względu na stosunek 


współczynników aktywności. Dla warunków startowych x = 0, T = 0, 1 (Ра 
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(п-1/2)є 
9 7, 
== = СА |== + СВ B-5 
12 е nn 
AL 
c= C „јех Се“ Г (р–1 је) /2] y B-6 
2 
(n-1)e 
ne Р? (1-пє) P 
AP ar B5 — ar uis 


Obliczenia rozpoczyna sie od wygładzenia danych prężności рагу w 
funkcji składu za pomocą np wielomianu. Umożliwia to interpolację 
wektora n prężności dla ne składów konieczną do numerycznego całkowania. 
Należy również wyestymować wartość granicznego współczynnika pierwszego 
składnika . Tao w swej pracy zaproponował wyrażenia podanego przez Gau- 


treaux i Goatesa (1955). 
1 
ES ] - 


PE ЕЗ 
1 p? Ox т| 
1 
Metoda ta została zaproponowana 30 lat temu, lecz nie zyskała рори- 


x со 
1 


larności wobec szybkiego rozwoju równań korelacyjnych i metod oblicze- 
niowych z jednej strony oraz jej niedogodności - z drugiej. Wymienione 
niedogodności to konieczność aproksymacji zmierzonego ciśnienia względem 
składu za pomocą jakiegoś wyrażenia matematycznego, konieczność 
estymacji granicznego współczynnika aktywności któregoś ze składników na 
podstawie danych dość dalekich od stężeń granicznych. Ponadto, obecność 
punktu azeotropowego powoduje rozbieżność w trakcie optymalizacji, co 


wymaga sztucznych zabiegów w celu przywrócenia zbieżności. 
Metoda MGC 


Metoda Mixona-Gumowskiego-Carpentera (Mixon i wsp., 1965), zwana 
dalej metodą MGC, powstała m.in. w calu, jak piszą autorzy, umożliwieniu 
obliczeń układów trój i więcej składnikowych bez konieczności stosowania 
explicite wyrażenia na с". Funkcja ta jest przyblizana za pomoca wyrazen 
111-1 i 111-2 dla przypadku ogólnego (układ wieloskładnikowy) . Wzór na 


ciśnienie całkowite nad roztworem N składnikowym przyjmuje postać: 
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ас = aG 


N-1 N-1 

P WSE + ax ча, 十 x Pex (6 - 3x | B-9 
gdzie G = G'/RT. 
Równanie to jest rozwiązywane, podobnie jak w metodzie Tao, przez 
obliczenie,w pierwszym etapie,wektora ciśnień nad roztworem z interwałem 
e, dla składów RM ne. Warto$ci pochodnych за przyblizane pochodnymi 
numerycznymi: 


G(ne + е) - G(ne -е) 
= = == —— A B-10 


дх 2є 
Wartości w nawiasach stanowią indeksy wektora funkcji G(n). Podstawiając 
B-10 do B-9 otrzymujemy: 


N-1 G(ne + е) - G(ne -е) N-1 С(пє + е) - ша a 


— 0 MA. ——— I ——Á 
Р(пе ) = posteo (s * 22 ne, ŻE 


im 
| 


B-11 


N-1 G(ne + е) - G(ne, -е) 
1 де | 


0 
+ x Pex [5 pu 


Wartość С jest obliczane iteracyjnie i n-ta poprawka jej estymaty g^ 


obliczana jest ze wzoru wiążącego  n-te przybliżenie ciśnienia 


całkowitego wzorem: 


(n) id (n) id 
p 


p p 
- 2—2 g(k-1) + P™ (k)g(k) + ———— g(k-1) 
= P(k) - P™ (k) В-12 


Wektor poprawek g z rozwiązania układu M równań tworzących macierz trój- 


przekątniową, której elementy są zdefiniowane następująco: 


о — 
Rozwiązanie układu równań: 
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(n) 


(n) (n) 

в. uy. © 0 0 g(1) PORE 

(n) (n) (n) (n) 

А, B e a 0 0 g(2) P(2) P esi 

(n) (n) (n) = _pm z 

0 А а) B a) Сүз) 0 g(3)| = |P(3) Fay B-14 
(n) (n) (n) 
А w) B d g(M) P(M) Р 


даје wektor g(i), który додаје sie do aktualnego wektora G(i). Iteracje 


przerywa się, gdy suma elementów wektora PP nie zmienie sie o 


określoną wartość. 


Metoda ta posiada podobne ograniczenia co metoda Tao aczkolwiek 
jest często stosowana przez grupę B.D.Smitha do układów wykazujących 
niewielkie odchylenia od doskonałości (np. Khurma i wsp.; 1983). 
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ZALACZNIK С 


W niniejszym zataczniku przedstawione sa wzory na pierwsza, druga, 
trzecią i czwartą pochodną względem składu wyrażeń na nadmiarowa energie 
Gibbsa stosowanych w tej pracy. 


aid Е: q'4/RT 
дх, 
ах = ах = 1, -1 i, ј=1,2 
1 
Q = Y x Ја x 


k dx k-1 
DAJĘ | k = 2,3,4 
k x 
дх i 


Q = СЕЛАТ 


© 
I 


ха - x) P A.G = 44) 5H 
i 


= E 1 3 w E 2 (1-2) 
ах = } A Y деј ерй =) ма 1)(1 - 2x) 


2 
ES E | JAA - еј M 
1 


+ } А|(ї-1)(1 - żyj 4 


- 2x(1 - x) ) А, (1-1)(1-2)(1 - zo 
i 
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3 

ô Q 4 | 7 Аүї(ї-1)(1 " 2x) 079 " 

i 

+ У A,(1-1)(1-2)(1 - 2017 + 
i 


"ЈА (i-1)2(1 - у“ + 
T 1 


- 2x(1 - x) Y A,a-1)a-2)0-3) - zam 
H 


ato (4-3) 
— = -8 4 Y A 1(1-1)(i-2) (1 - 2х) + 
4 1 
дх 1 
* Y A (1-1)(4-2)01-3) = 2x) 079 + 
1 
+) A,(1-1)(1-2)7(1 - 297 + 
i 
+ Y A a-0*a-20 - zo? + 
i 
=2х(1 = x) } A GHZ - Xd 
i 
NRTL 
Q = G'/RT 
ЕТ = e BT ur ехр[-ат, |] 
$1 = ) хт G $2 = Y xG 
JAJA EE 
S1 
„ję ) Хи SZ 
53 = ) dx T G 54 = ) dx G 
„ии p 
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20 _ 51 53 _ 51:54 
=) [ах Š + x, ip. )| 


2 2 
99 2 ах 53 _ 51-54 T 53.54 _ 51-54 
1| 52 2 i 


ах“ 7 52 s2? 523 

3 2 2 3 
29-6) | ах |- а BA | - х (- 83-84 ‚ 51 ч JI 
ax S2 S2 S2 S2 


3 4 


4 2 3 3 4 
29 ИСЕ ) - (82:8 _ 51:54 )] 
: 52 52 


UNIQUAC 
Q = G/RT 
T „м Ча 
$, z 8, 
е Пар F j хл. + 2 L q,x,1n ry 
1 1 1 
Ч = = ) an e] 
XT x q 
* ы 9, š 一 
X 
) jj ја 
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ZALACZNIK D 


Preznošci pary uktadu (benzen(1) + cykloheksan(2)) w funkcji tempe- 
ratury dla poszczególnych próbek. Kolumna trzecia zawiera odchylenia od 
równania Antoine'a, którego wspólczynniki wraz 2 odchyleniem standardo- 
wym podane sa w tablicach. Punkty oznaczone gwiazdką zostały usunięte z 
korelacji. 
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Tablica D-1. 


Cykloheksan 


bc0968 bc1342 


A=13.852715 
B=2828. 446 
C=47. 529 
RMSD=0. 009 kPa 


. 366 -. 
. 374 
. 380 
. 053 
. 051 -. 
. 053 
. 260 -. 
. 262. 
. 266 
. 129 
. 131 
. 133 
. 677 
. 678 
. 679 
. 081* 
. 069 
. 064 
. 318 
«307 =. 
.307 - 
637“ 
. 623 
. 627 
. 984 
.982 
.990 
‚ 665 
. 670 
. 680 
‚701 
‚702 


009 


. 004 
. 003 
‚001 


002 


. 000 


001 


. 002 
. 005 
. 005 
. 006 
. 009 
. 006 
. 008 
. 008 
. 017 
. 006 
. 001 
. 008 


001 


. 003 
. 017 
. 003 
. 007 
. 003 
. 005 
. 004 
.014 
. 007 
. 005 
. 001 
. 006 


А=13. 984061 
B-2893. 883 
С=43. 650 
RMSD=0. 006 kPa 


A=14. 093450 
B=2950. 982 
C=40. 788 
RMSD=0. 005 kPa 
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Tablica D-1. (c.d.) 


А=14. 110607 А=14. 236172 А=14. 299423 

В=2959. 034 B=3021. 716 B=3050. 100 
=40. 328 С=37. 179 С=35. 801 

RMSD=0. 005 kPa RMSD=0. 004 kPa RMSD=0. 004 kPa 
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Tablica D-1. (c.d.) 


bc4222 bc4730 bc5114 


А=14. 330348 А=14. 323835 А=14. 349177 
В=3062. 973 В=3057. 272 В=3069. 584 
С=35. 245 С=35. 550 С=35. 066 
RMSD=0. 006 kPa RMSD=0. 004 kPa RMSD=0. 005 kPa 
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Tablica D-1. (c.d.) 


bc5502 bc6484 bc7536 


A=14. 295028 A=14. 299545 A-14. 228745 
B=3037. 837 B=3037.213 B=2995. 389 
C=36. 603 C=36. 950 C=39. 396 
RMSD=0. 003 kPa RMSD=0. 005 kPa RMSD=0. 005 kPa 


T/K P(kPa) ӘР 
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Tablica D-1. (c.d.) 


А=14. 146185 
B=2946. 530 
C=42. 336 
RMSD=0. 008 kPa 


. 023 
‚021 
‚021 
. 881 
. 881 
. 881 
. 293 
‚292 
«291 
. 403 
. 403 
.401 
.222 
.220 
. 223 
. 989 
. 961 
. 962. 
. 626 
. 629 
. 629 
. 460 
. 464 
. 462. 
. 425 
. 418 
. 422 
. 791 
. 793 
. 794 
. 608 -. 
. 620 -. 
. 606 -. 


. 002 
. 001 
. 001 
. 002 
. 001 
. 001 
. 001 
. 000 
. 000 
. 001 
. 000 
. 002 
. 002 
. 004 
. 003 
. 004 
. 005 
. 004 
. 009 
. 008 
. 008 
.012 
.015 
. 013 
.014 
. 008 
.011 
. 003 
.004 
. 006 


016 
011 
015 


А=14. 000250 
=2859. 963 
С=47. 545 

RMSD=0. 007 kPa 


A=13. 883063 
B=2785.780 
C=52. 537 
RMSD=0. 008 kPa 
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ZALACZNIK E 


Tablica E-1. 


Preznosci pary w funkcji temperatury dla poszczególnych próbek. Kolumna 
trzecia zawiera odchylenia od równania Antoine'a. Punkty oznaczone 
gwiazdką zostały usunięte z korelacji. 


T/K P/kPa dP/kPa T/K P/kPa dP/kPa T/K P/kPa dP/kPa 
ә | sm NN NN 


293.223 2.351 .001 | 293.171 2.386 .003 | 293.194 2.400 .001 
2. 293.170 2.386 .003 | 293.194 2.399 .000 
2. 293.170 2.383 .000 | 293.193 2.399 .001 
2. 298.196 3.232 -.002 | 298.216 3.253 -.003 
3. 298.196 3.231 -.003 | 298.216 3.257 .001 
3. 298.200 3.233 -.002 | 298.216 3.253 -.003 
3. 303.177 4.326 -.001 | 303.194 4.356 -.003 
4. 303.177 4.327 -.000 | 303.194 4.361 .002 
4. 303.177 4.325 -.002 | 303.194 4.358 -.000 
4. 308.200 5.740 -.002 | 308.216 5.787 .002 
5. 308.201 5.742 -.000 | 308.216 5.783 -.003 
5. 308.200 5.742 -.000 | 308.216 5.783 -.003 
5. 310.213 6.391 -.021 | 310.228 6.460 -.002 
7. 310.211 6.410 -.002 | 310.228 6.459 -.002 
7. 310.211 6.413 .002 | 310.229 6.462 .001 
s 310.211 6.413 .002 | 311.229 6.828 .006 
9. 311.212 6.771 .002 | 311.229 6.828 .006 
9, 311.212 6.771 .002 | 311.228 6.828 .006 
9. 311.211 6.770 .001 | 312.230 7.208 .007 

12. 312.213 7.142 -.002 | 312.230 7.203 .002 
12. 312.213 7.147 .002 | 312.230 7.203 .002 
12. 312.213 7.147 .003 | 313.202 7.585 .000 

7.528 .003 7. 

7.530 .004 7. 

7.528 .003 8. 

8.364 -.005 8. 

8.368 -.001 8. 

8.370 .000 9. 

9.775 -.011 9. 

9.786 .000 9. 

9. 
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Tablica E-1. (c.d.) 


mw0551 mw0801 mw0897 


O O O @ 0 чычым 0 0 9 9 9 сиир P P > O W O N N N 
v$00000ududududoooooouuunhehhODOONNNN 


2: 
2. 
2. 
2. 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
5: 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
T: 
7. 
T^ 
T: 
T. 
7. 
8. 
8. 
8. 
9. 
9. 
9; 
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Tablica E-1. (c.d.) 


mw1104 mw1417 mw1772 


оосоеоволмлаллалтататтолиллврвосомм N 


2. 
2. 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
5. 
S. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
d^ 
Т» 
T. 
7. 
7. 
7. 
8. 
8. 
8. 
9. 
9. 
9: 
12. 
12. 
12. 


анионов 
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Tablica E-1. (c.d.) 


mw1933 mw2139 mw2163 


2. 
2% 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
7. 
7. 
i: 
7T. 
7. 
7. 
8. 
8. 
8. 
9. 
9. 
9. 


OO00WMNNNNNNONONADUNURAROÓWONNK 
ооофов ааа 3 9 0 0 9 О лл @ P P P O O O N N N 


=> ә н 
горо Ро 
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Tablica E-1. (c.d.) 


mw2334 mw2374 mw2399 


293.180 2.413 .010*| 293.205 2.420 .004 | 203.168 2.396 .001 
293.179 2.403 .001 2. 293.166 2.393 -.002 
293.178 2.403 .001 2. 293.165 2.395 -.000 
298.204 3.256 -.001 3. 298.190 3.249 .001 
298.203 3.259 .002 3. 298.190 3.249 -.000 
298.203 3.257 .000 3. 298.190 3.249 .000 
303.185 4.354 -.001 4. 303. 171 4.347 .001 
303.185 4.355 -.000 4. 303. 171 4.347 .001 
303.184 4.354 -.002 4. 303.170 4.346 .001 
306.199 5.165 -.001 5. 306.184 5.154 -.001 
306.199 5.166 -.001 5. 306.184 5.154 -.001 
306.199 5.166 -.001 5. 306.184 5.156 .001 
308.207 5.773 -.003 6. 308. 191 5.761 -.003 
308.207 5.777 .001 6. 308.192 5.761 -.003 
308.208 5.778 .002 6. 308. 191 5.763 -.001 
309.212 6.108 .004 6. 309.170 6.084 .001 
309.213 6.102 -.002 6. 309. 171 6.084 .000 
309.213 6.104 -.000 6. 309. 171 6.084 .001 
310.220 6.448 -.001 6. 310.205 6.440 .003 
310.220 6.448 -.001 7. 310.204 6.438 .002 
310.220 6.448 -. 001 T. 310.205 6.438 .002 
310.739 6.633 .000 T, 310.684 6.606 -.000 
310.739 6.632 -.001 T; 310.684 6.606 -.000 
310.739 6.633 .000 T. 310.684 6.607 .000 
311.216 6.808 .001 LE 311.195 6.790 -.002 
311.216 6.808 .001 8. 311.195 6.792 -.000 
311.215 6.807 .001 8. 311.195 6.792 .000 
311.376 6.866 .000 8. 311.358 6.853 .001 
311.376 6.864 -.001 9. 311.358 6.855 .003 
311.376 6.864 -.001 9. 311.358 6.854 .002 
311.749 7.005 .000 9. 311.624 6.950 -.001 
311.749 7.006 .001 12. 311.625 6.951 .000 
311.749 7.006 .001 12. 311.625 6.951 .000 
312.206 7.178 -.001 12. 312.189 7.165 -.001 
312.206 7.181 .002 6. 312.189 7.167 .001 
312.205 7.181 .002 7. 

313.178 7.563 .001 7. 

313.178 7.563 .002 T. 

313.178 7. 560 -.001 T. 

314.184 7.978 .002 7. 

314.184 7.975 -.000 Y 

314.184 7.973 -.002 7. 

315.187 8.407 -.001 8. 

315.187 8. 409 .000 8. 

315.187 8. 407 -. 001 8. 

318.198 9.831 .001 9. 

318.198 9.831 .001 9. 

318.198 9.830 .000 9. 

323.180 12.638 -.002 

323.180 12.640 .000 


mà 
N 
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Tablica E-1. (c.d.) 


mw2556 mw2656 mw2823 


2. 
Z: 
2. 
3. 
3. 
3: 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
7. 
7. 
7. 
7. 
7. 
T« 
8. 
8. 
8. 
9. 
9. 
9: 


POOOWWNNTANNNOADODODUMURARCÓGONNN 
OOOBWOWNNANNANNONONOUNUUWNARACECGNNN 
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Tablica Е-1. (c.d.) 


mw2987 mw3146 mw3322 


$0090puulddudoooooounuumRRRODONNN 
000000 4 O O O O O. O (Л (Л (Л i> i> > Q Q Q N N ƏN 


к> кэ къ кэ кэ xn рз рз рэ раз раз рз рз рз рз IH рз үз рз рз рз | p ра ка ра рэ ра рз 
0000 -1 -1-Ј6ОСОС 0 Ф AO 0 O ‹о O о хо м) -Ј) N (л (л (л м: N NAMM 
NNNNUURA AN eNOOOOOCOONNNAANANAN Aeree 


2. 
2. 
2. 
3: 
3. 
3. 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
ТА 
7. 
Ta 
7. 
7. 
7. 
8. 
8. 
8. 
9. 
9. 
9. 
12. 
12. 
12. 
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Tablica E-1. (c.d.) 


2: 
2. 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
T. 
T. 
7. 
Ts 
7. 
7. 
8. 
8. 
8. 
9. 
9. 
9: 
12. 
12. 
12. 
2. 
2. 
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Tablica E-1. (c.d.) 


= N N N = = = 
@ O O O ч JJ 
Q => O = O O MJ 


оо O % @ @ ааа 3 3 3 0 0\ 9 9 9 9 (л (лл e > г о QO O N N N 


2. 
2. 
2. 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
S. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
Т» 
Те 
T5 
7. 
7. 
7. 
8. 
8. 
8. 
9. 
9. 
9. 


M - 
оодофеовеоллаллалеттгтттиллрввосомм N N 


eeNNNNN = N N N N N N N N N N N N 
N О = = = ОООМОО ВБ нь [ук = к ҥ к кю Ооо О 
Ф Y.N Q Q P к = N N Q O O O со со ОДУ 502» 
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Tablica E-1. (c.d.) 


ооо — 3 3 0 0 9 о\ тоо D DNNN NNN ооо @ Q N N N 


2. 
2. 
2. 
2. 
3. 
3. 
3. 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
7. 
7. 
7. 
8. 
8. 
8. 
9: 
9. 
9. 


OOOWOONNJDDNNDDNAADUUUARAGOWEONNNNN 
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Tablica E-1. (c.d.) 


293.174 2.119 -.002 | 293.174 1.836 -.011 | 293.174 .193 .003 
2. 293.175 1.837 -.010 
2; 293.174 1.838 -.009 
2. 298.198 2.463 .010 
2. 298.198 2.463 .009 
2. 298.198 2.461 .008 
2. 303. 181 3.217 .005 
3. 303. 181 3.214 .002 
3. 303. 181 3.219 .007 
3. 308.170 4.156 -.003 
4. 308.174 4.159 -.001 
4. 308.176 4.162 .002 
4. 310.215 4.610 .001 
5. 310.215 4.609 -.001 
5. 310.215 4.609 -.001 
5. 311.216 4.848 .004 
5. 311.216 4.845 .001 
S. 311.216 4.842 -.002 ç 
S. 312.217 5.087 -.001 1. 
5. 312.217 5.087 -.001 1. 
5. 312.217 5.089 .001 1. 
5. 313. 191 5. 334 -.002 1. 
6. 313.191 5.332 -.004 1. 
6. 313.190 5.335 -.001 15 
6. 315. 201 5.877 -.002 1. 
6. 315.201 5.878 -.001 1. 
6. 315.200 5.877 -.002 1. 
6. 318.209 6.775 -.002 
7. 318.208 6.775 -.002 
7. 318.209 6.774 -.003 
7. 323. 191 8.514 .003 
8. 
8. 
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Tablica E-1. (c.d.) 


| Ко с 6 КУ КУ гок ою кю кю ә 
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ZALACZNIK F 


Tablica F-1. podane na nastepnej i kolejnych stronach zawiera zes- 
tawienie najważniejszych układów dwuskładnikowych badanych pod kątem 
zjawisk w otoczeniu punktu krytycznego mieszalności. Przedstawione tutaj 
informacje nie stanowią kompletu bibliografii dla podanych układów do 
początku lat 90-tych, lecz stanowi obraz bogactwa zagadnień badanych w 
tej dziedzinie. W kolumnie 3 podane są temperatury krytyczne, które jak 
widać, stanowią raczej wąski zakres stwierdzonych eksperymentalnie CST 
dla dużo więcej układów dwuskładnikowych (Francis, 1963); wynika to z 
wymogów badań CST podanych we wstępie tej pracy. W kolumnie 4 podano 
identyfikatory własności badanych w okolicy tego punktu. Znaczenie ich 
jest następujące: 

cpex - ciepło właściwe pod stałym ciśnieniem; 

CXC - krzywa współistnienia i związane z nia wykładniki a i В; 

diff - współczynnik dyfuzji; 

elect - własności elektryczne t. j. МОЕ czy efekt Kerra; 


expan - oznaczanie współczynnika rozszerzalności termicznej; 


gamm - rozpraszanie i dyfrakcja promieniowania gamma; 
hex - entalpia mieszania; 
nmr - badania za pomocą JRM; 


nucl - badania krzywej spinodalne jj; 
refr - współczynnik refrakcji; 

scat - badania rozpraszania światła; 
sound - pochłanianie ultradźwięków; 


visc - lepkość; 


vex - nadmiarowa objętość mieszania, objętość molowa; 
vle - pomiary prężności pary, ewaluacja g 
wett - zjawisko anomalnej zwilżalności; 


W kolumnie 5 podane sa numery odno$ników, których lista podana jest bez- 
pośrednio za tablicą F-1. Niektóre pozycje pokrywają sie z wymienionymi 
w spisie literatury lecz ten system pozwala utrzymać niniejszy załącznik 


jako integralną całość. 
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Tablica F-1. 


1 
тео тт —— WT 
Bp |» pe 


1| Woda + 3-metoksypropionitryl 311. 143 


| 
+ eter izobutylowy glikolu 298.95 
etylenowego 


+ kwas izobutylowy (masłowy) 299.536 
318.054 


+ tetrahydrofuran 


к ЖЕ = 


| + trojetyloamina 


E + М, N-dwumetylobutyloamina 


10| anilina + heksan 330. 05 CXC 


Tablica F-1. (c.d.) 


12| metanol * cykloheksan 


318.624| cxc 1, 33,41,50, 
87,96 
20,21 
ем 33,50 
13 + heptan 323.88 срех 46 
схс 88 


| 


15| etanol + heksadekan 325. 25 ћех 113 
уех 113 
у1е 113 


nitrometan + dwusiarczek wegla 335. 133| еже |28, 32, 102 


nitroetan + heksan CXC 68 
hex 89 
36 


| + 3-metylopentan 


+ 2,2,3-trojmetylopentan |303. 18 
= б 
Sound |62 


| 23| nitrobenzen + репїап 297.16 [vise |во | 
+ heksan 


293. 11 


+ cykloheksan 
+ heptan 292.18 | gamm |2 
Visc 18,80 


+ | + 2,2,3-trojnetylopentan 2, 3= | + 2,2,3-trojnetylopentan 302.31 | sound |43 | 
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Tablica F-1. (c.d.) 


28| nitrotoluen + alkany 


29| 2-nitrofenol + alkany 


30| 3-metoksypropionitryl 


* weglowodory C5 


31| N-metylopirolidon + cyklooktan 


N-metylopirolidon * nonan 117 


32 xc 
3-fenylopropanol-1 + heksan elect |60,75 


33 


34| alkohol benzylowy 


+ metylocykloheksan elect |76 


alkohol decylowy + acetonitryl | elect |76 
+ nitroetan 76 


+ fenyloacetonitryl |287.20 


35 


36 


C 


Q 


+ ftalan dwumetylowy 
+ dwuacetat etylenu|291.6 
+ rozpuszczalniki polarne 


41| acetonitryl + dwusiarczek wegla 


心 
~ 


+ cykloheksan 


43| 2-metoksyetanol + alkany 


anilina * cykloheksan 


chlorowcopochodne weglowodorów 
* weglowodory 


perfluorometan + metan 


Tablica F-1. (с.а.) 


=== perfluoroheksan + heksan 
| 50 perf luoroheptan + heptan SEE 


+ czterochlorek wegla|330.84 


52| perfluorometylocykloheksan 
+ metylocykloheksan 


+ alkany 


perfluorobenzonitryl + cykloheksan 一 | 
+ 2,2,3-trojmetylopentan|291. "=E — 
vex 


57| perfluorocyclooxide 
+ czterochlorek wegla 


Е polistyren + cykloheksan 
62 * trans-dekalina 108 
| 108 
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